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INTRODUCTION 


The  purpose  of  this  report  is  to  discuss  a user- 
oriented  computer  routine*  to  analyze  wire-type  antennas 
directly  in  the  time  domain.  Instructions  for  program 
input  data  and  Interpretation  of  the  output  are  presented 
in  concise  form  to  permit  use  of  the  program  with  minimum 
theoretical  understanding.  Sample  cases  are  presented 
which  v;alk  the  reader  through  a typical  problem  solution. 

The  computed  results  are  critically  compared  with  published 
experimental  results  to  establish  the  validity  of  the  computer 
routine . 

The  time  domain  approach  is  the  most  direct  way  to 
obtain  the  response  of  an  antenna  to  a transient  excitation. 
The  output  of  the  time  domain  solution  is  a time  history  of 
the  current  distribution  on  the  antenna  caused  by  a transient 
voltage  or  field  excitation.  Integration  of  this  current 
provides  the  time  varying  far  field  or  the  receiving  response 
of  the  antenna.  By  using  Fourier  transform  the  results  can 


*The  comouter  routine  was  developed  by  Associates  and  has 
been  implemented  on  corrmuters  at  RADC,  NRL,  and  'IAVSFC  under 
contract  N0001 4-73-C-0099  and  RADC  F30602-72X0008. 
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be  expressed  in  the  frequency  domain  as  the  wide  band 
response  of  the  antenna.  A single  time  domain  solution 
yields  data  over  the  bandwidth  associated  with  the  time 
waveform  of  the  excitation. 

This  report  is  divided  into  three  parts.  First  a 
reneral  description  of  the  formulation  is  provided.  Then 
program  input  and  instructions  for  usinp,  the  computer 
routines  are  piven,  includinp  preparation  of  input  data 
and  interpretation  of  output  for  a snecific  example. 

Finally  the  output  is  discussed  and  critically  connared 
with  experimental  results.  A listing  of  the  Fortran  source 
nro<-ran  is  included  as  Appendix  A. 


2 


SECTION  I Formulation  of  the  Time  Domain  Equations 

1.0  For  the  sake  of  completeness  and  providing  a 
self-contained  report,  this  section  presents  without 
rigorous  lerivation  the  integral  equations  which  were 
developed  and  used  in  the  Tine  Domain  Program  by  MS 
Associates  under  government  contract-*-. 

This  computer  routine  is  based  upon  the  time  domain 
equivalent  of  Packlington ' s Integral  equation  which  is 
widely  used  for  frequency  domain  numerical  solutions. 

In  the  time  domain  the  equation  becomes  an  integro- 
differential  equation  and  the  unknown  current  is  a 
function  of  two  variables,  position  and  time.  The  time 
domain  equation  can  be  derived  by  using  Maxwell's  time 
varying  equations  and  the  representation  of  the  electric 
field  in  terms  of  the  vector  and  scalar  potentials  and 
Introduction  of  the  retarded  potential.  Consider  the 
general  antenna  structure  depicted  in  Figure  1.1. 

The  electric  field  at  an  arbitrary  point  P (r,  e,  $) 
outside  the  volume  current  and  charged  source  V'  is 
expressed  as 
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E (r,t) 


v * (r,t) 


d 
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(r , t) 


where  A(r,t)  = u J(r',t-  R)  dv' 

4 ir  / C 
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(r,t)  = 


_ 1 


4it  f 


(r  , t - 
£L  dv' 


(3) 


V' 


are  the  retarded  vector  and  scalar  potentials  at  point 
P(r,  0,  0)  respectively  and 

J (rj  t ' ) = current  density  in  Volume  V' 

p (r',  t')  = charge  density  in  volume  V' 

C = velocity  of  light 

t'  = t-  5 = retarded  time 

C 

In  the  above  expressions,  the  primed  quantities  represent 
source  coordinates  and  the  unprimed  represent  the  field 
point. 

If  one  considers  a thin  wire  structure  the  current 
and  charge  are  now  confined  to  a filamentary  path  on  the 
wire  axis  and  the  integration  represented  in  equations 
(2)  and  (3)  reduce  to  contour  integrals  along  the  wire 
axis . 
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Therefore,  the  equations  representing  the  potentials  as 
a function  of  current  and  chnrpe  distributions  can  be 
rewritten  as 


A (r, t) 

= U r IMS', 

41.  * R 

±L± 

ds  ' 

(4) 

♦ (r,t) 

= 1 r q (S'  , 

4tt  c f i R 

JL 1 

ds  ' 

(5) 

3 I 

3S  ' 

(s\t*)  = - 3 

3 1 ‘ 

q (S' 

f) 

(6) 

where  t'  = t - R/C  retarded  time  and  *■  is  the  contour  of 
the  wire,  s’  is  the  source  position  vector  at  point  s' 
(3ee  Ri^ure  1.1b).  I (s',  t')  is  the  unknown  current 
distribution  which  is  a function  of  position  and  retarded 
time  and  q is  the  linear  charge  density. 

Performing  the  operations  on  the  potentials  as  re- 
quired in  eq  (1),  the  electric  field  can  now  be  expressed 
in  retarded  time  t'  as: 
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E (r ,t)  = - 


3 

Zt1 


I (s' ,t' ) + 


CR 

7 


_3  I(s',t’) 

3s' 


+C 


2 


3 

3s  ' 


/fcl  I (s' ,t' )dt' 
o 


ds' 


(7) 


where  s'  is  unit  vector  tangent  to  wire  at  s'. 


If  now,  R -*■  ® and  only  those  terms  are  retained  which 

decrease  as  1/R,  equation  (7)  reduces  to  the  farfield 
expression 


E (r , t)  = - us1 
4ir 


r I 3 

1 r TtT 


i (s' ,t' ) + c R 


_3  I(s’ ,t')]  ds'  (8) 
3s’ 
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Since  R » R can  be  replaced  by  RQ,  the  vector  from  the 

point  of  observation  to  the  origin  (See  Figure  lc) , and 
using 

3 I(s',t')  = d I(s',t')  - R°  ' S JL  Ks',t')  (9) 

3s'  ds7"  c 3t' 


The  far  field  formulation  can  be  written  as 

E (r,t)  = - Ml_  / 3 I(s\t')  [S'  - (S'.  Rq)  P0]  ds'  (10) 

4"Ro  ■ l*  TF 


Equation  10  is  the  final  form  of  the  lntegro  differential 
equation  for  the  radiated  field.  It  is  obvious  that  if 
the  current  on  the  antenna  is  known  both  as  a function  of 
position  and  time,  it  is  a relatively  simple  matter  to 
compute  the  field  or  received  response.  The  crux  or  the 
computer  routine  herein  is  the  computation  of  the  induced 
currents.  :7hile  no  attempt  will  be  made  in  this  report 
to  develop  the  numerical  techniques  employed,  a general 
outline  of  the  concept  follows. 

The  method  of  moments(l)  is  used  to  reduce  equation  (7) 
(with  appropriate  boundary  conditions)  to  a set  of  linear 
equations  to  solve  for  the  unknown  current,  involving 
both  position  and  time  rather  than  the  single  variable 
positions  as  occurs  in  the  frequency  domain,  fubsectional 
collocation  methods  are  applied  to  solve  the  linear 
equations  for  current  coefficients.  This  requires  that 
the  antenna  structure  be  divided  into  a number  of  seg- 
ments and  time  be  divided  into  equal  increments.  The 
details  of  the  numerical  techniques  will  not  be  presented 
here  because  of  their  complexity.  However,  the  effects 
of  the  numerical  techniques  are  to  compute  the  current 
(position  and  time)  at  the  center  of  each  subsection 
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due  to  impressed  excitation.  Once  these  currents  are 
computed  the  total  field  is  determined  through  equations 
( 8 ) throuph  ( 10) . 
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SECTION  II  Sample  Case  Solution 

2.0  The  preceding  section  presented  the  development 
of  the  integro  differential  equations  used  in  the  solution 
of  the  time  domain  response  of  antennas,  while  the 
development  of  the  expressions  is  rather  sketchy,  they 
suffice  to  outline  the  concepts  and  mathematical  procedures 
employed.  It  is  felt  that  from  a user  point  of  viev.> 

Section  I provides  sufficient  development  of  the  time 
domain  equations  with  minimum  theoretical  unders tandinp . 

In  this  section,  the  computer  routine  will  be  applied  to 
a simple  antenna  geometry  to  emphasize  user  procedures  and 
data  input  preparation. 

2.1  Description  of  Problem. 

The  problem  to  be  addressed  in  this  section  consists 
of  a cylindrical  mononole  antenna  over  a perfect  ground 
plane.  This  antenna  was  addressed  both  theoreticallv  and 
experimentally  by  Schmitt,  et.  al.,2  and  provides  n mean' 
of  comparison  against  commuted  performance.  The  actual 
"■eometry  is  depicted  in  Eigure  2.1.  "'he  monopole  is  ren- 
resented  by  a 3 meter  wire  with  length  to  radius  ratio 
h/a  = In  Schmitt's  experiments,  pulse  lengths  were 

used  which  varied  from  0.2  to  2.0  h/c  -where  h is  the 
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GENERAL  FIELD  POINT 


kv  kilkuLL 


length  of  the  wire  and  c Is  the  velocity  of  light.  The 
rise  time  for  the  pulses  Is  0.05  h/c.  The  excitations 
considered  are  voltage  sources  to  obtain  the  radiated 
field  response  and  field  excitations  to  examine  receiving 
time  domain  performance,  ^he  excitations  are  short  base 
band  time  waveform  voltage  and  field  excitations.  mhe 
time  characteristics  of  the  excitations  are  3 db  width  of 
0.16  h/c  and  2.0  h/c  with  rise  time  on  the  order  of  0.05 
h/c . 

2.2  Data  Input. 

In  order  to  Input  this  problem  Into  the  Time  Domain 
(TIMDOM)  program  the  problem  must  be  cast  into  input  data 
format.  The  input  data  can  be  separated  Into  six 
categories . 


1. 

Comments 

4. 

Loadinr 

o 

Data 

5. 

Excitation 

Y. 

Geometry 

6. 

Output  Data  requests 

The  data  is  input  as  nunched  cards  and  will  be 
described  in  the  following  paragraphs  as  applied  to  the 
example  problem. 

2.2.1  Comment  Cards. 

The  data  card  deck  must  begin  with  one  or  more 
comment  cards  which  nay  contain  a description  of  the 
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2.2.2  First  Data  Cards. 


This  card  does  not  contain  a mnemonic  identifier. 

The  contents  of  the  card  are  divided  into  7 usable  fields 
from  column  1 through  45.  The  parameters  included  on  the 
card  are:  Time  increment  for  solution  (TT);  Maximum  time 

for  which  currents  are  computed  (TMAX);  Number  of  time 
steps  for  solution  (NS);  Symmetry  planes  if  any;  Structure 
loading  if  any  (SL);  and  type  of  loads.  See  card  below 
for  specific  example. 

Fields  1-10  are  blank,  thus  not  specifying  TI.  How- 
ever, its  value  is  computed  automatically  by  making  it 
slightly  larger  than  the  maximum  segment  length/c.  The 
integer  number  occupying  field  11  - 20  specifies  TMAX. 
Either  TMAX  or  NS  must  be  input.  NS  which  occupies  field 
21  - 25  is  left  blank  and  is  therefore  computed  automati- 
cally. Columns  26  and  27  are  unused  on  this  data  card. 
Columns  28,  29  and  30  specify  the  planes  of  symmetry.  For 
the  case  being  considered  the  E in  column  30  specifies  an 
electric  conducting  plane  located  at  the  Z = 0 plane. 

Fields  31  - 35  identify  structure  loading,  in  this  case 
the  one  (1)  signifies  impedance  loading  on  the  structure. 
Fields  36  - 40  and  4l  - 45  are  used  to  identify  and  tag 
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those  segments  of  the  wire  antenna  for  which  currents  are 
to  be  printed.  In  this  case  where  they  are  left  blank, 
currents  will  be  printed  for  all  segments  as  a default 
condition . 
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2.2.3  Geometry  Cards. 

The  purpose  of  the  geometry  cards  is  to  input  the  antenna 
structure  geometry.  Specifically  these  cards  (mnemonic) 

GW  generate  lengths  of  straight  wire  between  two  points  in 
space  and  divide  them  into  an  equal  number  of  segments. 
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The  card  is  divided  into  10  fields  which  are  used  to 
identify  the  end  points  of  the  wire  (in  x,  y,  z space), 
identify  the  number  of  segments,  and  provide  the  radius 
of  the  wire.  Any  number  of  wires  can  be  generated  in 
this  fashion,  however,  only  a single  wire  can  be 
generated  on  one  card. 

See  specific  example  below  for  detailed  description 
of  the  data  fields.  Fields  3-5  are  used  to  identify,  by 
tag  number,  a segment  of  the  wire  for  later  use  if  needed. 
In  this  case,  no  segment  is  tagged  since  the  antenna  is 
composed  of  a single  wire.  Fields  6-10  are  used  to 
specify  the  number  of  segments  into  which  the  wire  is 
to  be  divided,  in  this  case  32.  Three  fields  are  dis- 
tributed between  columns  11  and  43  which  are  used  to 
specify  the  x,  y,  z coordinates  of  one  end  of  the  '"ire. 

In  this  case  the  wire  end  is  positioned  at  the  origin. 
Columns  4l  - 70  are  used  to  specify  the  other  end  of  the 
vrl  re.  The  data  in  column  64  places  the  wire  3 meters 
plonr  the  Z axis.  See  Figure  2.1.  Fields  71  - 80  are  used 
to  specify  the  radius  of  the  wire  in  this  case  0.00332 
meters.  Since  this  case  only  involves  a sinrle  wire 
antenna  only  one  Geometry  card  is  used.  If  a more 
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complex  structure  Is  being  modeled,  as  many  cards  as  re- 
quired would  be  input.  If  the  units  input  are  other  than 
meters,  then  the  units  must  be  converted  to  meters  through 
the  use  of  a structure  scaling  card  (OS).  However,  none 
is  needed  in  this  case  deck. 

The  geometry  input  data  must  end  with  a geometry  end 
card.  This  is  represented  by  mnemonic  OK  and  all  fields 
are  blank. 
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2.2.4  Impedance  Loading. 

If  an  antenna  structure  is  impedance  loaded,  allowance 
must  be  made  for  inputting  the  data  into  the  program.  This 
is  accomplished  by  a group  of  impedance  loading  cards. 

This  card  is  divided  into  six  usable  fields  which  ore 
used  to  identify  the  wire  segments  to  be  loaded  and  specify 
the  complex  impedance  of  the  load.  See  card  below  for 
format  details.  Fields  1-5  are  used  to  specify  the 
number  of  impedance  cards.  If  non  zero,  the  program 
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In  this  case  that 


searches  for  succeeding  loadinp  cards, 
field  is  zero  inplyinp  this  is  the  last  Impedance  loadin 
card.  Fields  6-10  and  11  - 15  are  used  to  identify,  by 
tap  nunber,  the  sepment  to  be  loaded.  In  this  case  sepment 
one  (1)  is  identified.  Fields  21  - 30,  31  - 40,  41  - 50 
are  used  to  specify  in  ohms  the  resistance,  inductance, 
and  capacitance  of  the  impedance.  In  this  case,  a 50  ohm 
pure  resistance  is  input  to  simulate  the  coax  feed  point 
impedance . 
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2.2.5  Excitation  Card. 

Having  described  the  geometry  and  loads,  the  excitation 
must  be  specified.  This  information  is  input  through  the 
excitation  card.  It  is  divided  into  8 fields  which  are  used 
to  specify  type  of  excitation,  i.e.  voltage  or  field  excita- 
tion; number  of  sources;  3 db  pulse  width  of  the 
excitation;  incidence  angle  for  field  excitation.  A 
discussion  of  the  various  fields  references  the  excitation 
card  of  the  example  case. 

Fields  1-5  are  used  to  identify-  the  time  variation  of 
the  excitation.  If  non  zero,  the  time  variation  is 
arbitrary  and  the  value  appearing  in  this  field  specifies 
the  number  of  values  to  be  read  from  succeeding  Source 
Function  cards.  In  this  case  the  value  zero  in  the  field 
specifies  Gaussian  time  variation.  The  fields  in  columns  5 - 
10  specify  the  number  of  sources.  In  this  case  one  voltage 
source  is  specified.  The  value  in  Fields  11  - 20  specifies 
the  3 db  pulse  width  of  the  excitation,  i.e.  1.6  ns  which 
is  equivalent  to  0.16  h/c  for  this  case.  Three  fields  in 
columns  21  - 30,  31  - 40,  41  - 50  are  used  to  specify  the 
0,  0,  n parameters  of  the  incidence  angle  for  field 
excitation.  The  field  source  position  is  specified  by 
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9,  and  0 (convention  spherical  coordinates),  and  n is 
the  electric  field  polarization  which  is  the  angle  in 
degrees  between  the  9 unit  vector  and  the  direction  of 
the  electric  field.  For  a voltage  source  excitation  as 
is  the  case  in  this  example  these  fields  have  no  meaning. 
Fields  51  - 60  are  used  to  reference  the  incident  field 
relative  to  the  antenna  structure.  It  has  no  meaning  and 
is  left  blank  in  this  case.  The  last  field  on  the  card 
specifies  the  time  at  which  the  Gaussian  pulse  reaches  Its 
maximum.  It  Is  chosen  such  that  the  pulse  has  a small 
percent  of  its  peak  value  at  t = 0.  If  left  blank,  the 
program  computes  an  appropriate  value  a 3%  of  peak  at 
t = 0).  In  this  case  it  was  specified.  This  card 
completely  specifies  all  the  parameters  required  to  input 
the  excitation. 
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2.2.6  Source  Voltage. 

This  card  Is  used  to  specify  the  peak  voltages  of  the 
individual  sources  used.  It  consists  of  three  usable 
fields  which  identify  the  segments  on  which  sources  will 
be  located,  and  input  the  peak  voltage  of  the  sources. 

On  this  card  (see  below)  of  the  input  data  deck,  fields 
1-5  and  6-10  are  used  to  specify  the  tag  number  of 
the  source  segment  and  the  segment  number  of  a set  of 
segments,  respectively.  In  this  case  since  there  is  only 
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one  tag  number  (only  one  wire)  the  tag  field  is  left  blank 
and  the  segment  number  1 is  specified.  That  is,  the  source 
is  located  at  the  base  segment  of  the  32  segment  wire 
monopole.  The  peak  value  of  the  source  is  unity  as 
specified  in  Fields  11  - 20. 

1. 
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cards  to  be  discussed  pertains  to  requesting  output.  This 
set  is  referred  to  as  Data  Request  Cards. 

2.3  Data  Request  Cards. 

A set  of  commands  is  required  to  identify  and  extract 
the  proper  data  from  the  output  solution.  For  example, 
the  program  computes  induced  current,  impedance,  admittance, 
field  patterns,  receive  response,  etc.  The  data  request 
cards  are  used  to  call  out  and  print  all  or  portions  of 
this  data.  Specific  requests  as  pertain  to  the  example 
run  in  this  report  will  be  addressed. 

2.3.1  Energy  Budget. 

This  card  contains  only  the  mnemonic  EB.  Its 
purpose  is  to  cause  the  computations  and  printing  of  the 
total  energy  input,  energy  lost--  in  the  loads  if  any  and 
the  total  energy  radiated.  It  results  in  the  connutation 
of  time  domain  efficiency  number.  These  nuantities  will 
be  discussed  in  the  output  data  section  of  the  report. 
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2.3.2  Radiated  Field. 

This  command  requests  the  computation  of  the  radiated 
field  and  outputs  appropriate  data  in  the  time  domain. 
Further,  this  card  is  used  to  transform  the  time  domain 
output  into  the  frequency  domains  if  desired.  The  card 
format  will  be  discussed  with  reference  to  the  card  in  the 
sample  case  data  deck  below.  The  format  of  the  card  is 


divided  into  9 usable  fields  including  the  mnemonic  RF 
Fields  3-5  are  used  to  control  the  computation  and 
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printing  of  the  Fourier  Transform  and  punch  cards  if 
desired.  In  this  case  0 means  no  transform  required. 

Fields  6-10  are  used  to  control  punching  cards  for  electric 
field  output.  In  this  case  the  electric  field  will  be 
printed  directly  and  not  punched  on  cards.  The  remaining 
six  fields  specify  the  direction  in  space  for  which  the 
field  is  to  be  computed.  Specifically,  fields  11-15 
and  16  - 20  specify  that  the  field  vrill  be  computed  for 
one  0 and  one  0 direction.  Fields  21  - 30  and  31  - ^0 
specify  these  directions  as  90°  and  0°  respectively,  i.e. 
alon>r  the  ground  plane.  Fi°lds  !il  ••  50  and  51  - 5o  specify 
the  h and  0 stepping  increment  for  field  computation. 

In  this  case  both  are  zero  since  only.one  direction  is  of 
interest . 
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2.4  The  preceding  paragraphs  describe  the  preparation 
of  a data  input  deck  for  a specific  case.  Recall  the 
problem  was  to  determine  the  radiated  field  of  a nonopole 
antenna  on  a perfect  conducting  ground  plane  fed  by  a 

50  ohm  coaxial  cable  (see  Figure  2.1).  The  sample  input 
data  deck  was  presented  to  input  the  problem  parameters , 
i.e.  geometry,  excitation,  loading,  etc.  and  request 
current  computation  and  far  field  transmitting  response. 

The  following  paragraphs  will  present  and  discuss  the 
output  of  the  samnle  case.  The  radiated  and  received 
responses  will  be  presented  and  compared  to  experimental 
results  in  Section  III. 

2.5  Program  Data  Output. 

The  purpose  of  this  portion  is  to  familarize  the 
reader  with  the  raw  output  of  this  program.  The  output 
presented  is  a result  of  the  input  problem  discussed  in 
paragraph  2.1.  Although  the  output  of  the  program  is 
automatically  identified  by  labels,  each  portion  of  the 
sample  output  has  been  identified  by  circled  numbers  for 
easy  reference  by  the  reader.  The  sample  outDut  which 
is  discussed  below  is  located  at  the  end  of  this 
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section . 


2.5.1  Output  Discussion. 

The  program  title  Time  Domain  Antenna  Modeling 
Program  appears  at  (l) automatically  with  each  output  as  the 
first  heading.  The  comment  cards  which  were  presented  as 
input  data  follow  the  heading  at  ■'2  and  can  contain  as  many 
statements  as  desired.  At  3 the  first  data  card  is  read 
and  printed.  This  card  contains  various  time  stepping 
information  as  well  as  any  problem  symmetry  (see  input 
data  card  discussion  paragraph  2.2.2). 

Next  the  geometry  input  cards  are  printed  at  4 . In 
this  case,  a single  monopole  is  modeled  by  32  segments 
with  wire  radius  0.0032  meters.  These  outputs  identify 
the  wires  by  tag  numbers.  The  parameters  specified  by 
the  geometry  cards  are  used  to  determine  the  segment 
centers,  lengths,  and  orientation  (a»  8 orientation 
angles)  which  are  printed  in  the  output  at  5.  This  print- 
out also  contains  segment  interconnection  data  and  tag 
numbers.  The  I-  1+  notation  is  used  to  reference  current 
flow,  i.e.  current  flows  from  negative  to  positive. 

The  wire  information  is  summarized  at  6. 

The  table  automatically  printed  out  at  7 summarizes 
the  time  stepping  Information  used  in  the  solution.  It 


30 


provides  the  number  of  time  steps  and  time  stepping  interval. 
If  any  one  of  these  quantities  was  not  input  as  data  (see 
first  data  input  card  paragraph  2.2.2),  it  is  automatically 
computed  and  printed  here. 

Any  options  such  as  symmetry  planes,  impedance 
loading,  etc.  are  printed  at  location  (B).  This  location 
also  Identifies  the  segment  for  which  current  will  be 
printed . 

The  parameters  of  the  excitation  are  provided  at  9. 

In  addition  to  printing  the  input  excitation  data  card, 
the  Gaussian  characteristic  time  information  is  provided. 

As  pointed  out  in  this  example,  the  source  has  11. 5%  of 
its  peak  value  at  t = 0.  This  value  was  specified  as 
data  input  otherwise  it  would  have  been  automatically 
commuted  and  printed  here. 

At  location  10,  all  voltage  source  information  is 
summarized.  It  specifies  the  number  of  sources, 
location  and  peak  voltages.  If  field  excitation  had  been 
used  this  section  would  be  replaced  by  angle  of 
incidence  information. 

At  this  point , the  nrogram  begins  the  solution  for 
the  currents  on  the  segments.  The  results  are  printed 
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out  as  a block  of  data  at  11.  The  currents  are  printed 
from  left  to  right  In  order  of  increasing  segment  number 
for  each  time  step.  In  this  case  there  are  160  time  steps 
of  = 0.3  ns  each.  The  value  of  the  source  strength  at 
each  time  step  is  the  final  entry  of  each  block,  labeled 
excitation . 

The  results  of  the  Energy  Budget  request  card  are 
printed  at  12.  ^he  Time  Domain  Efficiency  is  defined  as 
the  ratio  of  the  total  energy  radiated  to  the  total  energy 
input  and  accounts  for  any  energy  lost  in  loads  on  the 
structure.  If  no  loads  were  present  in  the  problem  this 
efficiency  would  be  100?. 

The  next  block  of  data  13  ih  the  output  is  the  result 
of  the  request  for  radiated  field  computations.  It  pro- 
vides the  electric  field  (volts/meter)  variation  for  each 
time  step  at  a specified  point  in  space  (theta  and  phi). 

It  also  provides  polarization  tilt  information  which  is 
also  a useful  performance  parameter. 

The  preceding  paragraphs  have  presented  a walk 
through  discussion  of  the  use  of  the  time  domain  program 
to  input  a problem  and  interpret  the  output . The  program 
was  written  to  emphasize  user-ease  with  a minimum  amount 
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of  previous  training  In  numerical  techniques  or  computer 
understanding.  In  the  next  section  the  results  of  the 
sample  problem  will  be  analyzed  and  compared  with 
experimental  results. 
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SECTION  III  Comparison  of  Computed  and  Measured  Results 

3.0  The  purpose  of  this  section  is  to  present  and 
analyze  the  response  of  a monopole  antenna  to  short  base- 
band pulses.  Both  the  transmitting  and  receiving  computed 
responses  will  be  presented  and  compared  to  measurements 
previously  published.2  The  computed  results  were 
obtained  by  plotting  the  output  of  the  TIM DOM  program 
previously  discussed. 

3.1  Transmitting  Response. 

The  cylindrical  monopole  and  its  zero  resistivity 
ground  plane  and  coaxial  feed  geometry  are  shown  in 
Figure  3.1.  The  specific  antenna  geometry  and  excitation 
used  in  the  computation  was  scaled  in  the  time  domain  to 
the  actual  measurement  geometry.2  Two  dc  pulse  excita- 
tions were  used  and  are  shown  in  Figure  3-2.  The 

excitations  are  plotted  with  respect-,  tn Jx/c.  the  one-way 

* 

travel  time  on  the  antenna  where  "h'  is  the  monopole 
length  and  "c"  is  the  free  space  velocity  of  propagation. 
The  narrow  excitation  is  a C.aussian  pulse  with  3db  pulse 
width  of  0.16  n/c  (1.6  ns)  since  h is  three  meters  in  the 
computed  case.  The  broader  excitation  represents  a 
rectangular  pulse  with  3db  width  of  2.0  h/c . 
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MONOPOLE  OVER  GROUND  PLANE 
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The  rise  time  for  the  pulses  is  approximately  0.05  h/c. 

For  the  transmitting  case  these  excitations  simulate 
voltage  sources  of  unit  amplitude  located  at  the  base 
of  the  monopole.  The  computed  response  of  the  mononole 
to  the  Gaussian  voltage  excitation  is  plotted  in 
Figure  3.3*  The  input  pulse  was  rescaled  and  is  included 
in  the  figure  for  direct  time  comparison.  The  plot  rep- 
resents the  time  history  of  the  radiated  field  in  volts/ 
meter  at  point  P (r,  ?0°,  0°),  i.e.,  at  the  ground  plane. 
It  is  evident  that  the  monopole  antenna  radiates 
essentially  replicas  of  the  input  excitation  from  the  base 
and  the  tip  of  the  monopole.  Since  the  pulse  length  is 
much  shorter  than  the  characteristic  time-length  of  the 
monopole,  the  sources  of  radiation  are  separated  In  time 
by  h/c  and  easily  identifiable.  The  first,  third,  fifth, 
etc.,  pulses  radiate  from  the  base,  while  even-numbered 
pulses  are  due  to  end  point  radiation.  A scale  normalized 
to  the  amplitude  of  the  first  pulse  is  included  on  the 
plot  so  that  relative  amplitudes  of  the  radiated  eulses 
can  be  compared.  The  o^p-erinental  results"3  for  this  case 
have  been  reproduced  apd  replotted  in  h/c  time  scale  in 
Figure  3*^* 
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Since  both  computed  and  measured  results  have  been  plotted 
in  h/c , direct  comparison  is  possible.  The  two  results 
have  been  reproduced  in  overlay  format  and  are  shown  in 
figure  3-5.  Direct  amplitude  comparison  between  the  two 
traces  is  not  possible;  however,  amplitude  ratios  between 
pulses  in  each  trace  will  yield  relative  results.  Tn 
general,  the  relative  amplitudes  compare  very  well  with 
experimental  results.  As  pointed  out  by  Schmitt,  et  al., 
Pape  122,  corrections  are  required  in  the  experimental 
results  to  account  for  the  1/R  variation  due  to  base  and 
end  point  displacement.  These  required  corrections  will 
bring  the  amplitudes  into  closer  agreement.  Note  the 
zero  crossovers  and  details  of  the  pulses  even  after  the 
third  pulse.  The  agreement  is  excellent  and  provides  a 
verification  for  the  computed  results.  The  disparity  at 
time  t=3  is  due  to  the  physical  path  difference  between 
the  base  and  tip  to  the  field  point.  In  the  computation, 
this  difference  is  neglected  since  R-*  “ for  far  field 
calculations . 
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COMPUTED  VS  MEASURED  RADIATION  FIELD 
Fi  gure  3 5 


The  computed  results  for  the  rectangular  pulse 
excitation  have  been  plotted  and  are  shown  in  Figure  3.^. 
Note  the  input  pulse  which  is  now  twice  as  long.  as  the 
monopole  time-length.  While  the  effects  of  isolated 
points  of  radiation  are  not  as  obvious  in  this  case  as 
the  previous  one,  the  radiation  field  can  be  relatively 
simply  described.  At  time  t = 0,  the  leading  edge  of  the 
voltage  pulse  emerges  at  the  base  of  the  monopole  and  a 
step  radiation  occurs.  At  time  t = h/c , the  leading  edge 
of  the  pulse  on  the  antenna  is  reflected  from  the  tip  and 
a step  radiation  of  opposite  polarity  and  twice  the  ampli- 
tude occurs.  Consequently,  v;e  note  negative  polarity 
radiation  from  h/c  to  2 h/c.  At  this  time  the  trailing 
edge  of  the  voltage  pulse  is  reflected  from  the  tip  while 
the  leading  edge  is  being  reflected  fro'  the  base  with 
smaller  amplitude.  These  two  sources  of  radiation  which 
occur  at  the  same  time  result  in  a negative  radiation  of 
approximately  twice  the  amplitude  of  the  initial  radiation 
and  last  anproxinately  n/c  seconds  the  time  for  the 
trailing  ed  -e  of  the  excitation  pulse  and  reflection  to 
arrive  back  at  the  base  of  the  nonopole.  Continuing  in 
t.hl3  fashion  and  noting  the  locations  and  reflections 
of  the  leading  and  trailing  edges  of  t.ie  voltage  excitation 
predicts  the  r Hated  field  from  3 h/c  to  5 h/c. 
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The  measured  results  from  Reference  2 have  been  rescaled 
In  h/c  and  plotted  In  Figure  3-7.  Note  the  zero  crossover 
and  similarity  of  the  radiated  field  compared  to  the  com- 
puted results  previously  discussed.  A composite  plot  is 
shown  in  Figure  3.8  which  overlays  the  computed  and 
measured  results.  The  agreement  between  the  plots  is 
excellent.  Correcting  the  tip  radiation  for  the  1/R 
variation  as  discussed  in  Reference  3 would  bring  the 
computed  results  into  even  closer  agreement  in  amplitude 
for  times  greater  than  3 h/c.  Note  again  that  absolute 
amplitude  comparison  between  measured  and  computed  results 
is  not  possible;  however,  normalized  amplitude  comparisons 
are  valid. 

3.2  Receive  Response. 

The  preceding  paragraph  has  presented  the  time 
history  of  the  field  radiated  from  a monopole  antenna 
when  excited  by  a baseband  voltage  pulse.  This  section 
of  the  report  will  present  the  time  domain  response  of 
a monopole  antenna  when  viewed  as  a receiving  antenna. 

,phe  antenna  geometry  for  receiving  antenna  is  shown  in 
figure  3-9-  The  physical  characteristics  of  the  antenna 
are  the  same  as  the  transmitting  case  except  a 50  ohm 
load  is  assumed  connected  at  the  ba ;e  of  the  monopole. 


INCIDENT  VOLTAGE  PULSE 
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RECEIVING  ANTENNA  GEOMETRY 


The  excitation  for  this  case  is  assumed  to  be  a plane 
incident  wave  with  time  shape  identical  to  the  previous 
transmitting  case  (see  Figure  3*2).  The  direction  of 
Incidence  is  broadside  to  the  monopole  as  shown  in 
Figure  3*9* 

The  computed  results  were  obtained  using  the  TIMDOM 
program  in  the  receive  mode.  The  raw  output  data  which 
was  not  included  in  the  report  for  brevity  represents  the 
current  induced  on  the  antenna  by  the  field  excitation. 

If  one  examines  the  current  as  a function  of  time  at  the 
base  segment  of  the  monopole,  he  has  the  time  domain 
receiving  response  of  the  antenna.  Figure  3-10  shows 
the  plot  of  160  time  steps  of  the  base  current  for  the 
narrow  Gaussian  field  excitation  whose  time  dependence 
is  shown  in  Figure  3*2.  The  raw  output  data  was  again 
plotted  in  terms  of  h/c.  Where  h is  the  length  of  the 
monopole,  this  plot  represents  the  current  appearing 
across  the  load  as  a function  of  time.  The  integrating 
characteristics  of  the  antenna  as  a receiver  are  not  too 

i 

obvious  for  the  extremely  short  excitation  field,  but 
will  be  more  obvious  for  the  rectangular  excitation  case 
to  be  discussed  later. 
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However,  the  current  that  appears  at  the  base  of  the  antenna 
is  seen  to  depend  on  the  entire  antenna  length  rather  than 
isolated  points.  The  alternating  polarity  of  the  current 
is  due  to  the  reflection  of  current  from  the  base  and  tip 
and  is  characterized  by  h/c  or  2 h/c  depending  on  the 
time  duration  of  the  excitation  field.  This  will  become 
more  obvious  in  the  next  case.  The  experimental  results 
were  reproduced  and  plotted  in  Figure  3- 11.  Again  the 
comparison  for  time  characteristics  of  the  computed 
results  is  excellent.  The  computed  and  measured  results 
have  been  overlayed  in  Figure  3*12.  A direct  comparison 
shows  the  excellent  agreement  in  relative  amplitude  and 
time  characteristics. 

The  receive  case  was  rerun  for  the  rectangular  pulse 
field  excitation.  As  before,  the  currents  across  the 
load  were  plotted  as  a function  of  h/c  time  and  are  shown 
in  Figure  3-13-  The  integrating  function  of  the  antenna 
as  a receiver  is  now  obvious  where  one  sees  the  triangular 
shape  of  the  response.  The  time  duration  of  the  waveform 
reveals  the  strong  dependence  on  2 h/c.  This  can  be 
explained  as  follows.  The  instant  the  leading  edge  of 
the  field  pulse  irradiates  the  monopole,  currents  are  in- 
duced over  the  full  length  of  the  wire  and  start  to  travel 
in  both  directions  along  the  wire  axis. 
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MEASURED  RECEIVE  WAVEFORM 


INCIDENT 

FIELD 


COMPUTED  VS  MEASURED  RECEIVE  WAVEFORM 


Those  currents  which  travel  toward  the  monopole  tip  are 
reflected  at  the  tip  and  travel  back  toward  the  base. 

Thus,  the  integrating  characteristics  of  the  antenna  last 
for  2 h/c  seconds.  At  that  time  the  trailing  edge  of  the 
field  reaches  the  monopole  and  the  current  polarity  changes 
resulting  in  the  negative  portion  of  the  curve.  Base  and 
tip  reflection  coefficients  cause  an  oscillating  current 
plot,  the  beginning  of  which  is  seen  at  time  slightly 
greater  than  4 h/c. 

The  experimental  results  were  plotted  in  Figure  3-14 
and  overlayed  with  the  computed  results  in  Figure  3-15- 
The  comparison  is  again  outstanding  and  worthy  of  mention. 
The  relative  amplitudes  are  well  predicted  and  the  time 
traces  are  exceptional. 
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SECTIOM  IV  Conclusions  and  Pecommendations 

A computer  routine  based  upon  direct  time  domain 
calculations  has  been  presented  and  discussed.  A walk- 
through description  of  lnputtinr  the  data  and  internretinp 
the  output  has  been  provided  through  an  example  case  to 
familiarize  the  reader  with  the  ease  of  user  Interface. 

The  computer  results  have  been  discussed  and  critically 
compared  with  experimentation.  The  comparison  is  excep- 
tional and  establishes  the  validity  of  the  program.  The 
progran  can  be  used  to  predict  the  transmitting  and 
receiving  properties  of  the  class  of  thin  wire  antennas 
directly  in  the  time  domain.  In  addition,  although  not 
done  in  this  report,  the  scattered  field  can  also  be 
obtained  through  the  program.  The  scattered  field  is  the 
field  reradiated  from  the  antenna  v;hen  irradiated  by  a 
field.  It  results  from  the  radiation  caused  by  the  currents 
induced  on  the  structure  due  to  the  field  irradiation. 

The  utility  of  direct  time  domain  solutions  for 
pulse  excitations  is  receiving  more  and  more  consideration 
in  antenna  and  scatterer  analysis.  Although  this  report 
has  not  attempted  to  establish  the  utility  of  direct  time 
domain  solutions,  its  imnortance  as  an  analytical  tool 
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should  not  be  overlooked.  It  offers  considerably  more 
information  concerning  the  broadband  frequency  performance 
of  an  antenna  than  does  c.w.  analysis.  In  addition,  the 
Fourier  transform  of  the  time  response  (an  option  in  the 
TIMDOM  program)  can  be  used  to  obtain  the  frequency 
response  of  the  antennas  modeled.  If  one  is  concerned  with 
short  pulse  excitations  and  time  waveform  distortions  the 
time  domain  offers  a direct  evalutlon  of  antenna  response, 
however,  if  characteristics  such  as  beam  width,  impedance, 
admittance,  etc.,  are  important,  then  the  frequency  domain 
offers  analytical  advantages.  Nevertheless,  it  is  apparent 
that  direct  time  domain  solutions  offer  a new  look  at 
antenna  performance  and  should  be  useful  to  the  antenna 
analyst  and  designer.  A Fortran  listing  of  the  TIMDOM 
program  is  included  in  this  report  as  Appendix  A.  It 
must  be  reemphasized  that  the  computer  routine  was  pre- 
pared for  user-oriented  ease  of  operation  and  can  be  used 
with  only  minimal  theoretical  or  numerical  techniciues 
understanding. 


58 


dcct  /•  //  / ■'  3 r rnnv 

ULj  1 / it  U..  . . . L./f  / 


SECTION  V References 

(1)  "Time  Domain  Antenna  Study1’  (Final  Report), 
NRL  Contract  J00014-73- C-00?9 

(2)  Harrington,  R.F.,  "Field  Ccmputations  By 
Method  of  Moments",  MacMillan,  JY , ly68 

(3)  Schmitt,  I.J.,  Harrison,  C.VJ.,  and 
.illiams,  Jr.  C.S.,  "Calculated  and 

Experimental  Response  of  Thin  Cylindrical 
Antennas  to  Pulse  Excitation",  IEEE  Trans, 
on  Antennas  & Propagation,  Vol.  AP-ltl, 

March  1066 


APPENDICES 

The  following  appendices  are  included  for  those 
interested  in  implementation  of  the  program  on  specific 
machines . 

Appendix  A contains  a Fortran  listing  of  the  TIMDOM 
program  which  was  written  by  MB  Associates  under  joint 
sponsorship  by  the  Navy  and  Air  Force  under  Contracts 
.J00014-73-C-0099  and  F30602-72X0008 . The  listing 
contains  all  subroutines  in  alphabetical  order  and 
all  are  written  in  Fortran. 

Appendix  B contains  a listing  of  all  machine 
independent  diagnostics  which  appear  in  the  program 
to  aid  the  user. 

Machine  dependent  subroutines  employed  in  the 
program  which  are  specific  to  the  Honeywell/GE635 
computer,  are  identified  and  listed  in  Appendix  C. 
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APPENDIX  A 


C 


1 


Tint  domain  antenna  MODELING  PROGRAM 

COMMON  /DATA/  N,NP,X(200),Y(200),Z(200),SIC200),PIC200),ALP(200), 
I BET  (ZOO)  . I CONI  (ZOO)  , I C 0N2  ( 2 0 0 ) , ITAG(200),IPX,IPY,IRZ 
COMMON  /SCOMP/  SX  (200  ) • SY  (200)  , SZI200) 

COmmun  /CONST/  CDT  # VEL « DT  # NTSTEP 

COMMON  /LOAD/  RES(200),ELD(200),CAP(200),ILOAD 

COMMON  / I OE  LG / NCP mx , ncUMX , IOC » IOCQ» ntma x , NBOUT , JP l , JP2 

COMMON  /maTPaR/  NDLO^S. NPbLK , NLAST , INX 

COMMON  /ARRAY/  CG(bROO) 

COMMON  /ESORC/  ESORC (102R),IFST,NTWaN,mtRan,NTNX,DFRO,ENIn,ENKD,I 
1 EE  N 

COMMON  /ASORC/  EMAG ( I 0 ) , ISRC ( 1 0 ) ,NSRC 
common  /cmat/  cure ( i boo ) . ipizoo j , ix (2oo) 

DIMENSION  CUM(1J,S),  ITST(B) 

INTEGER  HX,HY, HZ 

DATA  MX , HY , MZ/ l MX , 1MY. 1HZ/ 

DATA  !TST/2HCM,2HCt,2MAT,2MRF,2HPC,2MNX,2MfcN,2MEB/ 

R T M 1 s 0 , 

RTM2sO  , 

RTMisO  . 

CALL  SECOND  ( R TM J ) 

NCR  mx  = 1 7600 

NCOMxsbNOO 

JRESRVa I 600 

VEL=2,99HE*8 

PI=3,lul5926 

T‘=PI/I80, 

TD-160./PI 
NTNX=102R 
PRINT  57 

C*LL  SECOND  («TM2) 

IFENaO 

IFSTaO 

JCO-sO 

READ  AND  PRINT  COMMENT  CARDS 

JCOMsJCOM*! 

IF  (JCOM.GT.S)  JCOM*5 

READ  58,  INA, (COM( I , JCOM) , Isl , 1 J) 

PRINT  59,  (COM(l,JCOM),Ial,lJ) 

IF  (INA.EO.ITST(D)  GO  TO  2 
IF  (1NA,E3,ITST(2)>  GO  TO  I 
PRINT  77 
STOP 

HEAD  AND  PRINT  FIRST  OaTa  CARD 

READ  6 J * DT, Tmax.nTSTER, IPX, IPY, IPZ, ILOAD, JPl , JP2 
PRINT  6 m , DT , TmaX.NTSTEP, IPx« IPY, IPZ . ILOAD, JP1 , JP2 
C*LL  IFEM  (IPX) 

C»LL  IFEm  { I P Y ) 

CALL  IFEm  (IPZ) 

CALL  DaTagn  TO  SET  UP  STRUCTU»E  GEOMETRY 


Tr 

1 

T O 

2 

Tr 

3 

T » 

R 

T R 

5 

Jo 

6 

T R 

7 

Jo 

S 

Jo 

9 

Tx 

10 

T o 

1 1 

Tw 

12 

Tr 

13 

T- 

1R 

Tr 

15 

T « 

16 

Tr 

1 7 

TR 

18 

T o 

19 

Jo 

20 

Tr 

21 

Tr 

22 

Tr 

23 

Tr 

2R 

Tw 

25 

TR 

26 

Jo 

27 

Tr 

28 

Tr 

29 

TR 

30 

TR 

31 

Tr 

32 

Tr 

33 

Tr 

3R 

Jo 

35 

TR 

36 

Tr 

37 

Tr 

38 

Tr 

39 

Tr 

<10 

Tr 

Rl 

Tr 

u 2 

Tr 

R 3 

Tr 

RR 

Tr 

RS 

Tr 

Rb 

Tr 

u 7 

Tr 

ue 

T« 

U9 

Tr 

50 

Tr 

51 

Tr 

52 

Tr 

53 

Tr 

5R 

Tr 

55 

Tr 

56 

A- 1 


C *IL  DaTAGN 

Tw 

57 

NP  : N 

Tw 

S8 

PRINT  faO 

Tw 

59 

PRINT  61 

Tw 

60 

SSUMsO, 

Tr 

61 

SMUaO, 

TR 

62 

00  9 Isl  ,N 

Tw 

63 

A0=4(.P(  I ) « TO 

7 * 

69 

B0s8t  T ( I ) • TO 

TR 

65 

SI 1 = SI  ( I ) 

T W 

66 

PRINT  62.  I,XtI),Y(I),2(I),SII,AD,BD.BI(I),ICONim,I,ICON2Cn# 

Tr 

67 

1 ITAG(I) 

Tw 

68 

SSUH  = SSUM*SI  I 

T w 

69 

IP  (SII.GT.SHAX)  SMAXsSII 

Tr 

70 

u 

CONTI NUfc. 

T i* 

71 

PRINT  b5,  N.SSIJH.SMiX 

T w 

72 

c 

T- 

73 

c 

DETERMINE  TIn£  STEPPING  INTERVAL,  OT 

Tw 

79 

c 

TR 

75  : 

IF  (OT.GT.l.E-25)  GO  TO  5 

Th 

76 

OUSRM/VEL 

TW 

77 

SII=10.*»fINT(ALOG10(0T))-3) 

Tw 

78 

OT  = AINT(PT/SI I ♦ , R99 ) • S I I 

Tw 

79 

5 

NTMA  Xi I NT ( T«A  X/DT ) ♦ 1 

TW 

80 

IF  (NTMAX.GT, 1 ) GO  TO  6 

TW 

81 

IF  CNTSTEP.EU.O)  GO  TO  39 

Tw 

82 

GO  TO  a 

Tw 

83 

fa 

IF  (NTSTEP.GT.O)  GO  TO  7 

TW 

89 

NTSTEPsNTRAX 

Tw 

85 

a 

o 

—4 

c 

3> 

Tw 

86 

7 

IF  (NTMAX.LT.NTSTEP)  NTSTEPs.NTRAX 

Tw 

67 

6 

T«aasOT*FlOAJ(nTSTEP-1) 

TW 

88 

NTHANSNTSTEP 

Tw 

69 

if  (nTRAn.GT.NTNX)  NTRaNiNTNX 

Tw 

90 

COT=OT»VEL 

Tw 

91 

PRINT  73,  OT , THAX,NTSTEP 

Tw 

92 

IF  (JPl ,IT,1 ) GO  TO  9 

Tw 

93 

IF  (JP2.LT,  JPl)  JP2SJP1 

Tw 

99 

GO  TO  10 

Tw 

95 

Q 

JP1*1 

TW 

96 

JP2  = N 

Tw 

97 

10 

PRINT  66 

Tw 

98 

IF  (IPX)  11,13,12 

Tw 

99 

i l 

PRINT  67,  HX 

Tw 

100 

GO  TO  13 

TW 

101 

12 

PRINT  68,  HX 

TW 

102 

13 

IF  ( I P V ) 19,16,18 

Tw 

103 

10 

PRINT  67,  HY 

TW 

109 

GO  TO  16 

Tw 

105 

15 

PRINT  68,  HY 

Tw 

1 06 

lb 

IP  l I PZ  ) 17,19,16 

Tw 

107 

IT 

PRINT  67,  H l 

Tw 

108 

GO  TO  1 9 

Tw 

109 

18 

PRINT  68,  HZ 

Tw 

1 10 

10 

IF  (ILO*D.NE.O)  PRINT  69 

Tw 

1 1 1 

IF  (ILOAD.tO.O)  PRINT  70 

Tw 

112 

PRINT  71,  JP 1 » JP2 

Tw 

113 

IP  (ILO»n,EJ.O)  GO  TO  2« 

Tw 

11« 

C 

Tw 

115 

C 

RE  a o InPfcOANCE  LOADING  CARDS 

A- 2 

Tw 

\\b 

w 

c 


PRINT  7 a 
00  20  I ® 1 « N 
HESdJsO, 

ELOdJsO, 

20  C*P(I)=0, 

21  HE *0  75.  J,ITG,II,PE31,ELDI,CaPI 

Is! SEGNO (I  TG » II) 

IF  (*8S(C*PI),LT.l,t-20)  CO  TO  22 
PHINT  76,  ITG,II,I,RESI,ELDI,CAPI 
GO  TO  2 1 

22  PRINT  76,  ITG,II,I,RE3I,EL0I 

21  RES ( I ) sRtSI 

ELDd  >=ELDI 
C*Pd)=C*PI 
IF  (J.Nt.O)  GO  TO  21 

24  C *LL  CONSET 
C 

c caul  esct  (entry  point  of  eincf)  to  read  excitation  data  cards 
c 

CALL  EPtT  ( RE S I ) 

CALL  FfiLOCK  (NHLOKS.NPBLK.NLAST,  IRESRV.N,  INX) 

CALL  COFS  (CURA , CURE, N) 

CALL  FACIO  (CUWF.CufiF.N, IX, IP) 

CALL  SECOND  CRTMI) 

RTmsRTh1-RTM2 
PRIM  46,  RTM 

c 

C SOLVE  FOW  STRUCTURE  CURRENTS 

C 

call  rsoL  (cu*hix,  ip.n) 
call  slcund  c r t m j 

HTMsRTm-RTHI 

PRINT  4 7 , RTM 

25  REaO  «0,  INA,ITMl,nN2,lTMJ,ITMu,TMI»TM2#TMJ#TH4 
C*LL  SECOND  (RTm) 

RTMsRTm-HTHI 

PRIM  48,  RTR 

PRINT  Ul,  I NA,  ITM1  , ITE'2,  I TMi,  1 TRR,  T«1  , T*2,  TMS,  TM4 
IF  dNA.tO.ITST(l))  CO  TO  26 
IF  dr,A,EQ.lTST(U))  GO  TO  27 
IF  (INa ,t  Q. I TST (5) ) GO  TU  28 
IF  (!NA,t!3,ITST(6))  GO  TO  1 
IF  (1NA.EQ.ITST (7))  STOP 
IF  (INa.EQ, ITST(b))  GO  TO  29 
PRINT  42 
STOP 
C 

c COMPUTE  ANTENNA  INPUT  ADMITTANCE 

c * 

26  lalSECNOI 1TM1 . IT«2) 

CALL  AnTRAN  (I,T“|,ITM1) 

GO  TO  26 

C 

C COMPUTE  RADIATED  FIELDS 

C 

27  CALL  REPAT  ( I TMJ , p M2, 1 ThI,  I TMu , TMl  , TM2 , TR J, T««) 

GO  TU  25 

C 

c PUNCH  CURRENTS 


TW 

117 

TH 

118 

Tw 

119 

T« 

120 

Tw 

121 

TH 

122 

T *i 

121 

Th 

124 

Th 

125 

Th 

126 

Tw 

127 

Th 

128 

Tw 

129 

Tw 

1 SO 

Tw 

111 

Th 

112 

Th 

111 

Th 

114 

Tw 

115 

Tw 

116 

Th 

117 

Th 

138 

Th 

119 

Th 

1U0 

Th 

1 A»  1 

Th 

142 

Th 

1 R 3 

Th 

144 

Tw 

145 

Th 

146 

Tw 

147 

Tw 

148 

Tw 

149 

T w 

ISO 

Tw 

151 

Tw 

152 

TH 

151 

Tw 

154 

Tw 

155 

Tw 

156 

TH 

157 

Th 

158 

Tw 

159 

Tw 

160 

T w 

161 

Th 

162 

Tw 

161 

th 

164 

Th 

165 

Tw 

166 

Tw 

167 

Tw 

168 

T w 

169 

Th 

170 

Th 

171 

Th 

172 

Tw 

1 71 

Th 

174 

Tw 

175 

Tw 

176 
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c 

TM 

177 

26 

1=ISEGnO(iT«1,ITM2) 

Tn 

170 

C*LL  CGET  ( l ,CURE ,nT»an) 

Tm 

179 

PUNCH  43,  (COR( J, 1 ), Jsl,  13) 

Tw 

1«0 

PUNCH  44,  NTSTEP, I , I TM1 , ITM2 

T » 

101 

PUNCH  45 , (CU«E(J),Jcl#NT3T£P) 

Tw 

102 

GO  TO  25 

Tm 

103 

C 

Tm 

1 04 

c 

COhPuTE  ENtRGT  0UOGET 

Tm 

105 

c 

Tm 

100 

29 

PRINT  49 

Tm 

167 

ITM2«0 

TM 

109 

ENINsO  . 

Tm 

1 SR 

ENLS=0, 

Tm 

IRQ 

IE  (NSWC.LT.  1 ) GO  TO  33 

Tm 

1 9 1 

IE  (IEST.NE.01  GO  TO  32 

Tm 

192 

I EE  Ns  1 

Tm 

1 9 3 

P W I N T 50 

TM 

194 

00  31  1=1 ,NSRC 

Tm 

195 

issisRcm 

Tm 

l R6 

vlt=emag(  n*si  ( is) 

Tw 

197 

CALL  CGET  ( IS/CURE.NTHAN) 

Tw 

190 

sum=o, 

TM 

199 

DO  30  J = 1 , nTWAN 

T « 

200 

PESI=ESOWC(J) 

Tm 

201 

IE  (PESI  ,LT  ,1  , E * 20  ) GO  TO  30 

Tm 

202 

SUm*SUH»we3I»CU«E(J) 

Tm 

203 

30 

CONTINUE 

TM 

204 

Suh=SUh*OT*VLT 

Tm 

205 

PRINT  SI,  I , IS, SUM 

Tm 

206 

31 

EN I Ni E N I N »SUM 

Tm 

207 

GO  TO  33 

Tm 

200 

32 

PRINT  56 

TW 

209 

33 

IE  (ILO*0PEQ,OJ  GO  TO  37 

TM 

210 

00  30  1=1, N 

Tm 

211 

IE  CRESm.LT.l.E-20)  GO  TO  30 

Tm 

212 

IF  CIW.EU.I)  GO  TO  3U 

Tm 

213 

PRINT  52 

Tm 

2 1 4 

I T m2 s 1 

TW 

215 

34 

CALL  CGET  (I,CUWF,NTRAN) 

Tm 

216 

SUM®  0 . 

Tm 

217 

00  35  Jsi.NTWAN 

Tm 

216 

55 

5UM35UMACURE ( J5*CURE (J) 

Tm 

219 

SUMsSUM»DT«WES(  I ) 

Tm 

220 

PRINT  53,  I , R£  £ ( I ) , SUM 

Tm 

221 

ENLS=ENLStSUM 

TM 

222 

36 

CONTINUE 

TM 

223 

37 

enwo=enin-enls 

TM 

224 

PRINT  5«,  EMN,tNLS,ENRO 

Im 

225 

IF  (ENWO.LT , 1 .1-20)  I E E N « 0 

Tm 

226 

IE  (ENIN.GT,  1 ,E-20)  GO  TO  38 

Tm 

227 

IFENsO 

Tm 

226 

GO  TO  25 

Tm 

229 

3« 

T * 1 = 1 0 0 , *E  NR  D/EN  I N 

Tm 

230 

PRINT  55,  Tm l 

Tm 

231 

GO  TO  25 

Tm 

232 

3* 

PRINT  72 

Tm 

233 

STOP 

Tm 

234 

C 

Tw 

235 

C 

Tm 

236 
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UO  FORMAT  <A2,I3, 3IS.6E10.J)  T R 

9 t FORMAT  ( //, IX, 15h*»»  DATA  C A 9D * * , A ? , J 3 , 3 I 5 , 6E I 2 . 5 ) TR 

92  FORMAT  ( //, 1 X , 3»HINV»LI0  Data  CARO  LABEL  AFTER  SOLUTION)  Tr 

9 3 FURMAT  (lJAb)  T R 

49  FORMAT  (15,21 MClIRRENTS  FROM  SEGMENT,  15, 6X,4HTAf.s,  I5,2x,  IOHInCREME  Tr 

INIs, IS)  T R 

«S  FORMAT  l to£ 1 2,5 ) Tr 

U6  FORMAT  ( lX,i7nTIME  FOR  INITIALISATION  OF  CONS T AN T S= , F 9 . J , SH  SEC  Tr 

1 , ) Tr 

9 7 FORMAT  ( I X , 2 6 u T IRE  FOR  CURRENT  SOLUTIONS, F9, 3, 5H  SEC.)  Tr 

FORMAT  ( / , 1 « , 1 JrR'JNNI  NO  TIMEr,,C9,3(«;M  SEC,)  Tw 

09  FORMAT  ( ///,3bX,25H-  - « ENERGY  8U0GET  - - -)  Tr 

50  FORMAT  ( //,33X, RESOURCES  - , // , 35 X , 1 SMSOUHC E SEGMENT, 3X,  12HENER  T r 

) G Y INPUT, /, 37X, 3MNO, , bX , JMNO, , 6X, «M< JOULES) ) Tr 

51  FORMAT  (3SX,  IS.uX, IS.SX.El  1 ,R)  Tw 

52  FORMAT  ( //,33<,7HLOaOS  - , / , iSx  , TNSEGMEN  T , 3 X , I OMRE  5 1 S T A a'CE  , 2X  , H Tr 

JmENERGY  LOSS, /, 3SX, 3hnu, , 7x , 6H (ohmS) ,5x, bM( JOULES) ) Tr 

53  FORMAT  (3bX, IS,2X,2E12,«)  TR 

5a  FORMAT  ( / , 36X , 2GMTGT  AL  ENERGY  INPUT  o,Elt,a,7H  JOULE S ,/, 3b X , 22M  T« 

lENtRGY  LOST  IN  LOADS  =,Eil,4,7H  JOUU  5 , / , 36  X , 23HTOT  AL  ENERGY  RADIA  Tr 
2 T E 0 =,t I 1 , a, 7M  JOULES)  TR 

55  FORMAT  ( 36X, 2RHTIME  DOMAIN  EFFICIENCY  =,F7,2.8H  PE»CtNT)  Tn 

56  FORMAT  ! I , SRM I NPUT  ENERGY  NUT  COMPUTED  SINCE  SOURCE  MAS  SEEN  TR  T* 

IANSFORMEO)  Tr 

57  FO«maT  (1h1,  /,2PX,a2n........ ......................... .........  t» 

I,/  /, 31  X, 36HT  Ime  DOMAIN  ANTENNA  MOUELING  PROGR A m , / / , as x . U?M« * » * » * • TR 
.....A.,....*..*....,....,////)  Tr 


58  FORMAT  ( A 2 , 1 3 A 6 ) TR 

59  FORMAT  (28X,13Ab)  TR 

60  FORMAT  (////  33X.33H-  ---  SEGMENTATION  DATA  - - - »,  //,U0X,  Tr 

1 2 1 hCOORD  I n a TE  S In  mETERS,  //,25X,  Tr 

2 50m I ♦ and  I-  INDICATE  TM£  SEGMENTS  BEFORE  and  AFTER  I,//)  Tr 

61  FORMAT  (?X,aHSEG, , 3X , 26MCOOHD I NA  TE  S OF  SEG.  CENTER , 5X , UHSEG . , Tr 

1 5 x , 1 8 HOR 1 E N T A T 1 ON  ANGLE S , a x , a«R I HE , 5 X , 1 SHCUNNEC T ION  DATA, ax,  Tw 

2 JmTAG*/»2X,Shg0,,Tx»IHX,9X,1hy,9x,1HZ,7A,6HL£NGTh,5X,5rAl2Ma,  Tr 

3 5X, awBETA, 6X,6‘iRADIUS, UX , 2m(., ax, lMt ,5X , 2HI ♦ , ax, 3HNO.  ) Tr 


62 

FORMAT 

(IX 

,15, OF  10. S, IX, 3F 10.5,316,2*, 15) 

TR 

63 

FORMAT 

(2110,3, 15, 2x,SAl,3IS) 

Tr 

64 

FORMAT 

( 

//,lx,2lM*».  FIRST  DATA  CaRD**,2E12.5.!5,2X,3A1,3I5) 

Tr 

65 

format 

( 

//, lx, ) 9HNUM01 R OF  SEGMENTS*, 17,/, IX. 18HT0TAL  Rlwt  LENG 

TR 

l T ms , F 10, 

b , 3 M M, , / , 1 X , 23MMAXIMUM  SEGMENT  LENGTHS, F10, 5, 3M  m,) 

Tr 

66 

FO-maT 

( 

/ / / , l x , 2 1 hurt  ions  selected-  - • ,/) 

Tr 

67 

FORma  T 

( 

5M  THE  , A 1 , 3 7m*0  plane  IS  a MAGNETIC  SYMMETRY  PLANE) 

Tr 

68 

format 

( 

5m  ThE  , a i , 30h=O  PLANE  IS  AN  ELECTRIC  SYMMETRY  PLANE) 

Tr 

69 

format 

( 

1 x , 27HMQDEL  MAS  IMPEDANCE  LOADING) 

T * 

70 

format 

( 

1 » , 1 ohmejoe L IS  NOT  LUAOEO) 

T « 

71 

FORMAT 

( 

1X,37^CUhREnTS  -ILL  0E  PRINTED  FOR  SEGMENTS,  15, 8W  ThROU 

T . 

IGH  ,13) 

T * 

72 

F ORMA  T 

(IX 

, 25wZE  RO  TIME  STEPS  HEOUESTtO) 

Tr 

73 

FORMAT 

f 

T * 

lx,  2 JMT  J ►'£ 

STEPPING  INTERVAL*, fc 12,5,5m  St C . , / , 1 X , 1 3h  . a x l HUM  TImEs, 

T* 

2 El?, 5 

C i_, 
• J • • 

SEC.,/,’.X,2lHNuMBEi5  Or  TIME  GTE  PC* , I a, // , a 1 U ......... 

* 

• n 

3* <««.».. 

* * * 

TR 

74 

FO-MAT 

( 

///,35x,33m«  • • * STPuCTu-E  LOADING  « - - -,//,23«, 

T R 

1 6 0 m T A G 

Increment  SEG.  RESISTANCE  inductance  capacitance 

Tr 

2,/  . 2 3 x , 

JMNOt,lOX,3MNO,,5X,6M(LlMMS),6X,8M{HF.:jRYS),bX,0H(FARAO8)) 

T * 

75 

FORMAT 

(315, 5«,3E 10,3) 

Tr 

76 

f * a T 

(21 

»,I3,3x,I5,2x,I5,3»,fcl0.9.iA,tlO,M,u*,El0.a) 

Tr 

77 

FORMAT 

( 1 x 

, 3«MJNC0RREC T LABEL  FOR  A CUM-ENT  CARO) 

Tr 

F NO 

Tr 

237 

236 

239 

290 

291 

292 

293 

29U 

295 

296 

297 
2«0 
299 

250 

251 

252 

253 
2Su 

255 

256 

257 

258 

259 

260 
261 
262 
263 
269 
265 
2*6 

267 

268 

269 

270 
27) 

272 

273 
27a 

275 

276 

277 

278 

279 

280 
281 
202 
203 
<rea 
205 
2*s 
287 
2G  8 
2«9 

290 

291 

292 

293 
XYU 

295 

206- 
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SUBROUTINE  anTHAn  ( ISEG,ANOH, IPCH) 

AN 

1 

AN 

2 

antran  computes  ano  prints  the  INPUT  ADMITTANCE 

AND  IMPEDANCt  AT 

AN 

3 

segment  ISEG 

AN 

4 

AN 

5 

common  /OaTA/  N, nP, x (200) . Y(200) , Z (200) , SI (200) 

,«!<200)»ALP(200), 

AN 

6 

1 BET  (20  0),  ICON!  (POO),  I 00*2  ( 20  0 ) , I T AG  ( 2 0 0 ) , I P X 

, IPY, IPZ 

AN 

7 

COMMON  /ESORC/  ESOPCl J02U) , IEST,NTHAN,MTRAN,NTNX,DFRO,tNIN,ENR0, I 

AN 

8 

1FEN 

AN 

9 

COMMON  /CHAT/  CURE ( 1600) . IP(200) , IX (200) 

AN 

10 

COMMON  /CONST/  CDT,VEL,DT,NTSTEP 

AN 

1 1 

CUmmun  /ASORC/  EMAG ( I 0 ) , ISRC ( 1 0) , NSRC 

AN 

12 

Common  /SCRATM/  PnW ( 6 1 2 ) 

AN 

13 

COMPLEX  TMIT,ZPED,ZN0R 

AN 

19 

DATA  TSMIN/S,/ 

AN 

IS 

I A (1EST.EU.1)  GO  TO  2 

AN 

1 6 

AN 

17 

FOURIER  TRANSFORM  SOURCE 

AN 

18 

AN 

19 

MTRANsALOG(ELOAT(NTRAN) )/.69il 9718M  .5 

AN 

20 

CAUL  ITOF  (ESU«C,NTRAN, MTRAN, 1 , 1 , NTNX.NTNR) 

AN 

?l 

I F S T = 1 

AN 

22 

DFRQ=l,/(OT*FLUAT(2**MTRAN)) 

AN 

23 

NTNRsNTNR/2 

AN 

29 

DO  1 I = l , NTNR 

AN 

2S 

PWH( I )r0. 

AN 

26 

CALL  CGET  (ISEG.CURE.NTHAN) 

AN 

27 

AN 

28 

FOURIER  TRANSFORM  CURRENT 

AN 

29 

AN 

30 

CALL  ITOF  (CURF.NTRAN.NTHANfO/O.NTNX^TNR) 

AN 

31 

I NOR “ 0 

AN 

32 

IF  (ABS(ANOR) ,GT.  1 ,E-20)  INORsi 

AN 

33 

IF  (NSRC.LT.l)  GO  TO  a 

AN 

39 

DO  3 I s l # NSRC 

AN 

35 

I I 3 I 

AN 

3b 

IF  (IS»C(I) ,EO, ISEG)  GO  TO  5 

AN 

37 

continue 

AN 

3G 

VL  T s 1 , 

AN 

39 

1 SPRsO 

AN 

UO 

GO  TO  6 

AN 

Ml 

VLTsFMAG(H)*SI(ISEG) 

AN 

92 

ISP«=  . 

AN 

93 

PRINT  10,  ISEGrVUT 

AN 

9U 

1*  UNO0. EG,!)  PRINT  1!,  A NOR 

AN 

“5 

PRINT  l 2 

an 

96 

OF  =DF  ROM  ,t-6 

an 

u r 

F = -0F 

AN 

U A 

FMAXs|,E.6/(TS«rN»0T) 

AN 

u9 

NTnR*ntnr/2 

AN 

SO 

if  (ipcm.nf.0)  Punch  15 

AN 

SI 

DO  8 1 al , NTNH 

AN 

52 

F cF *OF 

An 

Si 

IF  (F.GT.fMAx)  GO  TO  9 

AN 

59 

12  = 2*1 

AN 

SS 

II "12-1 

AN 

S 6 

A- 6 
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10 


1 1 
12 


13 

*• 

15 


r«l ToCMPLX (ESORC C 1 1 ) , ESORC (12) ) *VUT 
SPC  T =C  A8S ( vm i T } 

YMI TsCMPLX(Cu«F(I 1 ).CUHF(12) J/YMIT 

2p£0s l , / YM J T 

IF  (ISPH.EQ.l)  P*R(1)sP*H(I)f.5*VLT*vLT*REAL(YMIT) 
IF  (IPCH.Mt.O)  PUNCH  l «*,  I,F,2PE0,YHIT 
IF  (INOW.EO.O)  GO  TO  7 
ZNUR  = Z»>FD/anoR 

P»INT  li,  I,F,ZPEO,YHIT,ZNOR,SPCT 
GO  TO  8 

PRINT  l<t,  I,F,ZPtO.VMIT,3PCT 

CONTINUE 

RETURN 


FORMAT  ( //,2UX,S0H.  - - ANTENNA  INPUT  impedance  and  admittance 

1-  - -»//#11X,23h50URCE  SEGMENT  NUMBER  = , 1 5 , / # 3 1 X , 2 1 MPE AK  SOURCt  VO 
ZLTAGE  = #fcll.«.t»M  VOLTS) 

FORMAT  ( 31 X | 32MIMPEDANCF  NORMALIZATION  FACTOR  s.Fll,«,5H  OHMS) 
FORMAT  ( /,  lX,MHSTEP.«X,0HFRtli0FNCY,6X.  16HIMREDAr.CE  CUrU'o  ) . 1 0 X , 

1 1 7MA0MJ TTANCE  CmmuS),7X,?OMnORMALIZED  I mreDa NCE , 6 x , 6HS0I'RCK / / , 

2 2X, 3MNO, , b*,5M(MRZ) # J(«X, UMHEALf  8X,SMIMAG. * 1 X) , :,8MSPECTRUM) 

FORMAT  (lX,Iu,2tl3.«.E12.R»2(Et«.R»tl2J<n#2X,En.«) 

FORMAT  (lx,Iu#<iElJ.«,tl2.u,Eia.«,Fl2.R,28x,Ell.a) 

FORMAT  (NX. RmFhEUUENCY.RX, RM IMPED A NCE# 1 8X , 1 OH  ADM  I TTANCE  ) 

END 


AN 

57 

AN 

58 

AN 

5R 

AN 

60 

AN 

61 

AN 

62 

AN 

6 3 

AN 

6 0 

AN 

65 

AN 

66 

AN 

67 

AN 

68 

AN 

69 

AN 

70 

AN 

71 

AN 

72 

AN 

73 

AN 

7 « 

A N 

75 

AN 

76 

AN 

77 

AN 

78 

AN 

79 

AN 

80 

AN 

81 

AN 

82- 
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function  atGn?  c * , y j 

*t&n^  15  * rc  t anger  r function  modified  to  return  o,  when  xtyoo 

1*  (*)  3,1,5 

1 IF  (Y)  3,?,3 

2 ‘TGNtJsO, 

RETURN 

3 ATGNc-aUn*  (X,Y) 

RETURN 
ENO 


o o o o 


SUBROUTINE  CGET  ( I SE  G , C SEG , NI.M ) CG  t 

CG  a 

CGET  Fries  4KRAV  CSEG  WITH  THE  CURRENT  ON  SEGMENT  ISEG  FOR  *11  CG  3 

TIME  STER3  CO  a 

CG  5 

common  /data/  N,Np,x(aoo),Y(aoo),zcaoo),stUoo).t!i(2oo),AiP(dOo),  c&  6 

1 8FT(20n),lCON|<20(i),ICON2(200)fIT4G(200),lPx,IPY,lPZ  CG  7 

COMMON  /ARRAY/  CG(oOOO)  CG  8 

COMMON  /IOEIG/  NCEMX, NCGMX, IOC, roCG*NTMAx, NSCJUT, JPI , JP2  CG  9 

DIMENSION  CSEG(102«)  CG  10 

II=N«2  CG  H 

1Bc2»ISEG-1  CG  12 

HIP-II  CG  13 

n0C0e999999  CG  19 

IF  (IOCG.EQ.O)  GO  TO  1 CG  15 

REWIND  11  CG  16 

READ  (11)  NOCQ, CCO(K),k=1,NOCU)  CG  17 

1 DO  3 Jsl.NLH  CG  18 

I a I ♦ 1 1 CG  19 

IF  CI.U. NOCQ)  GO  TO  2 CG  20 

READ  (11)  NOCQ. (COCK) ,K=1 , NOCQ)  CG  21 

I=IB  CG  22 

2 CSEG(J)-CQCI)  CG  23 

3 CONTINUE  CG  29 

IF  (IOCO.NE.O)  REWIND  11  CG  25 

RETURN  CG  26 

END  CG  27- 
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c 

c 

c 


1 

2 

c 

c 

c 
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SUBROUTINE  COES  ( C URF , CURL , HO  I » ) CF 

CF 

COFS  SETS  UP  THE  ARRAYS  OF  INTERACTION  COEFFICIENTS  FOR  BOTH  CF 

PRESENT  TIME  AND  RETARDED  TIME  INTERACTION  CF 

CF 

COMMON  /DATA/  N,NP,x{200),Y(200><Z(200),SI(200),dI(200).AI.P(200),  CF 
1 8ET(200),lCONtC200),ICON2C200),ITAG(200),IPX,IPY,IPZ  CF 

COMMON  /SCOMR/  SX(200),3Y(200),3Z(200)  CF 

COMMON  /CONST/  COT, VEL,DT,NTSTEP  CF 

COMMON  /array/  CJ(bROO)  CF 

COMMON  /AWRAx/  Oll,QI2,QI3,EC(200,5).EO(200,5),IHET(200),HFR(lSu00  CF 

1 ) CF 

COMMON  /EmaTs/  EMXC3, J)#EMy(J, J),EmZC3,1),qmx(j,3),qmy(J,j),(!MZ(s,  cf 
13)  CF 

COMMON  /load/  RESC200  ) .ELDC200) ,CAP(200) , ILOAO  CF 

Common  / 1 nF  L D / nCFmx,ncomx» IuC, IOCO/ nTma  x , nRuUT  » JP 1 » JP2  CF 

COMMON  /maTPaR/  NBlOkS, NPBL* # NlAST , INT  CF 

DIMENSION  CURFINDP’.NDI-),  CURE(I)  CF 

DIMENSION  CUE  C 17600) * ICOF(17hOO>,  SHT(tuOO),  ISRTI6900)  CF 

DIMENSION  ALF  a ( J , } ) , 8E  T a ( 3 » 3 1 CF 

EQUIVALENCE  C COF , i C ) > (IC0F,EC),  (S»T,CQ),  (ISRT.CO)  CF 

LPx=l  CF 

LPY  s l CF 

LPZ*1  CF 

IF  (IPX.NE.O)  LPX=2  CF 

IF  (IPY.NE.O)  LPY=2  CF 

IF  IIPZ.NE.O)  LPZ=2  CF 

ICnTCsO  CF 

I C NT  5 = 0 CF 

I C TUT  a 0 CF 

IOC=0  CF 

1 0 S s 0 CF 

ICOF (NCFMX)=9B9BR89  CF 

ISRT(NCQMX)=9898989  CF 

IST=2»N*LPX*LPT«LPZ  CF 

ICMAX=NCFMX-I8r  CF 

I$mAxbnCOMX»IST  CF 

ISTsl  CF 

N T M A x a o CF 

IF  (INT.FQ.O)  CO  TO  1 CF 

REMIND  13  i CF 

nObl*  SsNBLUKS  CF 

I2aNPBLR»N  CF 

1 T =npbl*  CF 

GO  TO  ? Cf 

NObLKS=l  CF 

I T a N CF 

1 » 0 CF 

CF 

b E G 1 N LOUP  OVER  FILED  EVALUATION  POINTS  CF 

CF 

DO  30  I UHL  * Sal#  NOBL  * S CF 

lF(jNt,tv,0)GOfui  CF 

IF  ( lUbU'xS.EQ.NOBtKS)  IT*nla3T  CF 

DO  2»  19*1  , IT  CF 

I«I*1  CF 


1 

2 
3 
a 

5 

6 

7 

8 
9 

1 0 
1 l 
12 
13 
19 

15 

16 

17 

18 

19 

20 
21 
22 
21 
29 

25 

26 

27 

28 

29 

30 

31 

32 

33 
39 
19A 

35 

36 

37 

38 

39 
UO 

91 

92 

93 
u a 

95 

96 

97 

96 
U 9 

50 

51 

52 

53 
59 
55 
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X I ■ X c I ) 

CF 

56 

V I * V ( I ) 

CF 

57 

ZlaZU) 

CF 

58 

8XI*SX(I) 

CF 

59 

SYIeSY(I) 

CF 

60 

SZXSSZ(I) 

CF 

61 

00  4 Jal  .N 

CF 

62 

CUKF (J, 1 8 ) a 0 , 

CF 

63 

IFL  = 0 

CF 

64 

CF 

65 

begin  loop  OVER  SOURCE  segments  including  IMAGES 

CF 

66 

C f 

67 

00  29  J X a j , LPX 

CF 

68 

RFx=FLO*T(3-JX»2) 

CF 

69 

SFXeRFx 

CF 

70 

IF  (IPX.LT.O)  SF  X a 1 t 

CF 

71 

00  29  JY=t,LPY 

CF 

72 

RFyrFLOXT ( 3-JY*2) 

CF 

73 

SFYsRFY 

CF 

74 

IF  (IPY.LT. 0)  SF  Ye  J , 

CF 

75 

00  29  JZ=1,LPZ 

CF 

76 

RFZ=FLO*T(3-JZ*2) 

CF 

77 

SFZsRF l 

CF 

78 

IF  (IPZ.LT.O)  SFZsl, 

CF 

79 

RFL=SFX*SFY*SF z 

CF 

80 

NTL  = 0 

CF 

81 

00  5 J=1 ,N 

CF 

82 

HXaxj-X(J)*RFX 

CF 

83 

RY=YI-Y(J)«MFY 

CF 

84 

RZ=Zl*ZC J)»RFZ 

CF 

85 

resort  chx*mx»hy»«y*rz*rzj 

CF 

86 

KsH/COT*.S 

CF 

87 

IF  (K.EQ.O)  Kal 

CF 

88 

IRET(J)»K 

CF 

89 

IF  (K.GT.NTU  NTLsK 

CF 

90 

00  5 Msi,5 

CF 

91 

EC(J,M)s0, 

CF 

92 

EQ(J,M)so. 

CF 

93 

00  21  J 3 1 * N 

CF 

94 

HX s X I “X ( J ) *RF  X 

CF 

95 

RYsYI»Y(J)*RFY 

CF 

96 

HZ=ZI-Z( J)*RFZ 

CF 

97 

r2=«x»rx*ry*RYtRZ*RZ 

CF 

98 

ReSQRT (R2) 

CF 

99 

KsIRET(J) 

CF 

100 

TAUeOT* FLOAT ( K ) 

CF 

101 

SX.)  = SX(J)*RFX 

CF 

102 

SYJ*SY(J)*RFY 

CF 

103 

SZJ3SZ(J)*RFZ 

CF 

109 

CALL  EmAT  (RX,RY,RZ#H2,TAU,SxJ,SYJ.SZJ»J) 

CF 

1 05 

00  6 L3 I i 3 

CF 

106 

DO  <>  M=i  , j 

CF 

107 

ALFA£L.m)=(EmX(l.M)*SXItEmY(L(u)*SY!tEf:U.u5‘SZI5*RFL 

CF 

103 

BfcTA(L,M)=(OMX(L,M)«SXI*OMV(L,M)*SYI»OMZ(L,M).sZI)*RH 

CF 

109 

JCt=IC0Nl ( J) 

CF 

1 1 0 

IF  CJCI.NE.O)  GO  TO  7 

CF 

1 1 1 

K KMs  0 

CF 

112 

GO  TO  11 

CF 

113 

IF  (JCl.LT. 19000)  GO  TO  6 

CF 

1 lM 

SIG1eFlOat(JCI-20000) 

CF 

115 

A-ll 


a 


JC  t » J CF 

CO  TO  10  CF 

IF  (IC0M2(JC1) ,NE , J)  CO  TO  9 CF 

S I C 1 3 1 , CF 

GO  TO  10  CF 

9 IF  (ICONt(JCt).Nt.J)  CO  TO  RR  CF 

SIC13-1.  CF 

10  K*=IWET(JC1)  CF 

KKM=KM-K  CF 

IF  (IAB3(KK«)  ,LE.  1)  GO  TO  It  CF 

PRINT  05,  I.J.JCl.KKH  CF 

STOP  CF 

11  JC2=ICON?(J)  CF 

IF  (JC2.NE.0)  CO  TO  12  CF 

KKPsO  CF 

GO  TO  lb  CF 

12  IF  CJC2.LT,  190001  GO  TO  15  CF 

SIG2sFLOAT( JC2-20000)  CF 

JC2sJ  CF 

GO  TO  15  CF 

13  IF  ( I CON 1 { JC2 1 , NE , J 1 GO  TO  IR  CF 

SIG2«1,  CF 

GO  TO  15  CF 

1 R IF  ( ICON2( JC2) ,N£ . J)  GO  TO  RR  CF 

SIG2=-1.  'n  CF 

15  KP=IR£T(JC2)  CF 

KKP=KP»K  CF 

IF  (IA8S(KKP),LE.l)  GO  TO  lb  CF 

PRINT  R5,  I,J,JC2,*KP  CF 

STOP  CF 

lb  00  20  N=l  , 3 CF 

K X X 3 M ♦ 1 CF 

KMX=KKM+KXX  CF 

kPX3*KP*KXX  CF 

KKs*-H*2  CF 

IF  (KK.EO.O)  GO  TO  18  CF 

IF  ( JC  1 ,t<3. 0)  GO  TO  1 7 CF 

EC(JCl,,<*x)sEC(JCl,KMX)»At.FA(l,M)*srGl  CF 

E0(JCl,KPx)=EQ(JCl,KMX)+bETA(l,F')*3IGl  CF 

17  EC(J,KXX)sEC(J,KXX)*ALF A(2,H)  CF 

EQ(J,«X* )=tU( J,«XX3*BETA(2,M)  CF 

IF  (JC2.E9.0J  GO  TO  20  CF 

tec JC2, KPX)=EC t JC2»KPX)*ALFA(3,H)»SIG2  ' CF 

tQ(JC2.*Px)itO(JC2»*PX)*bETA(i,rt)*$lG2  CF 

GO  TO  20  Cr 

18  IF  (JC1.EO.O)  GO  TU  1 9 CF 

CURF  (JC  1 . Its)  sCUHf  (JC1  , IB)-(ALF  A ( U J)  ♦ BETA  ( 1 , 31*0133  *SIG1  CF 

kOCJCl  ,K"Xl=tU(JCl  ,K*iX)*BETA(  1 , MJ.SIG1  CF 

19  CUHF  {j,  IHISCURF  (J,  JP.J-ALFA  ( 2 , 3 ) -BE  T a ( 2 , 3)*QI  J CF 

EO<J.FXX)rta(J.Kxx)*RflA(a,M)  CF 

IF  (JC2.EO.O)  GO  TO  20  CF 

CURF  CJC2. Ib) sCURF (JC2# IB)- ( ALF  a ( J, J) *BtTA{ 3, 31*0131*3102  CF 

EO(JC2,KPX)stO{JC2>KPX)*(lETA(3,M)*slG2  CF 

20  CONTINUE  CF 

21  CONTINUE  CF 

C CF 

C ADO  TERNS  HUE  TO  LOADING  CF 

C CF 

IF  (ILOAD.tU.Ol  GO  TO  22  CF 

IF  ( IF L .EO.  1 1 GO  TO  22  CF 


1 1 b 
1 1 7 

i ie 

119 

120 
121 
122 
1?) 
1 2R 
125 
12b 

127 

128 
129 
1 30 

131 

132 

133 
1 3R 
1 35 
1 3b 
137 
1 38 
139 
1 R 0 
1 R l 
1 R 2 
1 R 3 
1 RR 
1 R5 
1 R b 
1 R 7 
1 R 8 
1 U9 

150 

151 

152 

153 
15R 
155 
15b 

157 

158 

159 
l b 0 
lbl 
162 
1 h 3 
1 b u 
1 65 
lbb 
167 
lb« 
Ib9 

170 

171 
1 72 
173 
1 7a 
1 75 
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EC(I,2)*tCCI,2)'ELD(I)«.5/(DT*Sl(I)> 

EC(I,  J)«EC(I#3)»ELP(I)«2./(0T*SItn) 

CURF ( I , ImJeCURF  (I , IR)»(RESCI }*ELD(1  )*1  ,5/DT)/SI  (I) 
IF  t C A P ( I ) ,IT,  1 ,E-20)  GO  TO  22 
CURFU,  IB)sCUHF(I,  IB)»OU/(CAP(I)*SI(I)) 
EQ(I,«)=t0(I,<O-l,/(CAP(I)«SI(I)) 

22  NTLbNTL+2 

IF  (NTL.GT.NTMAX)  NTMAXsNTL 

C 

C SOP T COEFFICIENTS  FOR  FIELD  EVALUATION  SEGMENT  I 

c 

DO  28  LL= l , NTL 

l=ntl-ll 

l HDG  = 0 
JaO 

23  J = JM 

IF  (J.GT.N)  GO  TO  28 

K=IRET( J) 

K0I; =K-L 

IF  ( iABS(KDIF)  ,GT.2)  GO  TO  23 
IF  (IHOG.EQ.l)  GO  TO  2a 
I HDG=  1 

I5«T(IST)s-(L»10000M> 

IST=IST+1 
2a  KKsJ 

KKPa 1ST 
ISTrlST*! 

25  SRT(lST)=ECCJ,KDIF+3) 

S»  T ( I S T ♦ 1 )=EO(J»kOIFp3) 

ISTbIST+2 

J*J*  1 

IF  (J.GT.N)  GO  TO  26 
KslMET(J) 

KDIF=K-L 

IF  (IAHS(koIF)  .LT.3)  GO  TO  25“ 

26  IS»T(KKP)=20000*(J-KK)+2*KK*t 
IF  (1ST, LT.  ISMAX)  GO  TO  23 

IF  UOS.tQ.t)  GO  TO  27 
lOSsl 

REWIND  || 

27  ISTslST-l 

WRITE  (11)  1ST, (SRT(K) ,K=1 , 1ST) 

IC  >!T9»ICNT5M 
I S T » 1 
GO  TO  23 

28  CONTINUE 
I F L a 1 

29  CONTINUE 

IF  (INI ,ED,0)  GO  TO  30 
WRITE  (13)  (CURL(I),I=1,I2) 

JO  CONTINUE 

IF  (INT.NE.O)  REWIND  13 
C 

C SORT  AUL  COEFFICIENTS  ACCORDING  TO  RETARDED  TIME 

C 

ISwTtlSTJaO 

IF  (103. EO, 0)  GO  TO  31 

WRITE  (11)  1ST, (SRT(J),Jsl,IST) 

ICNTS3ICNTS*1 
REWIND  II 


CF  176 
CF  177 
CF  178 
CF  179 
CF  180 
CF  181 
CF  182 
CF  183 
CF  iea 
CF  185 
CF  186 
CF  187 
CF  188 
CF  189 
CF  190 
CF  191 
CF  192 
CF  193 
CF  1 9U 
CF  195 
CF  196 
CF  197 
CF  198 
CF  199 
CF  200 
CF  201 
CF  202 
CF  203 
CF  20a 
CF  205 
CF  206 
CF  207 
CF  208 
CF  209 
CF  210 
CF  211 
CF  212 
CF  213 
CF  21  a 
CF  215 
CF  216 
CF  217 
CF  218 
CF  219 
CF  220 
CF  221 
CF  222 
CF  223 
CF  ??U 

CF  225 
CF  22 6 
CF  227 
CP  228 
CF  229 
CP  230 
CF  231 
Cr  2 >2 
CP  233 
CP  2 Ja 
CF  235 
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31 


icr  = i 

DO  90  Llsl ,NT*AX 
UrNJf'AX-Lt 
ISTsl 

I SL»9999999 

1 F (IOS.EO.O)  CO  TO  3« 

ISLeO 

RE-INO  1J 
GO  TO  39 

32  IE  (•INOX/IOOOO.EG.L)  GO  TO  33 
!HCG  = 0 

GO  TO  3“ 

33  IHOGel 
ICOF(ICF)sINDX 
ICFsICF+l 

39  1 F (IST.LE.ISL)  GO  TO  35 

HEAO  (11)  ISL»CSHT(J),J=1,ISU) 

I STa  l 

35  INDxalsRTdST) 

ISTsISTM 

IF  (INOX)  3a, uo, 35 

36  KKsJNDX/l 0000 

IF  (IHOG.ES.0)  GO  TO  39 

1COF(ICF)sINOX 

ICFsICF'M 

KKPalST 

00  37  J=1,KK 

COF(ICF)sSRT(KKP) 

I C F = I C F ♦ 1 

37  KKPsKKP*! 

IF  (ICF.LT.ICMAX)  GO  TO  39 
IF  dOC.EQ.l)  GO  TO  38 
iOCal 

REWIND  12 
36  ICFalCF-l 

-HITE  (12)  ICF,U, (COF(J), J«l, ICF) 

I C TOT* 1 CTQT ♦ ICF 

1CNTC*ICNTCM 

ICFsl 

39  ISTsIST*KK 

GO  TO  39 
90  CONTINUE 

ICOF ( ICF )c0 

IF  (IOC.EO.O)  GO  TO  91 

-RITE  (12)  ICF.L, (COF(J),Jal,ICF) 

ICTOTUCTQTdCF 

1CNTC*1CNTC*1 

KE-Inu  13 

«l  IF  ( I 05  , F fj  , l ) Hf.JNn  11 

IF  (ISHT(NCG“X),NE, 9896969)  GO  TO  «2 
IF  ( ICOF (NCF-« ) ,Nt ,989*969)  GO  TO  9i 
IF  (ICMC.EG.O)  ICTuT  = lCF 
PRINT  99,  ICTCT, ICNTS, ICNTC 
RETURN 

92  PRINT  96 

STOP 

U3  PRINT  97 

S TOP 

99  PRINT  98,  J 

STD** 


CF  23b 
CF  237 
CF  238 
CF  239 
CF  200 
CF  291 
CF  292 
CF  293 
CF  209 
CF  295 
CF  29b 
CF  297 
CF  296 
CF  299 
CF  250 
CF  251 
CF  252 
CF  253 
CF  ?59 
CF  255 
CF  25b 
CF  257 
CF  258 
CF  259 
CF  260 
CF  26  1 
CF  262 
CF  263 
CF  269 
CF  2b5 
CF  266 
CF  267 
CF  268 
CF  269 
CF  270 
CF  271 
CF  272 
CF  273 
CF  279 
CF  275 
CF  276 
CF  277 
CF  27e 
CF  279 
CF  280 
CF  281 
CF  282 
CF  283 
CF  289 
CF  285 
CF  28b 
CF  287 
CF  258 
CF  289 
CF  290 
CF  291 
CF  292 
CF  293 
CF  299 
CF  295 
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f 


c 

af> 

<17 

«8 

a9 


FORMAT  ( 1 * » 27HHC  T A B 0 F 0 TIMES  FROM  SEGMENT , I S, I 2h  to  SEGMENTS, IS, 
1 ««  AND,I5,10M  0IFFE“  f?  Y » 1 5 , 1 1 H TIME  STEPS) 

FOP  mat  ( l y,2rhgvehf LO*  IN  FILLING  array  $9T> 

F OP  MAT  ( 1 X ,29HOVERFLO>*  IN  FILLING  ARRAY  COF  ) 

F OP M A T ( 1 X,29msEGm£NT  CONNECTION  ErrOR,  Js#I5) 

FORMAT  ( //,29h  LENGTH  OF  COEFFICIENT  A RR A Y 3 , I 1 0 , / , 1 X , I 1 0 , l X , 

l 2SHBLOCKS  Output  to  file  i l # » i x, i 1 o, ix , bphulocks  output  tu  file 
22,/) 

END 


CF  296 
CF  297 
CF  298 
CF  299 
CF  100 
CF  J01 
CF  i 02 
1 CF  JOJ 
CF  J0« 
CF  J05- 
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SUBROUTINE  CONECT 

CT 

I 

c 

CT 

2 

c 

connect  SETS  UP  segment  CONNECTION  Data  IN  ARRAYS  ICONI  and  ICON,? 

CT 

3 

c 

BY  SEARCHING  FOR  SEGMENT  ENDS 

THAT  ARE  IN  CONTACT, 

CT 

<1 

c 

CT 

S 

COMMON  /data/  N,np,X(200),Y(200).Z(200),S1(200),BI(20O),»LP(200), 

CT 

6 

1 HET(200), ICON  1(200), ICON ?(200),1TAG(200), IPX, IPY.IPZ 

CT 

7 

PlMfNSIUN  x2  < 1 ) , Y2 ( | ) , / 2 ( | ) 

CT 

8 

EQUIVALENCE  UZln.SKD),  ( Y2  (I  ) , ALP  ( I)  ) , ( Z 2 ( 1 ) , BE  T ( 1 ) ) 

CT 

9 

SMINsl ,E-3 

CT 

10 

JNOiO 

CT 

1 1 

DO  1 1=1, N 

CT 

12 

ICONI ( I )=0 

CT 

1 3 

J 

I CON2 ( I ) = 0 

CT 

1 9 

DO  33  1*1, N 

CT 

IS 

XI 1 = X C I ) 

CT 

16 

YI1=Y(I) 

CT 

1 7 

ZI1=Z(I) 

CT 

18 

X 1 2* X2 ( I ) 

CT 

1 9 

YI2=Y2(I) 

CT 

20 

ZI2  = Z2(  I ) 

CT 

21 

SLENsSURT((XI2*XI1)**2+(YI2-YI1)*a2aCZI2-ZI1)»»2) 

CT 

22 

c 

CT 

23 

c 

DETERMINE  CONNECTION  data  for 

END  1 OF  SEGMENT, 

CT 

29 

c 

CT 

25 

IE  (IPX. EC. 0)  GO  TO  3 

CT 

2b 

SEPsxIl/I.LfcN 

CT 

27 

SER2=XI2/SLEN 

CT 

28 

IF  (SEP.GT.SMlN)  GO  TO  2 

CT 

29 

IF  (SEP, LT. -Shin)  GO  TO  39 

CT 

30 

IF  (SER2.LE.SMIN)  GO  TO  39 

CT 

31 

ICON!  ( I )=200C!)EIPX 

CT 

32 

GO  TO  20 

CT 

33 

2 

IF  (SEP2.GT.SMIN)  GO  TO  3 

CT 

39 

IF  (SEP2.LT. "Shin)  GO  TO  39 

CT 

35 

ICUN2(I)=20000»IPX 

CT 

36 

GO  TU  7 

CT 

37 

i 

IF  (IPy.EQ.O)  GO  TO  S 

CT 

38 

SEPsV  I 1 /slln 

CT 

39 

S£P2*  Y I 2/SLEN 

CT 

90 

IF  (SEP.GT.SMIN)  GO  TO  9 

CT 

ui 

IF  (SEP.LT.-SMIN)  GO  TO  3<I 

CT 

92 

IF  (SEP2.LE.5MIN)  GO  TO  39 

c r 

«3 

ICONj (I)=?0000*IPY 

CT 

99 

rn  T n ■>  A 

tf  U 1 U L V 

CT 

95 

<1 

IF  (SEP2.GT ,S”IN)  GO  TO  5 

CT 

96 

IF  (SEP2.LT. -Smin)  GO  TO  39 

C T 

«7 

ICUN2(l)s20000«IPY 

CT 

98 

GO  TO  7 

CT 

99 

6 

IF  ( JP/.EQ.O)  GO  TO  7 

CT 

50 

SEPsZIl/SLEN 

CT 

51 

SEP2*ZI2/SLEN 

CT 

52 

IF  (SEP.GT.SMIN)  GO  TO  6 

CT 

53 

JF  (SEP.LT.-SMIN)  CO  TO  JU 

CT 

5u 

IF  (SEP2.LE.SMIN)  GO  TO  39 

CT 

55 

ICON) (I)=20000*IPZ 

CT 

56 

A- 16 


CO  TO  20 

CT 

57 

6 

1 F (StP2.CT.3MN)  GO  TO  7 

CT 

58 

IF  (SfcP2.LT, -SHIN)  GO  TO 

34 

CT 

54 

ICON2C  ) = ? 0 0 0 0 ♦ JPZ 

CT 

60 

7 

IF  ( ICON  1 ( I ) , NE , 0 ) GO  TO 

20 

CT 

61 

DO  4 1 C 3 1 , N 

CT 

62 

IF  (IC.EG.I)  GO  TO  4 

CT 

63 

IStGsic 

CT 

64 

IF  (ICClNi(IC).NE.O)  GO  TO 

S 

CT 

65 

SEPs(AoS(«Il-*(IC))»A8SCYll-V(IC))tAaS(ZIl-Z(IC)))/SLEN 

CT 

66 

IF  ($E*>.LT.SMJN)  go  TO  10 

CT 

67 

S 

IT  ( I Cu^'c  C X C ) • N'c  t 0 ) GO  TO 

4 

CT 

63 

SEP5(APS(xIl-x2(lC))*A8S(YIl-Y2(Icn*ABS(ZIl-Z2(IC)n/SLEN 

CT 

64 

IF  (SEP. LT, SHIN)  GO  TO  15 

CT 

70 

4 

CONTINUE 

CT 

71 

GO  TO  20 

CT 

72 

10 

IF  (ICON  1(1 SEG))  12,11,13 

CT 

73 

11 

ICON  1 ( I ) s I SfcG 

CT 

74 

ICON1 (ISEG)al 

CT 

75 

GO  TO  20 

CT 

76 

1? 

ICONi ( I ) = I CONI ( I SEG ) 

CT 

77 

GO  TO  20 

CT 

78 

13 

JNO=JNO»l 

CT 

74 

ICON  1(1) e JNO 

CT 

80 

1 X = I CON  1 ( I SE  G ) 

CT 

81 

ICONI ( I SEG ) 3 JNO 

CT 

82 

IF  (ICONi ( IX) ,EQ. ISEG)  GO 

TO  14 

CT 

83 

ICON2 (IX) sJNO 

CT 

84 

GO  TO  20 

CT 

8S 

1« 

ICONI (IX)bJNO 

CT 

86 

GO  TO  20 

CT 

67 

15 

IF  ( ICON2(ISEG)  ) 17,16,18 

CT 

88 

16 

ICON! ( I ) = I SEG 

CT 

84 

IC0n2(I5EG)3I 

* 

CT 

40 

GO  TO  20 

CT 

41 

17 

ICONI (I )=ICON2(ISEG) 

CT 

42 

GO  TO  20 

CT 

43 

18 

JNOsJNO-1 

CT 

4<I 

ICONi ( I )=JNO 

CT 

45 

IXsICON2( ISEG) 

CT 

46 

lCON2(ISEG)=JNO 

CT 

47 

IF  (ICnNl(IX),fcO,ISEG)  GO 

TO  14 

CT 

48 

ICON2 ( 1 X ) : JNO 

CT 

<J4 

GO  TO  20 

CT 

1 no 

14 

ICONI (IX)=JNO 

CT 

101 

C 

CT 

102 

C 

DETE»MNt  CONNECTION  DATA 

FOR  END  2 OF  SEGMENT, 

CT 

103 

C 

CT 

1 04 

20 

IF  ( ICONS (I).Nfc.O)  GO  TO 

33 

CT 

105 

DO  22  I C - 1 i N 

CT 

1 06 

IF  (IC.EO.I)  GO  TO  22 

CT 

107 

ISt&eIC 

CT 

108 

IF  (lCONl(IC).NE.O)  GO  TO 

21 

CT 

104 

SEPs(4rtS(xI2-x(IC))*AilS(YI2»Y(IC))**8S(ZI2-Z(IC)))/3LEN 

CT 

1 1 3 

IF  (SEP. LI ,SMIN)  Gu  TO  2) 

CT 

1 1 1 

21 

ix  ( I CQN2 ( I C ) , nx , 0 ) GO  TO 

22 

CT 

112 

SEP3(AaS(Xl2-X2(IC))AA05(Yl2-Y2(IC))AA9S(ZI2-Z2(IC)))/SLEN 

CT 

1 1 3 

IF  (SEP.LT , 3 4 I N ) GO  TO  28 

CT 

114 

22 

CONTINUE 

CT 

115 

GO  TO  33 

CT 

116 
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23 

IF  ( ICONl ( I3EC) ) 25,20,26 

CT 

1 1 7 

24 

ICON2( I )=ISEC 

CT 

118 

ICONl ( ] St  C ) * I 

CT 

119 

GO  TO  33 

CT 

120 

25 

lCON2(l)=ICONi (ISEGl 

Cl 

121 

GO  TO  33 

CT 

122 

26 

jNOsJNO-1 

CT 

123 

ICON2(I)=JNO 

CT 

124 

IX  = ICOM(lSEG) 

CT 

125 

I CUM  ( ISEG)=JNO 

CT 

125 

IF  (ICOMUX),EQ,ISEG)  CO 

TO 

27 

CT 

127 

1C0N2( I<)*JNO 

CT 

128 

GO  TO  33 

c r 

129 

?7 

ICONl (IX)sJNO 

CT 

l 30 

GO  TO  33 

CT 

131 

28 

IF  ( ICUN2  < I SEG ) ) 30,29,31 

CT 

1 32 

29 

I CON2 ( I ) s I SEG 

CT 

1 33 

I CON2 ( I SEG  ) a I 

CT 

134 

GO  TO  33 

CT 

1 35 

30 

ICON2(I)sICOn2(ISEG) 

CT 

136 

GO  TO  33 

CT 

137 

3 1 

J Nils  J NO*  1 

CT 

1 38 

1Cun2(I)sJN0 

CT 

1 39 

IX  = KOn2(1SEG) 

CT 

140 

1 CON 2 ( I SEG I s JNU 

CT 

141 

IF  (ICONl ( IX) ,t0, I SEG)  GO 

TO 

32 

CT 

142 

ICON2( IXJsJNO 

CT 

143 

GO  TO  33 

CT 

144 

32 

ICON  1 ( 1 X ) = JNQ 

CT 

145 

33 

CONTINUE 

CT 

146 

RETURN 

CT 

147 

34 

PRINT  35,  I 

CT 

143 

STOP 

CT 

149 

C 

CT 

150 

35 

FORMAT  (1X,7hSEG^ENT,IU,33M 

LIES  IN  OR  BEHIND  SYHHETRY  PLANE) 

CT 

151 

END 

CT 

152* 
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SUBROUTINE  CONSET 

CO 

1 

CO 

2 

CON'SET  COMPUTES  INTERPOLATION  CONSTANTS  FOR  USE  IN  EVALUATING 

CO 

3 

T«t  HELDS 

CO 

a 

CO 

5 

COMMON  /DATA/  N,NP,X(200)»Y(2CO)>Z(200),3I(200)«BI(200)«ALP(200), 

CO 

6 

1 BE T(200)  , ICON  1(200), ICON2C200) , ITaG(200),IPx,Ipy,IPZ 

CO 

7 

CORRON  /InTERP/  aT(3#200),9T(3,200)#CT(3.200)*E3(3).FSC3),GS(3), 

CO 

3 

t E ( 3 ) > h ( 3 ) 

CO 

9 

COMMON  /array/  GI1,QI2,QI3,EC(200,5),EQ(200#5),IPET(200).BFWUS400 

CO 

1 0 

1 ) 

CO 

1 1 

COMMON  /SCORP/  SxC200),SYC200),3Z(200) 

CO 

12 

COMRON  /CONST/  Ci>T,vEL,DT,NTSTEP 

CO 

1 3 

00  5 1=1 ,N 

CO 

la 

CALF=COS(ALP( I ) ) 

CO 

IS 

SX(n=CALP*CDSCnET(l)) 

CO 

16 

SY ( I )=CALF*SIN(8ET(I ) ) 

cu 

17 

szm  = siN(ALP(D) 

CO 

16 

1 M* ICON  1(1) 

CO 

19 

1P= ICON2 ( I ) 

CO 

20 

IF  (1M.GT. 19000  ) I R = I 

CO 

21 

IF  (IP, GT. 19000)  1P=I 

CO 

22 

S - 3 1 C I ) 

cu 

23 

IF  (IM.EQ.O)  GO  TO  l 

CO 

2a 

OM=,«5»(S»SI  C I”)  ) 

CO 

25 

GO  TO  2 

CO 

26 

DM= , 5*S 

CU 

27 

IF  (IP.EQ.O)  GU  TO  3 

CO 

28 

DP=.5*(StSI(IP)) 

CO 

29 

GO  TO  R 

CO 

30 

DP  = .5»3 

CO 

31 

D 1 =D  M * ( DmaOP ) 

CO 

32 

02=-OM«OP 

CO 

33 

03=0Po(DM*0P) 

CO 

3a 

A T C 1 , I)  = l , / D 1 

CO 

35 

AT (2, I )= 1 ,/D2 

CO 

36 

AT(3,I)=l./03 

CO 

37 

BT ( 1 i I ) s-OP/Ol 

CO 

36 

B T { 2 i I ) = (OM-OP)/D2 

CO 

39 

bT(3, I JSDR/D3 

CO 

U0 

CT ( 1 . I ) =0 , 

cu 

«1 

C 7 ( 2 # I )=-0P*DM/02 

CO 

U2 

CT( J, I ) =0 , 

cu 

a 3 

ES'l >=1 ,/C2.«0T»0T) 

CO 

aa 

ESC2)=-2.«f Sll  ) 

CO 

as 

E 3 ( 3 ) = E S ( 1 ) 

CO 

Ub 

F3(3)=.S/0T 

CO 

a 7 

FS( 1 )=-FS(3) 

CO 

a? 

F S ( 2 ) - 0 * 

CO 

aa 

GS(1 )»0. 

CO 

50 

G S ( 2 ) = 1 , 

CO 

51 

G S ( l ) » 0 . 

CO 

52 

C2=l ,/CVEL‘VEL) 

CO 

53 

E ( 1 )=ES( l )»C2 

CO 

5a 

E (2)=E3(?)»C2 

CO 

55 

t ( 3 ) =£S  C 3 ) *C2 

CO 

56 
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C?a-2,/VEl. 

H(n=Escn»ca 

M(3)=ES( 3)*C2 

OIl=((tS(i)*0T/J,*.5*FS(l))*0T*G3(i)).oT 
QIt'e(([_S(2)*OT/3,*,5«FS(£))«oT»G5(2))*DT 
OIi=((ES(})»OT/j,»,5*F3Ci))*OT*GS(l))*l)T 
HE  TUHN 
t NO 


CO  5 7 
CO  58 
CO  59 
CO  60 
CO  61 
CO  62 
CO  63 
CO  6« 
CO  65- 


6 
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SUBROUTINE  CONVRT 

CV 

1 

c 

CV 

2 

c 

CONVERT  Changes  GEOMETRY  04TA  from  The  form  stating  X,Y ,z 

CV 

1 

c 

COORDINATES  OE  EACH  SF&ment  ENO  TO  X,Y,Z  OF  the  SEGMENT  CEnT£R 

Cv 

a 

c 

PLUS  SEGMENT  LENGTH  and  ORIENTATION  ANGLES  AS  REQUIRED  IN  MAIN 

CV 

5 

c 

program, 

cv 

6 

c 

cv 

T 

COMMON  /OATA/  n,nP,X(200),Y(200),ZC20C).SI(200).&I(200),alPC200), 

Cv 

A 

1 BET (200) # ICunI (200) , 1CUN21200) , I TAGC200) , IPX, 1PY, IPZ 

cv 

9 

DIMENSION  X2  ( 1 ) , Y2(l),  Z2(l) 

CV 

J 0 

EQUIVALENCE  (X2(1),SI(D),  ( Y2  ( l ) , ALP  ( l ) ) , ( Z2 ( 1 ), BE T ( 1 ) ) 

cv 

1 1 

DO  1 JsJ.N 

cv 

12 

X A s X ( I ) 

cv 

1 5 

YA=Y ( I ) 

cv 

1 R 

Z A B Z ( I ) 

CV 

15 

X8=X2 ( I ) 

CV 

1 6 

YB  = Y2(  I ) 

cv 

17 

Z8sZ2( I ) 

cv 

18 

X(I)=(XA*XB)*,5 

Cv 

1R 

V ( I ) * ( V A * Yri  ) « , 5 

CV 

20 

Z(I)  = (ZA»ztJ)*,b 

cv 

21 

XAlXB»XA 

CV 

22 

YA=YB*YA 

cv 

2 3 

ZA=ZB-ZA 

cv 

24 

SI ( I ) =SQRT ( X A • X AtYA*YA*ZA*ZA ) 

cv 

25 

ALPU)  = »SIN(ZA/3I  ( I) ) 

Cv 

28 

J 

6ET(I)~-Arc.N2(VA,xA) 

CV 

27 

RETURN 

CV 

23 

END 

CV 

29 
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SUBROUTINE  OATAGN 

OA 

1 

c 

OA 

2 

c 

OATAGN  IS  THE  main  ROUTINE  FOH  INPUT  OF  GEOMETRY  OATA. 

OA 

3 

c 

Da 

4 

COMMON  /DATA/  N,NP,Xt200),Y(200),Z(200),SI(200),BI(200),ALP(200). 

DA 

5 

1 BET (200) , ICON  1 (200  ) , ICON  2(200  ) ,ITAG(200),IPX,IPY,IPZ 

OA 

6 

DIMENSION  x 2 ( 1 ) » Y2( l ) , 2 2(1) 

OA 

7 

DIMENSION  a T 3 T ( 4 ) 

D M 

S 

INTEGER  GM.ATST 

DA 

9 

EQUIVALENCE  ( *2( 1 ) , SI ( 1 ) ) , (Y2(l),ALP(i)),  ( 22  (1 ) , BE T ( 1 ) ) 

DA 

10 

Data  aTST/2mG*,  2HC,S,2HGE,2HGM/ 

DA 

1 1 

OATA  TA/, 017453292/ 

DA 

12 

N « 1 RE  s 0 

OA 

1 3 

N = 0 

DA 

1 4 

PRINT  8 

DA 

IS 

PRINT  9 

DA 

lb 

c 

DA 

1 7 

c 

READ  GEOMETRY  Data  Card  AND  ORanch  TO  SECTION  FOR  OPERATION 

DA 

18 

c 

REQUESTED 

DA 

1 9 

c 

DA 

20 

1 

READ  10,  GM,  I TG, ns, xwl , YM , Zwl , XW2, YW2, ?w2, RAD 

OA 

21 

IF  (GM.EQ,ATST( 1 ) ) GO  TO  2 

DA 

22 

IF  (GM,EQ.ATST(2) ) GO  TO  3 

DA 

2 3 

IF  (GM,EC1,ATST(3)  ) GO  TO  6 

OA 

24 

IF  (GM.EQ.aT3T (4) ) GO  TO  5 

DA 

25 

GO  TO  7 

DA 

2b 

c 

DA 

27 

c 

GENERATE  segment  Data  FOR  STRAIGHT  WIRE, 

DA 

28 

c 

Da 

29 

j 2 

NMIRE=nmIRc*1 

DA 

30 

1 

I l=Nt I 

DA 

31 

J2eN*NS 

Da 

32 

PRINT  11,  NRIRE,X^1,Ywi,Z-.1,Xw2,Yw2,2m2,RAD,nS,I1,I2,ITG 

Oa 

33 

CALL  WIRE  ( X*i  , Y„1  , 2*1 , X*2, Y*2, Zw2, RAD, NS, I TG) 

DA 

34 

GO  TO  1 

DA 

35 

c 

DA 

3b 

c 

scale  STRUCTURE  DIMENSIONS  BY  FACTOR  XWl, 

DA 

37 

c 

DA 

38 

3 

DO  u I * 1 , N 

OA 

39 

X(I)tT(I)»X*t 

OA 

40 

*(l)«»(I)»«»l 

OA 

4 1 

Z(I)sZ(I)*XRl 

Oa 

42 

X2(I)=x2(I)*X«l 

DA 

43 

Y2(I)SY2CI)*XWI 

DA 

44 

Z?(I3=Z2(I)*XM1 

DA 

45 

4 

erm«am>*x«u 

0 A 

4fe 

PRINT  12,  X A I 

DA 

47 

GO  TO  1 

DA 

48 

c 

0» 

UN 

c 

MOVE  STRUCTURE  OR  REPRODUCE  ORIGINAL  STRUCTURE  IN  NEh  POSITIONS, 

OA 

SO 

c 

OA 

51 

* 

PRINT  13,  ITG,NS,x»1,y«1,ZR1,x*2,vm2,Zw2,RAD 

ru 

52 

I«1»XWI*TA 

DA 

S3 

Y R 1 5YW1 • T A 

DA 

S 4 

ZwliZ-l *T» 

DA 

55 

CALL  MOVE  ( XRl  , Y-l  , Z*1  , x«,2,  y«2,  Z*2,  INT  (RADA.S)  , NS,  ITG) 

DA 

Sb 

L _ .. 
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9 


10 

11 

12 

11 

»« 

IS 


CO  TO  1 

TERMINATE  structure  geometry  INPUT, 

C *LL  CONECT 
CALL  CONVRT 
RETURN 
PRINT  i« 

PRINT  IS,  GM,ITG,N3,X-i,YRl,Z«l.XR2,YR2,Z*2.R*D 
STOP 

FORMAT  (////, 33X, 35H.  . - STRUCTURE  SPECIFICATION  - - «,//,37X, 

1 28mC00R0 I NATES  MUST  BE  INPUT  lN,/,37X, 

2 29HMf;rtRS  OR  BE  SCALtO  TO  METERS, /,37X, 

3 UhBEFORE  STRUCTURE  INPUT  IS  ENDED,  // ) 

FORMAT  ( 2x,«H*IrE , 79X,oHUO.  OF , « X , 5HF I R S T , 2 X , UHL A ST , 5 x , 3h T AG , 

1 /, 2X  * 3MNO. » 8X, 2«X 1 » 9X , 2mY l , 9X, 2MZ1 » 1 OX , 2Hx2*9x»  2MY2, 9X  » 2HZ2, 6X, 

2 6MRADiUS,3X,«HSECt,5X,UMSEG,,iX/i|HSEG,,5x,jMNO,) 

FORMAT  (42,11, IS, 7F10.S) 

FORM*!  (tX,IS,iFlltS,lX,<lFll.S,2X,IS,RX,IS,lX,lS,3XiIS) 

FORMAT  (6X.26RSTRUCTURE  SCAtEO  BY  F AC  TOR , F l 0 . 5 > 

FORMAT  (oX,«9MTHE  STHUCTORE  HAS  BEEN  MOVED,  MOvE  DATA  CARO  IS  -/ 
1 bX,13,15,7F10,5) 

FORMAT  (2SM  GEOMETRY  OATA  CARO  ERROR) 

FORMAT  <1X,A2,I3,I5,7F|0.5> 

END 


OA 

57 

DA 

58 

04 

59 

DA 

60 

04 

61 

DA 

62 

D4 

63 

DA 

6<J 

DA 

65 

04 

66 

DA 

67 

DA 

68 

OA 

69 

DA 

70 

04 

71 

DA 

72 

DA 

73 

OA 

7« 

DA 

75 

OA 

76 

DA 

77 

OA 

78 

DA 

79 

DA 

80 

OA 

81 

DA 

82- 
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SUBROUTINE  E1NCF  (T) 

El 

1 

c 

El 

2 

c 

EINCF  1$  CALLED  AT  EACH  TIME  STEP  TO  FILL  ARRAY  EINC 

hITH  Trt£ 

El 

3 

c 

APPLIEO  FIILO  ON  EACH  SEGMENT  AT  TIME  T,  THIS  ROUTINE  USES  A 

El 

4 

c 

GAUSSIAN  PULSE  OF  FORM  E XP ( » ( A * »2 ) • ( T -T M A X ) « »2 ) OR  GENERAL  TIME 

El 

5 

c 

OEPENOENCE  FROM  INPUT  TABLE) 

El 

6 

c 

El 

7 

COMMON  /DATA/  N,NP, X (200 )* V(200), 2(200), SI (200) *BI (200), ALP (200), 

El 

9 

I BE T (20  0), I COM (200  5 , ! COM2 ( 200 ) , I T AG ( 200 ) , IPX , IPY , 

IPZ 

El 

9 

COMMON  /SCOMP/  SX(200), 3^(200), SZ(200) 

El 

10 

COMMON  /CONST/  COT, VEL,DT,NTSTEP 

El 

tl 

COMMON  /EINC/  EINC(200),ESHC 

El 

12 

COmmun  /ASIJRC/  EMAG(  10  ) , ISRC  ( 10)  ,NSRC 

El 

13 

COMMON  /ENSET/  A, TMAX,EIX,PX,PY,PZ,tX,EY,EZ,VLT(200) 

,TR, IAR9, IAMX 

El 

1U 

DATA  TA/0. 017453292/, X*/0,1/ 

El 

15 

IF  (NSRC.EU.O)  GO  TO  6 

El 

16 

c 

El 

17 

c 

SET  UP  APPLIED  FIELD  FOR  TRANSMITTING  ANTENNA 

El 

18 

c 

El 

19 

DO  1 1=1, N 

El 

20 

1 

EINC ( I ) =0, 

El 

21 

IF  (IAR6.NE.0)  GO  TO  2 

El 

22 

eshc=a*(t-tmax> 

El 

23 

ESRC=EXP(-ESRC*ESRC) 

El 

24 

GO  TO  6 

El 

2S 

2 

TDsT 

El 

26 

3 

IF  (TD.LT.O.,OR,TO,GT,TAMX)  GO  TO  « 

El 

27 

IS*TD/Tw*j ,5 

El 

29 

IF  (IS.LT.2)  I S=2 

El 

29 

IF  (1S.E3.IARB)  IS»IS-1 

El 

30 

GO  TO  5 

El 

31 

a 

ESRC=0, 

El 

32 

RETURN 

El 

33 

5 

TObTD-TM«FLOaT( IS-I  ) 

El 

34 

ESRC«( ,5«TD»( ( TO-Tm)»VLT(IS-1 )»(TD*Tn)*VLT( IS»l ) )-(TD+Tw)«(T0-Tw)« 

El 

35 

IVLT(IS)  )/(T***T«) 

El 

36 

IF  (NSRC.EO.O)  RETURN 

El 

37 

6 

DO  T IS=),NSRC 

El 

33 

I > I SRC ( I S ) 

El 

39 

7 

EINC(I)=tMAG(lS)*ESRC 

El 

40 

RETURN 

El 

41 

C 

El 

42 

C 

SET  UP  APPLIED  FIELO  FOR  INCIDENT  PLANE  PAVE 

El 

43 

C 

El 

44 

8 

TCUncT-EIX 

El 

45 

DO  13  I «1 ,N 

tl 

40 

TDsTCON«(x(  I)*PX»Y(  I)«PY*Z(D»PZ) 

El 

«7 

IF  (JARB.Nt.O)  GO  TO  9 

El 

48 

TD»TD*  A 

El 

49 

ES"C«cX?(-TO*TO) 

El 

50 

GO  TO  12 

El 

51 

9 

IF  (TD,LT.O,.OR,TO,GT,TAMX)  60  TO  10 

El 

52 

I SsTn/T«*|  , 5 

FI 

53 

IF  (IS.LT.2)  IS«2 

El 

54 

IF  (IS.E0.1aR9)  I3«I3«1 

El 

55 

GO  TO  11 

El 

56 

A-24 


U I 

*** 

— 

to 

ElNC(I)aO, 

El 

57 

GO  TO  li 

El 

58 

It 

TOsTD«TR»FLOA  T ( I3»l ) 

El 

59 

ES«Ca{,5»T0*((T0-T*O»VLT(IS*n*{T0*Tw)«VLT(l9»|))-.(T0»TW)*CTD-Ta)* 

E I 

60 

IVLTC IS))/(T**T») 

ei 

61 

12 

EINCU)atS><C»($X<I)«EX»SY(l)*EY*SZ(l>*EZ) 

EI 

62 

1) 

CONTINUE 

EI 

63 

TDsTCON 

LI 

64 

1 F (lARQ.NE.O)  GO  TO  3 

EI 

65 

E3HCs»*TCnN 

EI 

66 

ESHCsEXPC-ESRCaESRC) 

EI 

67 

RETURN 

EI 

68 

ENTRY  ESET(T) 

EI 

69 

C 

EI 

70 

C 

READ  DATA  CARDS  AND  INITIALIZE  CONSTANTS  FOR  EXCITATION 

EI 

71 

C 

EI 

72 

READ  25#  IA«8,NSRC.Tw,TMET.PHl,ET.RZER0,TrtAX 

EI 

73 

PRINT  31,  lARb,NSRC,TN,THET,PM!,ET#RZERO»TMAX 

EI 

74 

PRINT  26 

EI 

75 

IF  (IAP9.NE.0)  GO  TO  15 

EI 

76 

ASSORT  (•AL0G(X#i)  )/(TR»,5) 

EI 

77 

IF  (TM*x,GT,t.E-20)  GO  TO  1« 

EI 

78 

E I NaO . 0 1 

EI 

79 

IF  (TM»x.LT.(-t,E-20))  EIN.-TMAX 

EI 

80 

TMAXsSORT(-*LOG(EIN))/A 

EI 

81 

THAXs(AlNT(TNAX/OT)*l ,)»DT 

EI 

82 

l* 

EIXsEXP(-A*A»TNAX«THAX)*100, 

EI 

83 

PRINT  27,  A,THAX,EIX 

EI 

84 

GO  TO  17 

EI 

85 

15 

IF  (IAOp.GT.200)  GO  TO  20 

EI 

86 

T AMXaTw»FLOAT l IARB-1) 

EI 

87 

read  21,  (VLT(I),Isl,lARB) 

EI 

88 

PRINT  23 

EI 

89 

Ta-T« 

El 

90 

DO  16  I ■ 1 , I AR8 

EI 

91 

TsTtTR 

EI 

92 

1 1* 

PRINT  22,  I,T,VLT(I) 

EI 

93 

* 

THAXaO. 

EI 

94 

17 

IF  CNSRC.EO.O)  GO  TO  19 

EI 

95 

PRINT  20 

EI 

96 

00  14  IS  = 1 « NSRC 

EI 

97 

READ  29,  1TG,II,VNAG 

EI 

98 

I si SEGNOC I TG, 1 1 ) 

EI 

99 

PRINT  30,  IS,ITG,II,I,VHAG 

EI 

100 

1 SRC ( 1 3 ) a I 

EI 

101 

ift 

E**AGCI9)aVRAG/SICI) 

EI 

102 

RETURN 

EI 

103 

l« 

PRINT  J2,  HZERO#TRET,PHT,ET 

EI 

104 

THuT=THt7*T« 

EI 

105 

PMlsPHI*TA 

EI 

106 

ETsET«TA 

EI 

107 

5T*SIN(THET) 

EI 

108 

C T sC'JS  ( T“ET) 

EI 

109 

SPsSIN(PhI ) 

EI 

110 

CPsCOS(PMl) 

El 

111 

SEsSIN(ET) 

EI 

112 

' | | 

CEsCOSCET) 

EI 

113 

PX«-ST«CP/vEL 

EI 

114 

PY:-3T»SP/VEL 

EI 

116 

5 

PZa-CT/VEL 

El 

116 

A- 2 5 

1 

EX*CT*CP*CE-SP*SE  El  117 

EY«CT*S»,»CLtCP*3E  El  118 

E7=-3T*CE  Cl  U9 

tlXiTMAXtRZERO/VEL  El  120 

RETURN  Cl  121 

20  PRINT  24  El  122 

STOP  El  125 

C El  124 

C El  125 

21  FORMAT  (6E12.5)  El  12b 

22  FORMAT  1 30  X , I5,E13.5,El4.S)  El  127 

23  FORMAT  ( J1X.52HTIME  OEPENCENCE  - GENERAL,  SET  BV  INPUT  TABLE  6 El  128 

1ELOR.//,31X.4M$TEP.7X,4HTImE,9x,  bHSOUKCE , / , 32X , 3HNO. , 6X , 6H ( SEC . ) , El  129 
2 7 X » 8HST  RENGT  M J El  130 

24  FORMAT  C52M  NUMBER  OF  EXCITATION  VALUES  EXCEEDS  ARRAY  OlMENSION)  El  131 

25  FORMAT  (2l5,6tl0,3)  El  132 

2b  FORMAT  ( //,38x,2bH-  ...  EXCITATION  ...  -,//)  El  JJJ 

27  FORMAT  ( 31X,7BHTIM£  DEPENDENCE-  GAUSSIAN,  PEAR  AT  T I ME  c T* AX  - - El  134 

1 EXP(-(a»«2)*(TIME-Tmax)*«2),//,31X,2has,E12,5,9M,  TMAXs,E12.5,  El  135 

2 //,31X, 10MSOURCE  H*S,F8,3,35M  PERCENT  OF  PEAR  VALUE  AT  TI"EsO.)  El  13b 

28  FORMAT  ( / / , 3 1 x , 17HVQLTAGE  SOURCES  • , // , 3 1 X , 2 3HSOURCE  SOURCE  SEG  El  137 
IMENT.bX, 12MPEAR  VOLTAGE, / , 3 5X , SHNO . , «X , JHT AG , 2X , 4H I NC , , 2X , 4HSE G , , El  138 

2 8X,7M(VOLTS),/,40X,3HnU.,9X,3hNO.)  El  139 

29  FORMAT  (215. CIO. 3)  El  140 

30  FORMAT  (SIX, 15, I7,2I6,5X,F10,4)  El  141 

31  FORMAT  ( //, 1 X , 26H*  * * EXCITATION  DATA  C ARD* * , 2 1 5, E 1 2 . 5 , «F 1 0 .5 , E 1 2 El  ia2 

l.S)  El  143 

32  FORMAT  ( //.31X.22HINCIDENT  PLANE  PULSE  31X.65HAT  TIMt  t II,  L El  144 

1EAOING  EDGE  OF  PULSE  IS  SHIFTED  BAC*  FROM  ORIGIN  BY.F10.5.7H  METER  El  145 
23,  //,  31 X, 26HINCIDENCE  ANGLES  - THE T A= . F 1 0 , 5 . 5M  DEG . , / , 5 1 X . 8MPH  I El  146 

3 *.Fl0,5,SM  OEG , , /, 5 1 X , bHE  T A s,F10,5,5H  DEG,,/)  El  147 

END  El  148* 
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SUBROUTINE  EM* T (Rx,FY,RZ,R2,T*U.SXJ,3YJ,3ZJ,J) 


SUBROUTINE  EM*  T (Rx , h Y , RZ , R2, T»U, 3X J , 3YJ , 3Z J , J ) 

EM 

I 

c 

EM 

2 

c 

EMxt  COMPUTES  The  COEFFICIENTS  WHICH  MULTIPLY  * 3 SEGMENT 

BY 

3 

EM 

3 

c 

T I ME  STEP  P*  TCH  OF  CURRENT  AnO  INTEGRAL  OF  CURRENT  TO  YIELD 

THE 

EH 

a 

c 

VECTOR  COMPONENTS  OF  ELECTRIC  FItLD 

EM 

! > 

c 

EM 

6 

COMMON  /DAT*/  N,NP,X(200),V(200),Z(200),3!(200),BI(200),ALP(2QO), 

EM 

7 

1 BET  (20  0),  I CONI  (200),  I COM2  (200),  lT*G(200),  IPX,  IPY,  IPZ 

EM 

0 

COMMON  /IntEhP/  XT(i,20u),Br(j,200),Ct(i,200),Es(3),Fs(i) 

, Gs  ( 3 ) , 

tH 

9 

1 E(3),m(3) 

EM 

10 

COMMON  /INTG/  XI , X2, X3» X4, X5, X6, X7, X8, X9, X 10, Xl  1 , X12 

EM 

1 1 

COMMON  /EMATS/  EMX(3,3),EMY(3,3),EmZ(3,3),OMX(3,3),OMY(3, 

3) 

, GHZ  ( 3 

EM 

12 

1,3) 

EM 

1 3 

COMMON  /CONST/  CDT,v£l,OT,NTSTEP 

EM 

14 

DIMENSION  P ( 3 ) , G ( 3 ) , M3) 

EM 

IS 

DO  1 L*l, 3 

EM 

16 

P(L)*2.*ES(L)*Tau*F3(L> 

EM 

17 

G(L)s(ES(L)*TAU*FS(L) >*TAUtGS(L) 

EM 

18 

1 

F (L)a-P(L)/YEL 

EM 

19 

BCON=-2,»tHX»SXJ»RV*SVJ*RZ»SZJ) 

EM 

20 

CCON=R2*OI(J)»BI(J) 

EM 

21 

SJ»S1 (J) 

Eh 

22 

CALL  INTEG  ( S J , 9C0N , CCON) 

EM 

23 

DO  2 Lsl,3 

EM 

24 

AL=AT(L, J) 

EM 

25 

HLsQT (L, J) 

EM 

26 

CL=CT(L, J) 

EM 

27 

DO  2 Ms), 3 

EM 

28 

EMsE(M) 

EM 

29 

FM=F(M) 

EM 

30 

GMaG(M) 

EM 

31 

TlsM(M)«(iL*XR*CL»Xl)*P(M)*(AL»XtO+BL»X6»CL*X2) 

EM 

32 

T2s2,*AL*(FM*X6+GH«X7)+fiL*(EM«xl+FH*x2*GH»x3) 

EM 

33 

T3a2,*AL*(EM*XR*FH*X10*GH*Xll )*BL*(FM«X<>4GM«X7) 

EM 

34 

TRs2,«*L*Ct*4*X6*FM*X7»GN*x8)*8L*(EI'oX2*FH*xJ*GM*X4) 

EH 

35 

T5s2.*AL*(tM*X10*FH*Xl 1*GM*X12)»6L*(EM«X6*FH*X7»GM«X8) 

EM 

36 

Tls-(Tl-VEL*T3)*t,E-7 

EM 

37 

T2»-.VtL*T2«l  ,E-7 

EM 

38 

c 

EM 

39 

c 

COEFFICIENTS  FUR  CURRENTS 

E M 

43 

c 

EM 

41 

EMX(L,h)bSXJ«TURX*T2 

EM 

42 

ENY(L,m)sstJ*TURY*T2 

EM 

43 

EMZ(L«M)ssZJ*T1*RZ«T2 

EM 

44 

T4s-vEL*VEL*T«*l .E-7 

EH 

45 

TSsVEL«YEL*T5*l,E-7 

EM 

46 

c 

EM 

47 

c 

COEFFICIENTS  FOR  TIME  INTEGRAL  OF  CURRENTS 

EM 

48 

c 

EM 

49 

Qmx(l«n)sPX»Tu*SxJ»T5 

EM 

50 

0Mr(L,M)aRT»Tu*STj*T5 

EM 

51 

2 

0mz(L,m)sRZ*T«*SZJ»T5 

EM 

52 

RETURN 

EM 

53 

end 

EM 

54. 

SUBROUTINE  F AC  10  (A,  aR,N,  IX, IP) 

FO 

l 

c 

FO 

2 

c 

8/P  NHICH  CONTROLS  I/O  FOP  FACTORIZATION 

FO 

3 

c 

FO 

4 

COMMON  ZmATPaR/  NBLCkS,NPBLK,NLAST,INT 

FO 

S 

DIMENSION  a(n,n),  JP(N),  ix(N),  ar(1) 

FO 

6 

IF  (NBL0KS.CT.2)  CO  TO  1 

FO 

7 

CAU.  F AC TR  (N,A,JP,N) 

FO 

B 

RETURN 

FO 

9 

1 

RE  « I NO  It 

FO 

10 

REWIND  13 

FO 

11 

RE w I NO  ta 

FO 

12 

ITbnPBIK*N 

FO 

13 

It»t 

FO 

14 

1 TwO*2 

FO 

15 

1 2s  I T 

FO 

16 

13*12*1 

FO 

17 

I4=2»IT 

FO 

IB 

c 

FO 

19 

c 

BUFFER  In  BLOCK!  AND  BL0CK2  F RQM  TAPE  1 

FO 

20 

c 

F(J 

21 

READ  (13)  (AR(I),IsIl,l2) 

FU 

22 

READ  (13)  (AR(1),IeI3,I4) 

FO 

23 

CALL  LFACTfi  (A,N,1,2,IP) 

FO 

24 

c 

FO 

25 

c 

BUFFER  OUT  BlOCKI  TO  TAPE2  (BLOCK!  FACTORED) 

FO 

26 

c 

FO 

27 

WRITE  (11)  (AR(I), 1*11,12) 

FO 

2d 

c 

FO 

29 

c 

BUFFER  OUT  BL0CK2  TO  F ILE3 

FO 

30 

c 

FO 

31 

WRITE  (1«)  (AR  (I), 1*13,10) 

FO 

32 

DO  2 1 X9LK2  = 3, NBLOKS 

FO 

33 

c 

FO 

34 

c 

BUFFER  In  B10CK2  FROM  TAPEl 

FO 

35 

c 

FO 

36 

READ  (13)  (AR(l),IaI3,14) 

FO 

37 

CALL.  LFACTR  ( A, N, 1 , I XBLK2, IP) 

FO 

38 

c 

FO 

39 

c 

BUFFER  OUT  BLOCK?  TO  F ILE3 

FO 

40 

c 

FO 

41 

WRITE  (1«)  (ARID, 1*13, 14) 

FO 

42 

2 

CONTINUE 

FU 

43 

IXHLKlal 

FO 

44 

3 

1 XBLK 1*1 XBLK 1 ♦ 1 

FU 

45 

I X8LK 2s  I XBLK  1 ♦ 1 

FO 

46 

c 

FO 

47 

c 

WITH  THE  EXCEPTION  OF  THE  FIRST  PASS,  IFILES  BECOMES 

I F ILE  4 AND  VI  FO 

48 

c 

FO 

49 

IFILE3*13 

FO 

50 

IFlLEwsl* 

FU 

51 

IF  (2*(IXBLK1/2).NE.IXBLK1)  GO  TO  4 

FO 

52 

IFILE3S14 

FO 

53 

IFILE4.13 

FO 

54 

4 

REwINO  IFILES 

FlJ 

55 

REMIND  IFILE4 

eU 

56 

A- 2 8 


c 

FO 

57 

c 

BUFFER  IN  0LOCM  * NO  BLOCKS  FROM  IFILE3  FO 

58 

c 

FO 

59 

REaO  ( IF  I LC 3 ) (AR(l),I«Il,I2) 

FO 

60 

RE *0  (IFILE3J  (AR(I),1*I3,I«) 

FO 

61 

CALL  LFACTR  ( A , N,  I XBL* l , I XBLK2 , IP) 

FO 

62 

c 

FO 

63 

c 

buffer  our  block i to  tares  (block t 

factored)  fo 

69 

c 

FO 

65 

WRITE  (11)  (AR(I),IeIliI2> 

FO 

66 

c 

FO 

67 

c 

BUFFER  OUT  BLOCKS  TO  If ILEtt 

Ft) 

68 

c 

FO 

69 

WRITE  (IFILER)  (AR< I ) , 1*13, 10) 

FO 

70 

IF  (IXBLK2.NE.NBLOKS)  GO  TO  S 

FO 

71 

c 

FO 

72 

c 

BUFFER  OUT  BLOCK?  TO  TAPES  (BLOCK! 

FACTOHEO—FACTORIZATION  FIUISHE  fo 

73 

c 

FO 

79 

WRITE  (ID  CARCI),I*I3,I«) 

FU 

75 

REWIND  11 

FO 

76 

REWIND  13 

FO 

77 

REwINO  1« 

FO 

76 

GO  TO  6 

FO 

79 

5 

IXBLK2=IXBLK2+I 

FO 

SO 

IF  (IX3LKS.GT.NBL0KS)  GO  TO  3 

FO 

ei 

c 

FO 

82 

c 

BUFFER  IN  BLOCKS  FROM  IFILE3 

FO 

83 

c 

FO 

eo 

READ  (IFILE3)  (AR(I),IaI3,I«) 

FO 

85 

CALL  LFACTR  (A,N, IXBLM, IX0LKS, IP) 

FO 

86 

c 

FO 

87 

c 

BUFFER  OUT  BLOCK2  TO  FREQ 

FO 

88 

c 

FO 

89 

WRITE  (IFILER)  (AR(I), 1*13, I«) 

FO 

90 

GO  TO  5 

FO 

91 

6 

call  LUNSCR  (A,AR,N,IX,IP) 

FO 

92 

RETURN 

FO 

93 

END 

FO 

99- 

SUBROUTINE  FaCTR  ( N . A . P . NO  I H ) 

FA 

1 

c 

FA 

2 

c 

SUBROUTINE  TO  FACTOR  A MATRIX  INTO  * UNIT  LO*ER  TRIANGULAR  MATRIX 

FA 

3 

c 

ANO  AN  UPPER  TRIANGULAR  matrix  USING  THE  GAUSS-OOOL I TTLf 

ALGuRI THM 

FA 

<1 

c 

PRESENTEO  ON  PAGES  UU-oib  OF  A,  RALSTON--A  FIRST  COURSE 

IN 

FA 

5 

c 

NUMERICAL  ANALYSIS,  COMMENTS  BELO*  REFER  TO  COMMENTS  IN 

RALSTONS 

FA 

b 

c 

TEST,  (MATRIX  TRANSPOSED) 

FA 

7 

c 

FA 

8 

OlMENSIUN  A(NDIM,NDIM),  P(NOIM) 

Fa 

9 

COMMON  /SCRATM/  0(200) 

Fa 

10 

INTEGER  H,P,HMl,RPt,PJ,PR 

FA 

11 

IFLGeO 

FA 

12 

DO  9 Rs 1 , N 

FA 

13 

c 

fa 

1 4 

c 

STEP  l 

FA 

15 

c 

FA 

lb 

DO  t X=1,N 

FA 

17 

0(K)SA(R,K) 

FA 

18 

1 

continue 

FA 

19 

c 

FA 

20 

c 

STEPS  2 ANO  3 

FA 

21 

c 

FA 

22 

RMJ«B-1 

Fa 

23 

IF  (RMl.LT.l)  CO  TO  « 

FA 

2a 

DO  3 Jcl(RMl 

FA 

25 

PJsP(J) 

FA 

2b 

A(R, J)sO(PJ) 

FA 

27 

DCPJ)sO(J) 

FA 

23 

JPJsJAI 

FA 

29 

DO  2 I s JP1 , N 

FA 

30 

DCI)6D(I)-A(J,X)»A(R,J) 

FA 

31 

2 

CONTINUE 

FA 

32 

J 

CONTINUE 

FA 

33 

a 

CONTINUE 

FA 

3a 

c 

FA 

35 

c 

STEP  0 

FA 

3b 

c 

FA 

37 

0MAXSAPS(0(R)) 

FA 

38 

P ( K ) »R 

FA 

39 

RP 1 sR 1 1 

FA 

ao 

IF  (RPJ.GT.N)  GO  TO  b 

FA 

«l 

DO  S I *RP  I * N 

FA 

02 

ELM AG*  * OS ( 0 ( I ) ) 

FA 

US 

IF  (ELMAG.LT.DMAX)  GO  TO  5 

FA 

aa 

OMax«ELMAG 

FA 

as 

P (R) • I 

FA 

Ob 

s 

CONTINUE 

FA 

07 

b 

CONTINUE 

FA 

as 

IF  (DMAX.LT. I. E-10)  IFLGaJ 

FA 

09 

PRaP(R) 

FA 

50 

A'  ( R . R ) aO  ( PR  ) 

FA 

51 

D(PR)aO(H) 

fa 

52 

c 

y a 

53 

c 

STEP  S 

FA 

50 

c 

FA 

55 

IF  (RPl.GT.N)  CO  TO  8 

FA 

5b 

A- 30 


oo  r i«Rpi»N 

A<R,l)«0<I)/A<R,K) 

7 continue 

6 CONTINUE 

IF  <IFLO.EO,0>  CO  TO  9 
PRINT  10,  R , DMA  X 
IFLCsO 

0 CONTINUE 

RETURN 
C 

10  FORMAT  (1X,6HPIV0T<»  I3.2H)s»EU.8> 

ENO 


FA 

FA 

FA 

FA 

FA 

FA 

FA 

FA 

FA 

FA 

FA 

FA 


57 

50 

59 

60 
61 
62 
63 
6U 

65 

66 
67 
68" 
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SUBROUTINE  Fame*  <NelOKS.NPBLK,NL*ST, IMAX,N, IhT) 

fb 

i 

F5 

2 

F BLOCK  DET{R«INES  BLOCK  SIZE 

AMO  NUMBER  of  blocks  when 

OUT-OF-CORE 

FB 

J 

MAfRJX  S10RAGE  15  RELUlRtt), 

InTso  PETLRNED  if  matrix 

FITS  IN  CORE 

F 9 

a 

FB 

5 

IF  (kan.LE. I^FA)  CO  TO  1 

FB 

b 

nP&L**INax/(2»n) 

fb 

T 

IF  (“PbLK.LT.l)  STOP 

FB 

B 

INTS1 

FU 

9 

nblo*s=<n-i )/np&lk+i 

fb 

10 

nL*STsn»(mblOKS-1)*AiP9LK 

fb 

1 I 

print  z,  NtlLONS.NPSLK.KLAST 

F8 

12 

RETURN 

FB 

11 

NBLOkS=2 

FB 

lu 

NPbLKS(NA( )/2 

FB 

15 

NLASTsN-NPBLK 

fb 

lb 

INT  = 0 

FB 

17 

RETURN 

FB 

18 

FB 

10 

FORMAT  (KiIIH  BLOCKING  * « 15/ ) 

FB 

20 

END 

FB 

21 
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DUCT 

Jlj  i 


t 

i 

i 


SUBROUT  I Nf  FURT  ( A , M , s , IFS#  I F £ R R ) 

FOURIER  TRANSFORM  SUBROUTINE.  PROGW*MMfo  IN  SYSTEH/360, 

BASIC  PROGRAMMING  SUPPORT,  FORTRAN  Iv.  F URN  C2B-6504 
TmJs  DEC*  SET  UP  FOR  1BSYS  ON  IBM  7094, 

DOES  EITHER  FOURIER  S YNTHE  S I S , I . £ . , COMPUTE  S COm“LEX  FUuPIER  SERIES 
G j v £ “■  i vfCTOR  OF  N C1’mPLEX  FOURIER  AhPLI  TL'OESiQR,  0!VE,J  A VECTOR 
OF  COMPLEX  U»It  % DOES  F liiiR  I E **  ANALYSIS,  COMPUTING  AMPLITUDES, 

A IS  A COMPLEX  VECTOR  OF  LENGTH  Ns****  CORPLtX  N'CS.  OR  2»M  HEAL 
NUMBERS.  A IS  TO  BE  SET  t)Y  USER. 
m 13  an  INTEGER  O.LT.m.LE, 1 3,  SET  Br  USER, 

3 IS  * VECTOR  S(J)=  SIn(2«PI«J/NP  ),  J=  1 NP /U  ■>  1 , 

C 0 " P U T E 0 BY  PROGRAM. 

IFS  IS  a PA«*ntTER  TO  BE  SET  H Y USER  AS  FOLLUR3- 
I F S = G TO  SET  NPs2**m  anD  SET  Go  SINE  TABLE. 

IFS=1  TO  SET  nsnps?**m,  SET  UP  SIN  TABLE,  AND  DO  FOURIER 

synthesis.  Replacing  the  vector  a by 

X(J)  = sum  over  KSC.N-J  OF  ACO*EXP(2*PIM/N)a*CJ*K), 

JSO.N-I,  WHERE  I=SORT(-J) 

TMt  X'S  A R a STORED  WITH  RE  X(J)  IN  CELL  2*J*1 
AND  IM  X(J)  IN  CELL  2*J*2  FOR  J = 0 , 1 , 2 , . , . , N- 1 , 

THE  A • S ARE  STOREO  In  T-E  S*mE  hannER, 

I F S = - 1 TO  3ET  NsNP  = 2»*H, SF T UP  SIN  TABLE,  AND  00  FOURIER 

ANALYSIS,  taking  the  input  VECTOR  A AS  X and 

REPLACING  IT  BY  THE  A SATISFYING  THt  ABOVE  FOURIER  SERIES, 

IFS=t2  TO  00  FOURIER  SYNTHESIS  ONLY,  WITH  A PRE-COmPUTED  S. 

IFS=-2  TO  DO  FuuRIER  ANALYSIS  ONLY,  WITH  A PRE-COMPUTED  S. 

IFESR  IS  SET  BY  PROGRAM  TQ- 
eO  IF  NO  ERROR  DETECTED. 

*1  IF  M IS  OUT  OF  range.,  UR,  "MEN  IFS=t2,-2,  THE 

pRf -COmPuTE 0 S Table  IS  NOT  large  ENOUGH, 

s-J  WHEN  IFS  *+1,-1,  MEANS  ONE  IS  RECOMPUTING  S TABLE 

UNNECESSARILY, 

NOTE-  AS  STATED  ABOVE,  THE  MAXIMUM  VALUE  OF  M FOR  THIS  PROGRAH 
ON  The  IBM  7094  IS  13,  FOR  360  MACHINES  HAVING  GREATER  STORAGE 
CAPACITY,  ONE  may  INCWIASE  This  LIMIT  BY  REPLACING  IS  IN 
STATEMENT  3 BELOW  mY  L0G2  N,  WHERE  N IS  THE  max,  NO,  OF 
COMPLEX  N U M 3 E R S 0K'E  Can  STORE  IN  hIGm-SPEEo  CORE,  ONE  MUST 
ALSO  *DD  -L'R£  00  STATEMENTS  TO  ThE  binary  SORT  ROUTINE 

F Gl  L U- i i,  STATEMENT  cm  AND  ChANGE  THE  EQUIVALENCE  STATEMENTS 
FOR  THE  K'S. 

DIMENSION  a(1),  SO),  k ( l u ) 

EGuI.-lEn C v (K(13),Ni)»  ( K ( t 21 , N2 ) , (N(J1  ),K3),  (K(10),Kr) 

EG'JI  V *'_£  'CE  (<(«>, «5),  («(B),wfc),  (<(Y),AY),  (K(b),K«) 

EOuIvalEnCE  («(5),kR),  (K(u),<lO),  (* ( 3) ,*  1 1 ) , (k(2),R12) 
EQUIVALENCE  (K(I).kIJ),  ( K ( l ) , n2 ) 

IF  ( m ) a, *,l 

IF  (m-U) 

IFERwoi 
RETURN 


\ 


FT 

1 

FT 

2 

FT 

3 

FT 

a 

FT 

5 

ft 

6 

FT 

7 

FT 

S 

FT 

9 

FT 

10 

FT 

1 1 

FT 

12 

FT 

13 

FT 

14 

FT 

IS 

FT 

16 

FT 

17 

FT 

18 

FT 

19 

FT 

20 

FT 

21 

FT 

22 

FT 

23 

FT 

20 

FT 

25 

FT 

26 

FT 

?7 

FT 

28 

FT 

29 

FT 

30 

FT 

31 

FT 

32 

FT 

33 

FT 

34 

FT 

35 

FT 

36 

FT 

37 

FT 

38 

FT 

39 

FT 

uo 

FT 

41 

FT 

<12 

FT 

43 

FT 

44 

F T 

45 

FT 

46 

FT 

47 

FT 

43 

FT 

49 

FT 

50 

FT 

51 

FT 

52 

FT 

53 

FT 

54 

FT 

53 

FT 

56 
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one 


a 


c 

c 

s 

6 

c 

7 

8 


10 


11 

c 

12 


1 J 

C 

14 


f 

15 

c 

16 

c 

c 

c 


IFFRR=0 

FT 

57 

N*2«»*« 

F T 

58 

IF  (I4C'SUcS)-n  25 » 25  * 5 

FT 

59 

wE  ARE  DOING  TnansfoR*  ONLY,  SEE  IF 

PUF. COMPUTED 

ft 

60 

S TABLE  IS  SUFFICIENTLY  LARGE 

FT 

61 

IF  <n-nP)  7,7, to 

FT 

62 

IFf RRal 

FT 

6 J 

GU  TO  25 

FT 

64 

SCRAkrtE  A , 6Y  SAnOE'S  METHOD 

FT 

65 

Ml  )=2*N 

FT 

66 

DO  d L=2,h 

FT 

67 

MD=ML-l)/2 

FT 

66 

DO  9 L=M,12 

FT 

69 

Ml*  1 ) a2 

FT 

70 

NOTE  EQUIVALENCE  OF  KL  AnD  M14-L) 

FT 

71 

to  I N A k Y SURT- 

FT 

72 

I J 32 

FT 

71 

on  11  J|s2,Kl,2 

FT 

74 

DO  11  J23J1.52.KI 

FT 

75 

DO  11  J4=J2.k3,K2 

FT 

7b 

DO  11  J4=J3,n4,k3 

FT 

77 

DO  11  J53J4,N5.*4 

FT 

78 

DO  11  JtosJ5,Kto,K5 

FT 

79 

DO  1 1 J7s.16,K7,K6 

FT 

80 

DO  11  Jd  = J7 , K 8 , < 7 

FT 

81 

DO  11  J9=J8,K9,K8 

FT 

82 

DO  11  J10  = J9,M0,K9 

FT 

83 

DO  11  JtlsJlO.Ml.MO 

FT 

84 

DO  11  J12sJll.M2.Kll 

FT 

85 

DO  11  Jl=Jl2,M3,M2 

FT 

86 

IF  (IJ-JI)  ICrll.ll 

FT 

87 

TsAUJ-U 

FT 

88 

A( JJ-1 )=4(JI-t) 

FT 

89 

ACJI-1 )=T 

FT 

90 

T s A ( I J ) 

FT 

91 

AtI,')  = A(JI) 

FT 

92 

A ( J I ) s T 

FT 

93 

IJ=I  J*2 

FT 

94 

IF  ( I F S ) 12.2.14 

FT 

95 

DOING  FOURIER  ANALYSIS, SO  DIV.  BY  N 

AND  CONJUGATE, 

FT 

9b 

F N 3 N 

FT 

97 

DO  13  I 3 1 , N 

FT 

98 

A(2»I-1)  = A(2«I-1  )/FN 

FT 

99 

a(2«I)3-A(2»I)/FN 

FT 

100 

SPtCIAU  CASE-  L=1 

FT 

101 

DO  IS  1 = 1, N, 2 

FT 

102 

TsA (2* I -1  ) 

FT 

103 

A ( 2 » I - 1 )sT*A(2*I*l  ) 

FT 

104 

A { 2*!*S)*!- A (2*1*1) 

FT 

l 05 

T s A ( 2 » I ) 

FT 

10b 

A(2*I)sT*A(2*I*2) 

FT 

107 

* ( 2* 1 ♦ ?) » T-A (2*  I ♦ 2) 

FT 

108 

IF  («-1 ) ?,4,1h 

FT 

109 

SET  FOR  L*2 

F T 

110 

L E * B 1 =? 

FT 

1 1 1 

L E *p  1 s2**(L-l ) 

FT 

112 

Lt  »p  = « 

FT 

1 1 3 

LExP«2»*(t*i) 

FT 

114 

NPL=2«*“T 

FT 

115 

NPL  3 nP*  2*a-L 

FT 

116 

A-  34 


DO  22  1 = 2, * 

FT 

117 

c 

SPECIAL  CASE-  J-0 

FT 

118 

DO  17  I=2,N2,LtXP 

ft 

119 

n = i *le  x»»  j 

FT 

120 

12=11 *LEYP1 

FT 

121 

1 3=I2*LEXP1 

FT 

122 

T*A( 1-1 ) 

FT 

123 

A(1-1)=Tta(12-1> 

FT 

129 

A(I2-t )=T-A(I2-1 ) 

FT 

1 25 

T = A ( I ) 

FT 

12ft 

A(i)=r*»(i2) 

FT 

127 

* ( 1 2 ) = T"* ( 12 ) 

FT 

128 

T=-A(I3) 

FT 

129 

TI  = 4( I 3-1 ) 

FT 

130 

ACI3-I )=A(I1-1)-T 

FT 

131 

A(li)=A(Il)-TI 

FT 

1 32 

A ( 1 1 - 1 ) = MIM)»T 

FT 

133 

17 

A(I1)=»(11 ) ♦ T I 

FT 

139 

IF  (L-2)  21,21,16 

FT 

13S 

18 

KLAST=n2-LEXP 

FT 

13  6 

JJ=nPL 

FT 

137 

DO  20  J=«,LEXP1,2 

FT 

136 

NPJJ=NT-JJ 

FT 

1 39 

UP»S(NPJJ) 

FT 

190 

UlrSCJJ) 

FT 

1 « 1 

IL*ST=J4KLAST 

ft 

m2 

DO  IP  1=  J,lL*ST,LEXP 

FT 

m3 

I 1 = 1 *Lt  XP  1 

FT 

199 

12=1 l»LEXPt 

FT 

195 

I3=I2+LfcXPl 

FT 

196 

T»A{ J2-1 ) «uS-A ( I 2) «UI 

FT 

my 

TI=»(J2-))*UI*A(12)*UR 

FT 

ms 

A(I2-1)=A(1-1)-T 

FT 

199 

A(12)=A(1)-TI 

FT 

150 

A(I-l)  = A(i-n»T 

FT 

151 

AtI)=AC)*TI 

FT 

152 

T*-A ( 13-1 ) *U1 -A ( I 3) »US 

FT 

153 

TI  = A(i  3-n»UR-A(l3)»ui 

FT 

159 

ACI3-1 )=ACI1-1)«T 

FT 

155 

A(I3)=A(11)-TI 

FT 

156 

A (I 1-1 )=A (11-1 )♦! 

FT 

157 

19 

A (1  ’.  ) = A U 1 ) ♦ T I 

FT 

156 

C 

END  OF  I LOOP 

FT 

159 

20 

JJsJJfMPi. 

FT 

160 

c 

END  OF  J LOOP 

FT 

161 

21 

LE  XP 1 =2*lE  XP1 

FT 

162 

LExP=2»LExP 

FT 

163 

22 

nPl=NPL/2 

FT 

169 

C 

FT 

165 

IF  ( I F s ) 2 J , 2 , 3 

FT 

166 

C 

DOING  F C'JP I fcP  ANALYSIS,  REPLACE  a BY  CONJUGATE, 

FT 

167 

23 

DO  29  I = l , N 

FT 

166 

29 

A (2« I ) s»  a (2* I ) 

FT 

169 

CO  TO  3 

FT 

170 

C 

hETUHN 

FT 

171 

c 

MAKE  TABLE  OF  S(J)=SIN(2*PI«J/NP)>Jel,2,,,,,NT-l,NTaNP/U 

FT 

172 

25 

NP«N 

FT 

173 

MpaM 

FT 

179 

NHN/9 

FT 

175 

MTfN-2 

FT 

176 
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If  CT)  31,31,26 

FT 

177 

26 

THETAs. 785398163 

FT 

178 

C 

THET*sBl/2«»(U*»)  FOR  1*1 

FT 

179 

JSTERsMT 

FT 

180 

C 

JSTEP  s 2«« ( HT-L+1  ) FOR  1*1 

FT 

181 

JOIFs.VT/2 

FT 

182 

C 

JO  I f = 2 • * ( h T *L ) FOW  L«t 

FT 

183 

SC  JOIF  )=S1M  THET*) 

FT 

189 

IF  (MT-2)  31,27,27 

FT 

185 

27 

DO  30  1*2, *T 

FT 

186 

THET1sthEt*/2, 

FT 

187 

jSTEkVsJSTER 

FT 

188 

JSTEPsJOIF 

FT 

189 

JPIFrJOlF/2 

FT 

190 

SCJOIE)=SIN(TH£TA) 

FT 

191 

JC 1 *M  »JDIF 

FT 

192 

SCJC1 JsCOSCTmETa) 

FT 

193 

JUASTSNT-JSTER2 

FT 

199 

IF  ( JlaST-JSTEP)  30,28,28 

FT 

195 

28 

00  29  JSJSTEP, JlAST, JSTEP 

FT 

198 

JCsNT-J 

FT 

197 

JOsJtJOIF 

FT 

198 

29 

s(jo)=scj)*scJcntscjoiFj«s(JC) 

FT 

199 

30 

CONTINUE 

FT 

200 

31 

IF  llfS)  7,3,7 

FT 

201 

END 

FT 

202 

* 
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SUBROUTINE  IK  EM  U) 

IF 

t 

c 

IF 

2 

c 

IFtM  CONVERTS  I FROM  E,  M,  OR  BLANK  TO  l,  *»,  OR  0, 

IF 

i 

c 

IF 

« 

OAT*  IE,  IH/1HE, 1HH/ 

IF 

S 

IF  CI.tG.lt)  CU  TO  I 

IK 

b 

IF  CI.EQ.IM)  CO  TO  2 

IF 

7 

1*0 

IF 

8 

RETURN 

IF 

R 

1 

I»1 

IF 

10 

RETURN 

IF 

1 1 

2 

I»-l 

IF 

12 

RETURN 

IF 

1) 

END 

IF 

t« 
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SUflHOUTINg  InTEG  (EL, 0,0 


C 

c 

c 

c 


I 


2 

i 


a 

5 


6 

7 


JNTEG  COMPUTES  12  INTEGRALS  NEEDED  To  SET  UR  THE  FIELD 
COEFFICIENTS 

Common  /InTG/  x 1 , Xl,  X3,  X'i , *5,  xb,  *7,  x9,  X9,  xl  0,  XU  , XI  2 
OOJHLE  PRECISION  52,Sl.PlS,W2S,Rl,R2,82,H9,0El,0E2.DIS,ALR,DlH, 
1 Of  t 9 , 0E2  S « Tl , T2, 02, 03,04,06 
S2=EL*.5 
31=-S2 

R 1 b= ( SJ  *m ) • S 1 *C 
P2S=CS2*8)*S2*C 
RlsOSORr(PiS) 

R2sDS0RT(R2S) 

02  = 8*0 

H9s , 5*6 

UElsSWHB 

OE 2=S2*M0 

l'lS*u,*C-02 

ALRsDL0G(R2/Rl) 

L IM  = 0 

IF  (AoSU>)  .LT.  1 .E-27)  GO  TO  2 
XLIm=DIS/AWS(0) 

IF  I»LlM.GT.l.t-S)  GO  TO  2 

L I ms  1 

01  ms  0,0*0 

IF  (OEl*DE2.GT,0.)  GO  TO  1 
PRINT  10,  e,c 

STOP 

OE 1 S = DE 1 *0E  1 
OE2S=DE2*Ot2 
GO  TO  3 

0 1 msDSGRT ( D I S ) 

*1=El 

IF  (OEl.LT.O.)  GO  TO  4 

02  = CLOG ( (R2*DE21/(RI*0E1  ) ) 

CO  TO  5 

D2=OtOG( (»1-Dtj )/(K2-DE2>) 

X2*D2 

IF  (LIM.EQ.1)  GO  TO  6 
T2=DATAS12,*0E2/DIM) 

T1=0*T*N(2. *DE1/DIM) 
ni=2,*(T2-Ti )/oik 
D«=a.* (DF2/R2-DE1 /R1J/0IS 
GO  TO  7 

DJ=(DIS/Cl2.*DE2«r>E2)-t ,)/0E2-CDIS/(12,*DEl*DEl)-l.)/DEl 
Oi3,5«DA0S{ (0E2/CE l -OE 1 /DE2 ) / ( P 1 *P2) ) 

*3=03 
xa=aa 
*5  = 0. 

Xt>  sH2-al-M0*D2 
X7=4LR-M3*0J 

if  (LiM.ro.n  go  to  0 

(,Es.2,«((M«S2*2,»C)/M2>(o*St»2,*C)/Rt)/0IS 
GO  Tu  V 

09*,25»(B/0E2S-9/0E 1$)*1 ,/OEl-l ./DE2 


IN  I 
IN  2 
IN  i 
In  4 
IN  5 
IN  6 
IN  7 
IN  6 
IN  9 
IN  tO 
IN  11 
IN  12 
IN  1 J 
IN  1U 
IN  19 
IN  16 
IN  17 
IN  19 
In  19 
IN  20 
IN  21 
IN  22 
IN  2J 
In  24 
IN  25 
IN  26 
IN  27 
IN  29 
In  29 
IN  50 
IN  31 
IN  32 
IN  33 
IN  34 
IN  35 
IN  36 
IN  37 
In  38 
In  39 
In  40 
In  at 
In  a2 
IN  43 
IN  44 
In  a 5 
In  46 
IN  47 
IN  46 
In  49 
In  50 
IN  51 
IN  52 
IN  53 
IN  54 
In  55 
IN  56 


¥ 
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F 


IF  (DEI. IT. 0.)  08«-08 

IN 

57 

9 

X8s08 

IN 

58 

x9=£l«EL»EL/12. 

IN 

59 

Xl0a,b*( ( sa- 1 ,5a»)*P2-(31M ,^>‘B)*ai)+,12S»(i,»B2-tt,*C)»D2 

IN 

60 

XI  laS2-Sl-B*ALR*,5a(B2-2.*C>*03 

In 

bl 

Xl2r02-3»D8-C»0<l 

IN 

62 

KETIINN 

IN 

63 

C 

IN 

btt 

C 

IN 

65 

10 

FOWMAT  (1X,«5MINTeGR*TJON  ATTEMPTED  OVErf  SINGULARITY,  B.CS/2E15.6 

IN 

66 

1) 

IN 

67 

END 

IN 

63 

FUNCTION  IStGNO  ( I T *GI , M ) 

IS 

1 

c 

IS 

2 

c 

I SEGNO  net uR  NS  T«t  SEGMENT  NUMBER  OF  The  mTm  SEGmEnT  HAVING  ThE 

IS 

3 

c 

T»G  M/mBeR  ITAGI,  IF  I T AG  1 2 J SEGMENT  NUMBER  H IS  RETURNED# 

is 

a 

c 

IS 

•> 

COMMON  /DATA/  N,NP,X(200),r(200),Z(200),Sl<200),Bl(200),ALP<200), 

IS 

6 

I BET(200)»  ICE'Nl  (200),  IC''N2<200),  ITaG(200),IPx.IPY,|RZ 

IS 

7 

IF  (m.GT.O)  GO  TO  1 

IS 

S 

PRINT  5 

IS 

9 

STOP 

IS 

to 

I 

I CN  TsO 

IS 

11 

IF  (ITaGI.NE.O)  go  TO  2 

IS 

12 

ISEGNOsM 

IS 

13 

RETURN 

IS 

1« 

2 

DO  3 1*1, N 

IS 

15 

IF  (ITAG(I),NE. ITAGI)  GO  TO  3 

IS 

16 

ICNTsICNTM 

IS 

17 

IF  (ICnT.EQ.m)  GO  TO  a 

IS 

18 

3 

CONTINUE 

IS 

19 

PRINT  6,  ITAGI 

IS 

20 

STOP 

IS 

21 

a 

ISEGNOsI 

IS 

22 

RETURN 

IS 

23 

C 

IS 

2« 

S 

FORMAT  (UX.91HCHECK  DATA,  PARAMETER  SpECIFYING  SEGMENT  POSITION  IN 

IS 

25 

1 A GROUP  OF  EQUAL  TAGS  mu$»  NOT  BE  ZERO) 

IS 

26 

6 

FORMAT  (///, 10X.26HNU  SEGMENT  MAS  AN  ITAG  OF  ,15) 

IS 

27 

END 

IS 

28- 

A-40 
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SUBROUTINE  I TOE  ( T RA N , N T , « , NOE , NOC , NTNX * N TNR ) 

IT 

IT 

1 

2 

ITt'F  PBOCFSSES  DATA  IN  ARRAY  TPAN  FOB  FOURIER 

TRANSFORMING  and 

IT 

3 

CALLS  SUhBUUT  I NF.  Fo»T  TO  DO  T*E  TRANSFORMING 

IT 

IT 

a 

5 

COMMON  / ARH  A X / f,It,tiI2,uM,ET(200,5,8),EP(200 

,3,8),IRET(200,B) 

IT 

6 

DIMENSION  TBaM 1 0 2 U ) , A(a09o),  S(5|l) 

IT 

7 

DIMENSION  TZERU(IOO),  AMAX(IOO).  TEXT(IOO) 

IT 

8 

EQUIVALENCE  (A, Til),  ( TZER3,  £P  ( 1 » 1 , 2 ) ) » (AMAX, 

EP (1,1.3)).  (TEXT ,£P 

IT 

9 

Kixi. «n 

IT 

10 

OATA  MM*X, ITZX/0. 100/ 

IT 

11 

PI  = S, 1B15926 

IT 

12 

TP  = 2 . *PI 

IT 

13 

lFEBBso 

IT 

la 

Ns?*‘M 

IT 

IS 

PRINT  2«,  NT.N 

IT 

16 

N2  = N*2 

IT 

17 

NTPsNTtl 

IT 

18 

00  1 I a 1 # NT 

IT 

19 

A(I)sTBAN<I) 

IT 

20 

IF  (NT.LT.N)  GO  TO  2 

IT 

21 

NTsN 

IT 

22 

GO  TO  12 

IT 

23 

IF  (NOE.NE.O)  CO  TO  12 

IT 

IT 

2a 

23 

EXAMINE  DATA  FOB  POSSIBLE  E X TN APOu A T ION 

IT 

IT 

26 

27 

NTMsNT »1 

IT 

28 

I T Z 3 0 

IT 

29 

ITMSO 

IT 

30 

DO  B IsJ.NTH 

IT 

31 

TMAsA C I — 1 ) 

IT 

32 

TBBsA(I) 

IT 

33 

TRCsA ( I ♦ l ) 

IT 

3a 

ATRAsACS(TRA) 

IT 

33 

IF  (ATNA.LT.1.E-3S)  GO  TO  3 

IT 

36 

IF  <1.EM0«TBA»TBB.GT,0.)  GO  TO  3 

IT 

37 

ITZ*ITZ*1 

IT 

38 

IF  (1TZ.GT.1TZX)  GO  TO  II 

IT 

39 

T Zt  R')(  1 T Z )=FLO*T  ( I»2)*ABS(  TR*/(TRA"TRB)  ) 

IT 

ao 

ATPAsAmS(TBA) 

IT 

a l 

ATRRs  *fi$(TR8) 

IT 

a2 

ATBCaAliS(TBC) 

IT 

US 

IF  ((ATHB.LE.ATBa).OB,(AT«O.LT,ATBC))  GO  TO  a 

IT 

uu 

TM*xisTB*»2,«T»B*TBC 

IT 

as 

IF  (ams(T“aXI),LT,I.E“3S)  GO  TO  a 

IT 

a 6 

ITPslThM 

IT 

u7 

IF  (ITM.GT.ITZX)  GO  TO  11 

IT 

U8 

|MA»is-,<j«(JBC-1*A)/TMAXl 

IT 

a 9 

TEXT  (I  T-.JsTMAXl  .FLOAT  ( 1-1  ) 

IT 

50 

AMA<(lTM)a(,S*T«A-TKB*,b»TBC)»TMAXl»TMAXlA(T«C 

•TRA)»,3*TmaxI»TRB 

IT 

51 

CONTINUE 

IT 

52 

IF  (1TM.LT, 3)  GU  TO  11 

IT 

53 

ALF  AsO, 

IT 

5a 

ICOL'nT*0 

IT 

55 

I53S«2 

IT 

36 
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AHAXL*AMAX(2) 

IT 

57 

00  7 Isi.iTH 

IT 

58 

ANAX 1 3AMAX ( I ) 

IT 

59 

IF  <»H*XI.A'UXl.GE.Ot)  GO  TO  5 

IT 

60 

icount=icount*i 

IT 

61 

AHAXSa-AMAXl/AHAXL 

IT 

62 

IF  ( IMmH  ,r>T  , 1 , ) GO  TO  5 

IT 

61 

AlFAeAuFA*ALOC(AMAX<»)/(TEXT(I-l)-TEXT(I)) 

IT 

64 

GO  TO  b 

IT 

65 

5 

ICOUNTeO 

IT 

66 

ISSS*I 

IT 

67 

ALFAsO. 

IT 

68 

6 

A«AXLSAHAXI 

IT 

69 

7 

CON  T 1 NUt 

IT 

70 

IF  UCOUNT.UT.il  GO  TO  11 

IT 

71 

ALFAsAlFa/FCOAT (ICOONT) 

IT 

72 

THAXJsTt'XTCISSS) 

IT 

71 

DO  8 1=1, ITZ 

IT 

74 

I COUNT  = 1 

IT 

75 

IF  (TZKHOUl.CT.TMAXI)  GO  TO  9 

IT 

76 

0 

CONTINUE 

IT 

77 

GO  TO  11 

IT 

78 

9 

PERIs2.*(TZER0(ITZ)-TZER0(IC0UnT))/FL0*T(ITZ-IC0UNT) 

IT 

79 

OMtGAs)P/P£Hl 

IT 

80 

TZEHX  = TZE*OUTZ> 

IT 

81 

TENOsFLO*T(NT-l  1 

IT 

82 

AMXsa(i.T)/SIn(ON£Ga»(TEND-TZEHX)) 

IT 

81 

print  21,  alFa, OMEGA 

IT 

84 

C 

IT 

85 

C 

E x T 9 A POL  ATE  -1TH  ATTENUATED  SINE  WAVE 

IT 

88 

C 

IT 

87 

00  10  I SNTP » N 

IT 

88 

TIHsf LOAT(I-l) 

IT 

89 

10 

A<  I)sanxaEXP(-auFA«(TIh-TEN0) ) »3IN (OMEGA* (TIM-TZERX) ) 

IT 

90 

NTX*N 

IT 

91 

GO  TO  19 

IT 

92 

11 

PRINT  22 

IT 

91 

12 

NTXsNT 

IT 

94 

C 

IT 

95 

C 

NO  EXTRAPOLATION 

IT 

96 

c 

IT 

97 

00  11  ISNTP#N 

IT 

98 

13 

A ( I ) S 0 . 

IT 

99 

l« 

IF  (NOC.NE.O)  GO  TO  16 

IT 

100 

c 

IT 

101 

c 

APPLY  HALF  COSINE  BELL  TAPER  TO  LAST  1/10  OF  POINTS 

IT 

102 

c 

IT 

101 

NTPsNTx/10 

IT 

104 

IF  (NTP.LT.l)  GO  TO  16 

IT 

105 

nTksnt X-NTP*  | 

IT 

106 

TRAsMI/FuOAT(nTP*1) 

IT 

107 

TRbsO, 

IT 

108 

00  15  I*NTM,NTX 

IT 

109 

TRdsTSs*  t o a 

IT 

110 

15 

A(I)aA(I)*(l,*C0S(TR8))»,5 

IT 

111 

16 

NTPs\*  1 

IT 

112 

00  17  1 1 1 # N 

IT 

111 

ITZ*NTP-I 

IT 

114 

I TMs/* I TZ 

IT 

115 

AUTmj.O, 

IT 

116 

A- 4 2 


l 7 

A(ITN-1)=A(ITI) 

IT 

117 

IF  (m.lE.hm**)  CO  to  ta 

IT 

110 

MMAXsM 

IT 

119 

)CI)iO, 

IT 

120 

IFS=-1 

IT 

121 

CO  TO  19 

IT 

122 

10 

IFSa-2 

IT 

123 

CALL  FORT  <4,H,S, IF3, IPEHR) 

IT 

12« 

IF  (IFEH9.NE.0)  PRINT  23,  IFEhR 

IT 

1 25 

IF  (IFE9H.E0.1)  STOP 

IT 

12b 

nTnHsn 

IT 

127 

IF  (N.GT.nTNx)  NTNRsNTNX 

IT 

120 

00  20  1=1 , NTNR 

IT 

129 

20 

THAN ( I } S A ( I ) 

IT 

130 

RETURN 

IT 

131 

C 

IT 

132 

C 

IT 

133 

21 

FORMAT  { 1X.04HINPUT  OATA  EXTRAPOLATED  USING  ATTENUATED  SINE  FUN 

IT 

134 

ICTIOn  EXP(-A*T)*SIN(**T)  HITH  A*DT=,E12.5, 10M  ANO  n *0 T = , E 1 2 . 5 > 

IT 

135 

22 

FORMAT  C1X,27HINPUT  OATA  NOT  EXTRAPOLATED) 

IT 

136 

21 

FORMAT  (//,  1X,46HERR0R  FLAG  RETURNED  br  SUBROUTINE  FORT,  IFeRr=,I 

IT 

137 

19,/) 

IT 

130 

24 

FORMAT  ( /,  lx,21«P0U»IE«  TRANSFORM  - - , / . 1 X , 23HNUMUER  OF  POINTS 

IT 

139 

1 INPUT*,  17,/, IX.29MNUM9ER  OF  POINTS  TH ANSFOHMEOa , 1 7 ) 

IT 

140 

ENO 

IT 

141- 
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o n o o n n n %j\  c:  w n ••  o o r* 


SUbHOUTlN  '.FACTS  (A,N,IXl,IX2,P>  IF  | 


SUbHOUT  l S ' '.FACTS  (A,  N, 1X1, 1X2, P)  IF  I 

C IF  2 

c I)/*.  -HJCM  PERFORMS  GAUSS-OOO;.  I T T u F.  MANIPULATIONS  ON  THE  T»0  BLOCKS  LF  J 

C UNFACTOHu  OS  PakTIALLV-FACTOREO  haTRIX  hhICm  IS  IN  THE  CCRfc  SIORA  LF  u 

C LF  5 

C THE  MANIPULATIONS  ARE  BASEO  ON  THE  C-AUSS-OOOl  I T TIE  ALGORITHM  PHESE  LF  6 

C PAGES  9 1 1 -9  1 6 OF  A.  RALSTON  — A FIRST  CO'JRSE  IN  NUMERICAL  ANALYST  LF  7 

C COMMENTS  RELO*  REFER  TO  COMMENTS  IN  RALSTONS  TEXT.  LF  8 

CCmmON  /maTPAR/  NBLL'kS.nPBLK.NLAST,  INT  LF  9 

COMMON  /SCmATM/  0(200)  LF  10 

0*  mens J ON  • ( n , N ) , P ( N ) LF  11 

TNTFGEF  H, R1 ,«2,P,PJ,PR  LF  12 

LF  1 J 

II  LGsO  LF  19 

LF  IS 

INITIALIZE  R 1 , P2 , J 1 , J2  LF  16 

LF  17 

LUU1.E0.1  ,ANU.  IX2.E0.2  LF  IB 

L2a(IX?-l).EQ.lXl  LF  19 

L TsIxP.EO.kBlOKS  LF  20 

IF  (LI)  GO  TO  1 LF  21 

C-0  TO  2 LF  22 

F l * 1 LF  21 

P2s2«NPBLK  LF  29 

JlPl  LF  28 

JM-l  LF  26 

G J TO  8 LF  27 

K.si»RbLK*t  L F 2 B 

v?.z..  'Knc.r  ' F 29 

,J  s ( ( * 1 * 1 J *NP8wK+  1 Lc  TO 

If  (L2 ) »0  TO  i LF  11 

GO  TO  0 LF  52 

J2=J1 +NPPLK-2  LF  31 

GO  TO  5 LF  39 

J2SJ1ANPSLK-1  LF  35 

Ic  (111  R2=NP6LK+NLAST  LF  56 

J?P2S=J2>2  LF  37 

00  16  H*n,R?  LF  38 

LF  59 

SUP  1 LF  90 

LF  «1 

00  6 Ks^i,N  LF  <12 

0(K)sA(K,R)  LF  91 

CONTINUE  -F  «« 

LF  US 

STEPS  2 AND  3 LF  96 

LF  .9? 

IF  f1.  •,0®,L2)  J 2 * J 2 M LF  as 

• • . , ' T . • p > Ml  T *1  9 i.F  UV 

l < 30 

t . A A J s J i , J i tF  Si 

1 •' 

Pvl9Pf.j)  ' 

AT J.P)rO(PJ)  L 

C J)*C(J) 

J J 1 » J M l F 1 •’ 


00  7 UJPl.N 

LF 

57 

P(I):P(l)-A(l,lXJ)*MJ,R) 

LF 

56 

7 

COM  T I Nut 

LF 

59 

6 

CONTINUE 

LF 

90 

9 

CONTINUE 

LF 

91 

C 

LF 

92 

C 

STEP  « 

IF 

8.1 

C 

LF 

69 

J*ri«J2M 

LF 

95 

I 7 (L1.0N.L2)  CO  TO  11 

LF 

66 

IP  (N.LT.^Pl)  CO  TO  1 o 

LF 

97 

00  10  l*J2u\,* 

LF 

63 

*(  j.wisom 

LF 

69 

JO 

CONTINUE 

LF 

70 

CO  W 16 

IF 

71 

!1 

0N»«  = »'JS(0(J2P1)) 

LF 

72 

PU2P1  >=J2P1 

LF 

7 3 

LF 

79 

If  (J2P2.GT.Nl  GO  TO  1J 

LF 

75 

oo  ia  i=J2Pa,N 

LF 

76 

ELN*G  = *t)  '5(0(1)) 

LF 

77 

IP  (tLN4C.LT.DMA*)  GO  TO  12 

LF 

76 

Oha«=Elm*C 

LF 

79 

P(  J2*»l  )«I 

LF 

60 

12 

CONTINUE 

LF 

81 

JJ 

CONTINUE 

LF 

62 

IF  (Onax.LT.I.E-10)  IFLCsl 

LF 

6 i 

PNsP(.)2Pl  ) 

LF 

6 9 

A(J2P1 ,N)=0(PH) 

LF 

65 

D(PP)=0( J2P1 ) 

LF 

66 

C 

LF 

S7 

c 

STEP  5 

LF 

86 

c 

LP 

89 

IP  (J2P2.CT.N)  CO  TO  15 

LF 

90 

DO  1«  I=J2P2.N 

LF 

91 

A(l.N)sO(I)/A(J2Pl,R) 

LF 

92 

M 

CONTINUE 

LF 

9 J 

is 

CONTINUE 

LF 

99 

IF  (IFLG.EO.O)  CO  TO  16 

LF 

95 

PRINT  17,  J2,DMAX 

LF 

9b 

IFLC=0 

LF 

97 

19 

CONTINUE 

LF 

98 

ketuhn 

LF 

99 

C 

LF 

100 

17 

FOHNAT  (1M  #'PIV0T(',I3,')B',E16.8) 

LF 

101 

END 

LF 

102 
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SUBROUTINE  ITSOIV  ( 4 , A R , N , I X , 8 ) 

LT 

1 

c 

LT 

2 

c 

3/R  TO  SOLVE  THF  maTRIx  EQUATION  Y ( R ) »LU  ( T ) a9  ( 9 ) wMERF  (R)  DENOTES 

LT 

3 

c 

HO-<  VECTOR  and  LO(T)  DENOTES  Tuf  LU  DE'COmoOSITIOn  QF  TmE  TRANSPOSE 

LT 

4 

c 

OF  The  ORIGINAL  COEFFICIENT  matrix,  Tne  LUtT)  DECOMPOSITION  IS 

LT 

3 

c 

STORED  ON  FILE  13  IN  BLOCKS  In  ASCtN0I~G  ORDER  and  ON  FILE  U In 

LT 

ft 

c 

BLOCKS  OF  DESCENDING  ORDt R , 

LT 

7 

c 

LT 

8 

c 

LT 

9 

COMMON  /KaTF'aR/  NBLOkS,NPULK#NLA3T«  INT 

LT 

10 

COMMON  /SCRATM/  Y (200 ) 

LT 

1 1 

DIMENSION  A(N,N),  e(N),  1 X ( N ) , A R ( 1 J 

LT 

12 

c 

LT 

13 

c 

FORWARD  SUBSTITUTION 

LT 

14 

c 

LT 

15 

11  = 1 

LT 

16 

I 2 = npdlk  »N 

LT 

17 

JeO 

LT 

18 

DO  « IXBLKl=l*NBLOKS 

LT 

19 

READ  ( 1 3 ) (AR (I), 1=11,12) 

LT 

20 

K2=NP8L* 

LT 

21 

IF  (IXBLK1 .EQ.NBLOKS)  K2sNLA3T 

LT 

22 

DO  3 «sl,K2 

LT 

23 

JMjeJ 

LT 

24 

J = JM 

LT 

25 

SUM=0. 

LT 

28 

IF  (JMl.LT.l)  GO  TO  2 

LT 

27 

DO  1 I=1,JM1 

LT 

29 

SUm*SUmaa(I,K)*Y(I) 

LT 

29 

1 

CONTINUE 

LT 

30 

2 

CONTINUE 

LT 

31 

Y(J)=(B(J)-SUM)/A(J,K) 

LT 

32 

3 

CONTINUE 

LT 

33 

4 

CONTINUE 

LT 

34 

c 

LT 

35 

c 

BACmmaRO  SUBSTITUTION 

LT 

36 

c 

LT 

37 

J»N  + 1 

LT 

38 

DO  8 I xBLK  1 a 1 * NBLOSS 

LT 

39 

READ  (14)  (AR(I),l  = n,i2) 

LT 

40 

K2=nPBLK 

LT 

41 

IF  (MHLKl.EU.l)  R2  = NLAST 

LT 

42 

DO  7 K a l , K 2 

LT 

43 

kPsA2-<a i 

LT 

U<4 

JR  1 a J 

LT 

45 

JaJ-l 

LT 

46 

SUMaO, 

LT 

47 

IF  (n.lT.JPII  CO  TO  6 

LT 

48 

DO  5 I a JP 1 , N 

LT 

49 

SUMaSUMAAd.RPjASd) 

LT 

SO 

s 

cuntinuE 

LT 

51 

e 

CONTINUE 

LT 

52 

B(J)aY(J) • SUM 

LT 

53 

7 

CONTINUE 

LT 

S 4 

8 

CONTINUE 

LT 

55 

C 

LT 

So 

A- 4 6 

J 

c 

UNSCRAMBLE  SOLUTION 

LT 

57 

c 

LT 

se 

DO  9 Isl.N 

LT 

59 

IXIsIX(I) 

LT 

60 

V(lXI)sBU) 

LT 

61 

9 

CONTINUE 

LT 

62 

DO  10  I a 1 » N 

LT 

63 

10 

b(I)aV(I) 

LT 

64 

R£*lNO  15 

LT 

65 

RO-INO  14 

LT 

6b 

RETURN 

LT 

67 

END 

LT 

ha- 
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SUBROUTINE  LUNSCR  ( * , AR , N , I X , IP ) 

LU 

1 

c 

LU 

2 

c 

3/R  "MICH  UNSCRAMBLES,  SCRAMBLED  FACTORED  MATRIX 

lu 

3 

c 

LU 

0 

COMMON  /mxTPaR/  NBLOKS.NPMLRjNLAST, I NT 

LU 

s 

DIMENSION  A ( N , N ) , I P C N ) , IX(N),  AR(1) 

LU 

6 

Il«t 

LU 

7 

I2=nPSL**N 

LU 

8 

NM| sN- | 

LU 

4 

00  0 l x8LK 1 s 1 # N8L0K S 

LU 

10 

READ  Ml)  (AR(I),taIl,l2) 

LU 

1 1 

K|s ( I XtlLK 1 " 1 ) *NPaL**2 

LU 

12 

IF  (NM1.LT.K1)  CO  TO  5 

LU 

13 

J2  = 0 

LU 

10 

00  2 KsK1,NM1 

LU 

1 5 

] F ( J2.LT. nPBLK)  J2-J2* 1 

LU 

16 

I»K3lP(K) 

LU 

17 

00  l Js|,J2 

LU 

16 

TEMPSA(K.J) 

LU 

19 

A(K,J)=A(IPK,J) 

LU 

20 

A (IPK, J)aT£MP 

LU 

21 

1 

CONTINUE 

LU 

22 

2 

CONTINUE 

LU 

23 

3 

CONTINUE 

LU 

20 

WRITE  (13)  (AR(I),I*I1,I2) 

LU 

25 

• 

CONTINUE 

LU 

26 

on  s l xbl* 1 *1 ,nrlOk3 

LU 

27 

BACKSPACE  13 

LU 

28 

IF  (lXBUKl.NE.l)  BACKSPACE  13 

LU 

29 

READ  (13)  (AB(I),JsIl,I2) 

LU 

30 

WHITE  (10)  (ARID, 1311, 12) 

LU 

31 

5 

CONTINUE 

LU 

32 

DO  6 I * 1 » N 

LU 

33 

I A ( I ) * I 

LU 

30 

6 

CONTINUE 

LU 

35 

00  7 I « 1 , N 

LU 

36 

I'lsiPd) 

LU 

37 

IXTstX(I) 

LU 

38 

I*CI)»IX(IPI) 

LU 

39 

I*(JPI)*I*T 

LU 

OO 

7 

CON . j nuE 

LU 

01 

RE  » I NO  11 

LU 

02 

RE » I NO  13 

LU 

03 

R f *•  l N 0 10 

LU 

06 

RETURN 

Lu 

06 

£N0 

LU 

08- 

A- 4 8 


c 

c 

c 

c 

c 

c 


1 

2 


SUBROUTINE  ( ROX , RO Y , ROZ , X$, YS , Z S , I T S, NRP T , I TC I ) 

SUBROUTINE  MOVE  MOVES  THE  STRUCTURE  WITH  RESPECT  TO  ITS 
COORDINATE  SYSTEM  OR  REPRODUCES  STRUCTURE  IN  N£*  POSITIONS, 
STRUCTURE  IS  ROTATED  ABOUT  X,Y,Z  AXES  BY  ROX, ROY, ROZ 
RESPECTIVELY,  Then  ShIETEO  BY  XS,vs,ZS 

COMMON  /DATA/  N,NP,X(200).Y(20C),Z(200), 31(200), 01(200),  ALP(200), 
1 BE  T(  200),  ICON  1 (200),  I CON2  ( 200  ) , I T AG  ( 20  0 ) , IPX,  IP  Y,  If*  Z 
DIMENSION  >2(1),  Y2(l),  22(1)  ' 

EQUIVALENCE  ( X2 ( 1 ) , SI  ( l ) ) , ( Y2 ( l ) , ALP ( 1 ) ) , < Z2 C 1 ) . BET < 1 ) ) 

SPS^SINIROX) 

CPSsCOS(ROX) 

STHsSIN(ROY) 

CThsCOS(ROY) 

SPh*SIn(H,j2) 

CPHsCOS(ROZ) 

X*rCPH«CTH 

XYsCPh*3)m»sp9-sph«CP3 

XZ=CPm*STh»CPS»SPh»SP3 

YX=SPh»CTH 

YY=SPH«STh*SPS*CPH«CPS 

YZ=SPh«STh*CPS-CPh*SPS 

ZX  = -STm 

ZY=CTH*SPS 

ZZbCTH»CPS 

1 1 s I SEGNO ( I TS,  1 ) 

IXsIl 

NRPsNRPT 

IF  (NRPT.EU.O)  NRPsl 
Ken 

IF  (NRPT.EO.O)  K*I1-I 
DO  2 IR=1,NHP 
00  1 I a 1 1 , N 
K sK  ♦ 1 
XI»X(I) 

Y 1 3 Y ( I ) 

ZI=Z1I) 

X(K)3XlNVY+YI*KY*ZI*XZ*XS 
Y(k)3*i«yx+yi*yy»ZI*y2»YS 
2(k)=xi*2x+Yi«ZYyZI*ZZa2S 
X 1 3X2 ( I ) 

YI=Y2( I ) 

ZI=Z2(I) 

X2(K)s*I*XXYYI*XYyZI*XZYXS 

Y2(W)sxI»yx*YI*YYaZI*YZ*YS 

Z2(w)»aI*ZX+tI*ZY»ZI*ZZyZS 

BI(K)sbki) 

ITaG(k)3I1aG(I)aITGI 

CONTINUE 

I1=N*1 

NSK 

CONTINUE 

return 

END 
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HO 
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MO 
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MO 
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MO 
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MO 
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MQ 
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MO 
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SUBROUTINE  RFLD  ( THE T , PH  1 , t TM , EPh , T $ TR T , KLM ) 

RF 

1 

RF 

2 

RFLD  COMPUTES  THE  HAOlATED  FIELD  IN  THE  DIRECTION  THETA, 

PhI  for 

RF 

3 

ALL  TIME  STtPS  POSSIBLE  nITm  Th£  AVAILABLE  CURRENTS 

RF 

a 

RF 

5 

CO«»ON  /DATA/  n,np,X(200),Y(200). 7(200), SI(200),3I(200),AL PC  200), 

RF 

6 

I BET (200) , ICON! (200) , ICON2 (200) , ITAGC200) , IPX, IPV, IP2 

RF 

7 

common  /sroMP/  ox (200 ) , sy (200 ) , S2 (20o ) 

RE 

a 

COMMON  /InTEMP/  AT( 3.200) ,BT (3,200) «CTC3. 200) ,ES(3) ,FS(3) 

t US  ( 1 ) t 

RF 

9 

1 E ( 3 ) » M ( 3 ) 

RF 

10 

COMMON  /CONST/  COT, VEL.DT.NTSTEP 

RF 

11 

COMMON  /array/  011,012, Oil, ET(200,5»«),EP  (2<>0, *»,«), I RE  T(200,«) 

RF 

12 

COMMON  /ARRAY/  C 0 ( 60  0 0 ) 

RF 

13 

common  /IOFLC/  NCFmx, NCOMX, IOC , IOCU , NT M A X , NBOU T , JP1 , JP2 

RF 

la 

DIMENSION  A ( 3 ) , 6(3),  C ( 3 ) , P(S),  t TH ( l ) , EPM(l),  ALFAC3, 

3) 

RF 

15 

N2  = N*2 

RF 

16 

STHsSIM  ThET  ) 

RF 

17 

CTHsCOS( TMtT  ) 

RF 

18 

SP=3IN(PmI ) 

RF 

19 

CP=CnS(PMI) 

HF 

20 

ERx=3Th*CP 

RF 

21 

ERYsSTH»SP 

RF 

22 

ERZ=C1m 

RK 

23 

THXsCTmacp 

RF 

2a 

THY=CTH*3P 

RF 

25 

TH/s-STH 

RF 

26 

PMXS-S2 

RK 

27 

PHYsCP 

RF 

28 

KMNS9999 

RF 

29 

IPLANSO 

RF 

30 

LPXsl 

RF 

31 

LPYsl 

RF 

32 

LP2  = 1 

RF 

33 

IF  (IPX.NE.O)  LPX  = 2 

RF 

3a 

IF  (IPy.NE.O)  LPYS2 

RF 

35 

IF  (IPZ.NE.O)  LPZ»2 

RF 

36 

RF 

37 

bEG In  LOOP  OVER  SEGMENTS,  INCLUDING  IMAGES,  TO  COMPUTE 

RF 

38 

COEFFICIENTS  FOR  RADIATED  FIELD 

RF 

39 

HF 

ao 

DO  17  JXsj.LP* 

RF 

«1 

RF  xsELUAT ( 3“ J * • 2 ) 

RF 

U2 

SFx=mF  X 

RF 

a 3 

IF  (IPx.LT.O)  SFXH, 

HF 

ua 

DC  17  JYsj.LPT 

Rr 

«S 

RFYsALOAT ( 3" J V • 2 ) 

RF 

Ub 

SFYiRFy 

RF 

a 7 

IF  (IPV.LT.O)  SFYsl, 

RF 

a8 

CO  17  JZ=1,LPZ 

Rf 

a 9 

RFZsA  luat ( J-JZ»2 ) 

RF 

50 

SFZsRFZ 

RF 

51 

IF  (1PZ.LT.0)  SFZ«1. 

RF 

52 

RFL=SFX»SFY*SFZ 

RF 

S3 

rfls-i  ,e-t««fl 

»F 

5a 

IPLANajPtANAl 

RF 

55 

00  1 J>t,N 

RF 

56 
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AL»-(ERX.X(J)»RFXtERV»Y(J)*RFY4EHZ»Z(J)»RFZ) 

RF 

57 

TRET  = A|./YEL 

RF 

53 

KsTBET/OTtSlGNC  .5,  TRET) 

RF 

59 

I»ET (J, IPLANjsK 

RF 

60 

IF  (K.LT.KHN)  KMNaK 

RF 

61 

00  1 Ms  1,5 

RF 

62 

ET(J,m, IPl*n)«0, 

RF 

63 

EP(J,m, lPLA*)sO, 

RF 

64 

DO  16  J* 1 , N 

RF 

65 

XjsX(J)»RFX 

RF 

6b 

YJ=Y(J)»RFY 

RF 

67 

ZJ»2(J)*RFZ 

nF 

fc* 

SXJsSX ( J ) • RF X 

RF 

69 

SYJ:SY(J).RFY 

KF 

70 

SZJ=SZ(J)«HFZ 

RF 

71 

*t=-(f»**xJ*E«Y«Yj+ERZ*ZJ) 

RF 

72 

BA=»(tRX«SXJ«ERY»SYjAERZ*SZJ) 

RF 

73 

KsIRET  (J, IPLAN) 

RF 

74 

T AUsDT *FlO*T (K) 

RF 

75 

00  2 L= 1 , J 

R? 

76 

A(L)=*T(L,J) 

RF 

77 

B(L)3BT  CL, J) 

RF 

78 

C ( L ) =C  T ( L » J ) 

RF 

79 

P(L)=2.«ES(L)*TAU*FS(L) 

HF 

eo 

EL*SI (J) 

RF 

81 

ELC  = EL*EL»EL/'12. 

RF 

82 

SDUTTsSXJ»Tmx4SYJ*TMY4SZJ*THZ 

RF 

83 

$Of?TRsSXJ«PMXfSYJ*PHY 

RF 

84 

DO  3 M a 1 , 3 

RF 

85 

T 1 SM (M ) » AU*P ( H ) 

RF 

86 

T2sh(M)«ba 

RF 

87 

00  5 Ls 1 , 3 

RF 

88 

ALFA(L,M)s( (A(U*Tl*d(L)»T2)»ELC»C(U)*Tl*El)»REU 

RF 

89 

JC 1 s ICON  J (J) 

RF 

90 

IF  (JCl.ME.O)  GO  TO  4 

RF 

91 

KKMsO 

RF 

92 

GO  TO  8 

RF 

93 

IF  (JC1.LT. 19000)  GO  TO  S 

RF 

94 

SIG1=FlO*T(JC1-20000) 

RF 

95 

JCl=J 

RF 

96 

GO  TO  7 

RF 

97 

IF  (Icn*2(JCl),NE.J)  CO  TO  6 

RF 

93 

SIGl»l, 

RF 

99 

GO  TO  7 

RF 

100 

IF  ( ICON) ( JC | ) ,NE , J ) GO  TO  23 

RF 

101 

5 1 G 1 = • 1 , 

RF 

102 

FMalRETUCl  r I RL  AR ) 

RF 

103 

RF 

104 

IF  (Ilo5(»‘!,),l.£il)  GO  TO  3 

R r 

105 

PRINT  24,  J.JCl.KKM 

RF 

106 

STOP 

RF 

107 

JC2s 1 C0N2 ( J ) 

RF 

103 

RF 

109 

KKPsO 

RF 

110 

GO  TO  13 

RF 

1 1 1 

IF  (JC2.LT. 19000)  GO  TO  10 

RF 

112 

Sir.2sFi.OAT(  JC2-20000) 

RF 

113 

JC2«J 

RF 

114 

GO  TO  12 

RF 

115 

IF  ( ICONi ( JC2) tNE , J)  GO  TO  11 

RF 

116 
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3 I G2a  l , 

RF 

117 

GO  TO  12 

RF 

lid 

It 

IF  (1C0N2(JC2).NE,J)  GO  TO  23 

RF 

119 

3IG2s-l  , 

RF 

120 

12 

APelRET(JC2. IPLaN) 

«F 

121 

A AP* AP«A 

RF 

122 

IF  U*RS(KKP),LE,1)  GO  70  13 

RF 

123 

PRINT  20,  J.JC2.KKP 

RF 

12a 

STOP 

FiF 

125 

13 

oo  15  m»i , i 

PF 

126 

AXX  = M*  | 

RF 

127 

IF  (.IC1  ,FOtO>  GO  TO  10 

RF 

123 

KMX  * AKM6 A X X 

RF 

129 

ETUCl ,Amx« IPL*N)seT(JC1,KpX,IPL*N)**UF*(i , h ) • SOOT  T *3 1 G 1 

WF 

130 

tP(JCl»K^X.lPL*N)rEP(Jci,KMX, IPLAN)*ALFA(1,M)*5D0TP*5IG1 

RF 

131 

l« 

tT(.I,KXX,IMLAN)si.T(J,KXX,  IPLAH)*aLFA(2,M)»SOUTT 

RF 

132 

EPlJ.KXX,IPL*N)=tP(J»KXX, IPLAN) ♦ALF* (2»H) *SOOTP 

RF 

133 

IF  (JC2.EC.0)  GO  TO  15 

RF 

1 3a 

KPXsKAP+AXX 

RF 

135 

ETUC2.KPX.  IPLAN)=ET(JC2,KPX, I PL  4*0  ♦ AIK  A ( 3,M)*SDQTT*3IG2 

RF 

136 

EP{  JC2.KPX,  IPLAN)stPlJC2»KPX,  IPLASO  + ALFAt  J,.-1)»S00TP*SIG2 

RF 

137 

15 

CONTI  NuE 

RF 

i 35 

16 

CONTINUE 

RF 

139 

17 

CONTINUE 

RF 

lao 

C 

RF 

1 a l 

c 

COMPUTE  RADIATED  FIELD  AT  ALL  TIME  STEPS  POSSIBLE  WITH  SUMMATION 

»F 

1 a 2 

c 

PHuCEEDlNG  IN  ORDER  OF  RtTAROEO  TIMES 

RF 

10  3 

c 

RF 

laa 

KMNsKMN-1 

RF 

105 

AKMs  1 

RF 

106 

KKPsNTSTtP 

»F 

107 

KMsKKM»KMN 

RF 

las 

APsAAP^KMN 

RF 

109 

X XX*WM-l 

RF 

150 

klmskp-xxx 

RF 

151 

TSTRT=DT*FL0AT(KXX) 

RF 

152 

DO  IS  K«1,*LM 

RF 

153 

ETM(K)«0, 

RF 

150 

l# 

EPrtlXJsO. 

RF 

155 

A TSs  1 

RF 

156 

N0CQ=9999999 

RF 

157 

IF  UUCG.EO.O)  GO  TO  19 

RF 

158 

REMIND  11 

RF 

159 

RE  * 0 (11)  nOcQ  » ( C Q ( J ) » J = 1 » NOCO ) 

RF 

160 

19 

DO  22  kCskkm.AAP 

RF 

161 

20 

ACt>S(XC*l)*N^»XTS 

RF 

162 

IF  ( ( xcB ♦ | ) . L T , NOCO ) GO  TO  21 

RF 

163 

ATSsaTS^.OCQ 

RF 

160 

READ  111'  f.OCD.  (CCt.M.Jsl.MOCO) 

sc 

165 

CO  TO  20 

RF 

166 

21 

kku3XC~axx») 

RF 

167 

DO  22  1PLNE=1 , 1PLAN 

RF 

163 

DO  22  J»!,n 

RF 

169 

ACC  *aCB*2* J 

RF 

170 

kskxm*I»ET(J. IRLNE) 

RF 

171 

DO  22  m»i,5 

RF 

172 

RPXsA-M 

HF 

173 

IF  (APX.LT. 1)  GO  TO  22 

RF 

170 

IF  (APX.GT.alm)  GO  TO  22 

RF 

175 

ETh(APx)sETh(ap*)*eT(J,m, IPLNE)«CO(ACC) 

RF 

176 
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kPM(KPX)BEHH(KPX)*EP(J#M»  IPlNE)<*CO<KCC)  «F  177 

22  CONTINUE  RF  178 

RETURN  OF  179 

23  PRINT  25  »F  180 

STOP  PP  181 

C PP  182 

24  FORMAT  ( /, 1X.26MHETAR0E0  TIMES  TO  SEGMENTS, IS, UH  AN0,IS,  »P  183 

I 10H  OlPPtR  br,I5,llM  TIME  STEMS)  RP  184 

25  FORMAT  </, JX,62H*»»ERROR  - ROUTINE  PFIP--  INCONSISTENT  SEGMENT  CON  HF  165 

1 NEC T ION  DATA)  RF  186 

ENO  PP  187» 
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SUBROUTINE  RFPAT  ( J TRAN, IPCM, NTH, NPH, THET . PHI 

1 

DTH, 

DPH) 

RP 

1 

RP 

2 

RFPAT  CALLS  RELO  TO  COMPUTE  THE  RADIATED 

FIELD 

AT 

A 

SET  OF  ANCLES. 

RP 

3 

THETA  A‘.'c  phi,  AND  PRINTS  The  RESULTS. 

THE  COMPUTED 

field  Can 

RP 

a 

ALSO  BE  FOURIER  TRANSFORMED  TO  The  FREQUENCY 

DOMAIN 

RP 

5 

RP 

6 

COH“ON  /ESORC/  ESOPC(1029),IFST,ntRaN,mtRan, 

NTHX, 

DFRQ.ENIN.ENRD,  I 

RP 

7 

1FEN 

HP 

a 

COMMON  /CHAT/  CURF<lbO0),IPC200),IX(200> 

RP 

9 

COMMON  /CONST/  CDT,VEL,DT,NTSTEP 

RP 

10 

COMMON  /SCRATm/  F*H(S12) 

RP 

1 1 

DIMENSION  ETMt 1029) , EPH ( 1 029 ) , E TC  C S 1 2 ) 

, EPC (512) 

# 

ESC (Sl2) . I POL 

RP 

12 

1(1) 

HP 

13 

complex  tpj,carg,faz,ethc,ephc,etc,epc 

, ESC 

RP 

10 

EQUIVALENCE  (ETh.CURF),  (ETC.CuRF),  (EPC 

,EPM) 

« 

(ESC, 

ESORC) 

HP 

IS 

DATA  TA,TSMIN,TPJ/,0l7O53292.5.,C0,,6, 

28 J 1 853 ) / 

RP 

16 

DATA  IPQL/6HL  InEaR, 5HRICHT,<iMLEFT/f  16LK/IM 

/ 

RP 

17 

TSTRTsO. 

RP 

18 

NFLOaO 

RP 

19 

ThsThET-OTH 

RP 

20 

DO  lb  IT=1,NTH 

RP 

21 

TH*Th+dTH 

RP 

22 

THAs  TH*  TA 

RP 

23 

PHsPHl -DPH 

RP 

29 

DO  IS  I Phs 1 , NPH 

RP 

25 

PmsPh+DPh 

RP 

26 

pma=ph*Ta 

RP 

27 

PRINT  is,  th.ph 

«p 

28 

call  rflo  CTHA,PMA,ETH,EPH,TSTRT,NFL0) 

RP 

29 

NPR  TsNFLD/2 

RP 

30 

NADDsNPRT 

RP 

31 

IF  (2*NPRT,NE,NFL0)  NAOD»NPRT»l 

HP 

32 

TIm.tSTRT-DT 

RP 

33 

TADDsDT‘FLOAT(NADO) 

RP 

39 

DO  1 11*1, NPRT 

RP 

35 

I2sII*NAOO 

RP 

36 

TIMaTIHfOT 

HP 

37 

TIh2sT1M*taDD 

RP 

38 

POL1s*TCn2(EPH(I l),FTH(It))/TA 

RP 

39 

POL2=ATCn2(£PH( 12) ,ETh(I2))/TA 

RP 

90 

PRINT  19,  II.TIM.EThUD.EPHUI  ).,P0L1,I2,TIM2 

,ETH(I2),EPM(I2),POL 

RP 

91 

12 

RP 

92 

IF  (NAOO.tO.NPKT)  GO  TO  2 

RP 

93 

P0l1s*TGn21£Ph(naDD).ETH(MADD) )/TA 

RP 

<IU 

TImsTSTRT»OT«FLOaT(NP«T) 

RP 

95 

PRINT  t<5,  NAOD,  TIf; .EThINADD)  ,£PH(NADO) , 

POL  1 

RP 

NO 

KP 

07 

COMPUTE  RADIATED  ENERGY 

RP 

98 

RP 

99 

F«0, 

HP 

50 

DO  3 I a l , nF LD 

HP 

51 

Fsf  ♦ETM(U»ETH(I)*£PH(I)aEPMCI) 

RP 

52 

FaF.DT/176,7 J 

HP 

53 

PRINT  f\,  F,Th,Pm 

OP 

59 

IF  (IFtN.EQ.O)  CO  TO  0 

RP 

55 

POL  1 *F /Enin 

RP 

5b 
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POL2*F/EnRD 

RP 

57 

IF  (POLl ,LT, l.fc-20)  GO  TO  fl 

RP 

58 

POL  1 s 1 0 , » *L0G1 0 (POL  1 ) 

RP 

59 

POL2=IO.»*LOGJO(POL2) 

RP 

60 

PRINT  ?0,  PULl,POL2 

RP 

61 

4 

IF  (IPCH.EQ.O)  GO  TO  5 

RP 

62 

c 

RP 

65 

c 

PUNCH  field  VALUES 

RP 

64 

c 

RP 

65 

PUNCH  26,  TSTRT,0T 

RP 

66 

PUNCH  22,  NFLD, TH,PM 

HP 

67 

PUNCH  fi,  (ETh( I ), Ib| ,NFLO) 

HR 

66 

PUNCH  24,  NFLO,  Tm,Pm 

RP 

69 

PUNCH  2 5 , (EPh(I),I  = 1,NFLO> 

RP 

70 

5 

IF  (ITWAN.tO.O)  GO  TO  15 

RP 

71 

c 

RP 

72 

c 

FOURIER  TRANSFORM  TO  OBTAIN  FREQUENCY  RESPONSE 

RP 

75 

c 

RP 

74 

IF  (IFsT.EO.l)  GO  TO  7 

RP 

75 

MTRANsALOGt FLOAT (NTH AN) )/, 6931471861 .5 

HP 

76 

CALL  ITOF  (ESOHC,NTRan,MTRaN, 1 , I,nTnx,NTNR) 

RP 

77 

IFSTsJ 

RP 

78 

OFRO=l ,/(DT«FLOAT(2**HTRAN)) 

RP 

79 

NTNRsNTNR/2 

RP 

eo 

DO  6 l=l,NTNK 

RP 

ei 

6 

PKRdlsO. 

RP 

82 

7 

CALL  ITOF  (ETH,NTRaN,mtR*N,0,0,NTNX,nTNR) 

RP 

85 

CALL  ITOF  (EPH,NTRAN,hTRAN,0,0,NTNX,NTNR) 

RP 

84 

NTNP=NTNO/2 

RP 

85 

FHAXsl ,/( TSHIN«0T) 

RP 

66 

CARG=-TPJ»TSTRT 

RP 

87 

PRINT  17 

RP 

80 

IF  (ITRAN.GT.l)  PUNCH  25,  TH, PH 

RP 

09 

Fs-OFRQ 

RP 

90 

DO  14  I s l , nTNR 

RP 

91 

F *F  tDFRO 

RP 

92 

IF  (F.GT.FHAX)  GO  TO  15 

RP 

95 

FOsF*l.E-fc 

RP 

94 

FAZ=CExP(CARG«F)/ESC(I) 

RP 

95 

ETHC=ETC(1)»FAZ 

RP 

96 

EPhC=EhC(I)*F*Z 

RP 

97 

EThM2=hEal (ETHC*CONJG(ETHC) ) 

RP 

98 

E'»hH2sREAL(EPHC«CONJG(EPHC)  ) 

RP 

99 

ETHMsSGMT(tTHH2) 

RP 

100 

EPhmsS'JHT  (EPHM2) 

RP 

101 

EThasatGN2(AIMac(ETHC) ,REAL(ETmC) )/TA 

RP 

102 

f PmAs*TGn2( aJhag(EPHC ) » HE  *L (EphC ) ) /T  A 

RP 

105 

IF  (ETHM2,GT.l,t-20,OR.EPHH2.GTf l.t-20)  GO  TO  S 

RP 

104 

TILTA*0, 

UP 

105 

EHAJR2sO, 

RP 

106 

tHINH2sO, 

HP 

107 

AX«AT=0, 

RP 

106 

ISENS=I8LK 

RP 

109 

GO  TO  15 

RP 

no 

a 

OPAZsEPHA-tTHA 

RP 

til 

IF  (EPha.LT. 0. J GO  TO  9 

RP 

112 

1)F  AZ2=Dc  A Z"i65  . 

RP 

115 

GO  TO  10 

RP 

114 

9 

DFAZ2*r)FAZAj60, 

HP 

115 

JO 

IF  (aHs(DFAZ),GT.A0S(DFAZ2))  DFAZ«DFaZ2 

RP 

116 
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. . _ V , ’ 

cofaz3cos<ofaz»ta) 

RP 

117 

fk 

T S 1 «E  T P«M2 

RP 

1 IB 

TS2  = 2.*EPnm»ETH,i.CDFAZ 

RP 

119 

TlLT4a,5*ATGN2(TS2»TSl) 

RP 

120 

8TILT*sS1m(1 ilia) 

RP 

121 

TSlsTSl •STll!A*STIlTA 

RP 

122 

TS2»TS2*STlLTA*COS(TIiT*) 

RP 

123 

£MAjR2a-T81 *732*6 1MM2 

RP 

120 

EMlMK2=lsi-T32*EPHH2 

RP 

125 

IF  (EM1MJ2.UI.0,)  EHINR2«0, 

RP 

126 

AXHAT3S0rlT(EHlM?2/EMAjR2) 

RP 

127 

TlLTArTlLTA/TA 

RP 

120 

IF  (AxRAT.GT.l.E-5)  CO  TO  11 

RP 

129 

I SE  NS  s I POL ( 1 ) 

RP 

130 

GO  TU  13 

RP 

131 

u 

IF  (OFAZ.GT.O.)  GO  TO  12 

RP 

132 

I SE  N$s 1 POL ( 2 ) 

RP 

133 

GO  TO  13 

HP 

13« 

12 

I SENS= I POL ( 3 ) 

RP 

135 

li 

KLAMsF 

RP 

136 

IF  (k-UAM.lT.  1 ,fc-25)  «L*m=1.E-25 

RP 

131 

KLA*is2.9<)Bt*e/i*LAH 

RP 

138 

PUR=P»fi ( I ) 

RP 

139 

IF  (POM.LT.l.E-25)  POR=.0013273**LAM*WLAM 

RP 

loo 

GTOTe.0  16t>8*(ETHM2*EPHM2)/POR 

RP 

101 

IF  (GTOT.LT, 1,6-20)  GTOT=t,E-20 

RP 

102 

GTOT2lO,*ALOG10(GTOT) 

RP 

103 

E Tmme£  thm/rL  AM 

RP 

100 

EPmmsEPhm/kLAM 

RP 

10S 

PRINT  16,  t,FQ,ETHM,ETMA,EPHMfFPHA,AXRAT,TlLTA,ISENS,GTOT 

RP 

106 

IF  (ITHAN.GT.l)  PUNCH  16,  I.FO,ETHH,£TMA,EPHM,EPHA 

RP 

107 

CONTINUE 

QW 

108 

»5> 

continue 

ftp 

109 

RETURN 

RP 

150 

C 

HP 

151 

c 

RP 

152 

16 

FORMAT  (lX,I5,El2,a,El3,a,F9.j,Ei5.«,F9.3,Fll,5,F9,2,2X#Afc,«lR,3) 

RP 

153 

17 

FORMAT  ( ///, 33X , U3H-  - - FREQUENCY  DOMAIN  RAOIATEU  FIELD  - - 

RP 

150 

1 //,2X,«H3TEP,2X,9MFRtOUtNCY,3X,20H-  - - E(ThETA)  - - -,ox, 

RP 

155 

2 IBM-  - - E (Pm I ) - - -,6X,20m-  - - POLARIZATION  - - -,SX, 

RP 

156 

J IOMPOhER  GAIN,/, JX,3MnO,,UX,5h(MM2),6X,9MMAGN1TuOE,«X# 

RP 

157 

a 5hphase,6X,9hmaonITUDE,«x,5hPhaSE,7x,5haxIAl,5X,uhTILT,3X, 

RP 

158 

5 5HSFNSE,9x,«H(n6),/,22X,7MVOLTS/N,«X,7HDtGHEES,6X,7MVULTS/R, 

RP 

159 

6 <*X  , 7M0EGRLES,  6X,5m«AT  IO,5X,  UHOFG,  ) 

RP 

160 

1« 

FORMAT  ( //,O2X,30h.  - - . RADIATED  FIELO  - - - - . // , 0 7 X , 6-THE T A 

RP 

161 

1».  FB,3'./,a7x,6MPhI  «,F6.  J,//,2(2X,(iHSTEP,6X,«HTIME,0x,20MELfcCTRI 

RP 

162 

2C  FIELO  (v/M),«x,UMPOL,,8X),/,2(3x,iHNU.,5x,bH(S6C,),ex,5MTHETA,9X 

RP 

163 

3,  3hpmI,5X,6M(0EG.),7X)) 

RP 

160 

19 

FOR-AT  (2(lX,I5,Eli.3,tl3,«,El2.U,F9.3,7X)) 

HP 

165 

20 

FORMAT  ( 1 x , 2omt I RE  nilMAIN  PruFP  GAIN  z,F7,2,3H  OB. /, IX .2BRT1ME 

RP 

166 

1 UC)M  * I N DIRECTIVE  CAIN  e,F7,2,3h  DR,//) 

RP 

167 

21 

FORMAT  ( //, 1X.23HTOTAL  ENERGY  RA01ATED  s,Ell.U,17H  JOULES/STER* 

RP 

168 

lDI*N,3X,7M(ThETAs,F8,3,l2M  Ot G , , PHIa,FB.3,6H  DEG,)) 

RP 

169 

22 

format  ( I5,2x,?6MR»0lATtD  FIELDS  FUR  THE T A a , F s , 3 , l 1 H OE  G , , PMl  = ,F 

RP 

170 

16.3. IBM  DEG, , THETA  POL. ) 

HP 

171 

23 

FOR-AT  (6E12.5) 

RP 

172 

2# 

FORMAT  ( !5.2X,2bMRADlATED  FIELDS  FOR  T HE T As , F 0 , 3 . 1 1 h DEG.,  PHla.F 

RP 

173 

16.3, 1 6M  OF G , , P«I  POL.) 

RP 

170 

2S 

FOR-AT  (23MRA01ATEC  FIELD.  THE T A a , H 0 . 5 , 6«  Pnla.FtO.S) 

RP 

175 

p 

26 

FOR-AT  (lUHSTARTING  T I mE s , E 1 2 . S , 3 X , 1 5-T I -E  1 NCRE m£ NT  a , E 1 2 . 5 ) 

RP 

176 

£NO 

R? 

177- 
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uni— 

SUBROUTINE  SECOND  (T)  SE  1 

C SE  2 

C RETURNS  CURRENT  RUNNING  TIME  IN  SECONDS  SE  I 

C SE  « 

CALL  TUSEOM  (I)  SE  5 

TsI/1000,  SE  A 

RETURN  SE  7 

END  SE  6* 


SUBROUTINE  SOLVE  (N, A,P.8,nDIM) 

SO 

l 

c 

SO 

2 

c 

SUBROUTINE  TO  SOLVE  THE  MATRIX  EQUATION  LU*X=B  WriERE 

L IS  A UNIT 

so 

J 

c 

LO~£R  TRIANGULAR  KaIRI*  AND  U IS  AN  uPPE* 

TRIANGULAR 

matrix  both 

so 

4 

c 

Of  WHICH  ARE  STORED  IN  A.  THE  HH$  VECTOR 

0 IS  INPUT 

AND  THE 

so 

5 

c 

SOLUTION  IS  RETURNED  THROUGH  VECTOR  0. 

(MATRIX  TRANSPOSED) 

so 

(l 

c 

so 

7 

DIMENSION  A(NDIH,NDIM),  P(NDIM),  b(NDIM) 

so 

9 

COMMON  /SCRATM/  Y ( 200 ) 

SO 

9 

INTEGER  P,PI 

so 

10 

c 

so 

11 

c 

EORRARO  SUBSTITUTION 

so 

12 

c 

so 

13 

00  3 I a 1 # N 

so 

tu 

PlaPCI) 

so 

IS 

V(I)80(PI) 

so 

16 

BtPDsBd) 

so 

17 

I P I a I ♦ I 

so 

la 

IF  (IPl ,GT.N)  GO  TO  2 

so 

19 

DO  1 Js  I p 1 » N 

so 

20 

B(J)80(J)-A(I,J)AV(I) 

so 

21 

1 

CONTINUE 

so 

22 

2 

CONTINUE 

so 

23 

3 

CONTINUE 

so 

24 

C 

so 

25 

c 

BACKWARD  SUBSTITUTION 

so 

?6 

c 

so 

27 

DO  6 K a l , N 

so 

28 

I aN-K ♦ I 

so 

29 

SUMaO, 

so 

30 

IP  1 * I ♦ I 

so 

31 

IF  (IPI,GT.N)  GO  TO  5 

so 

32 

00  « JaIPl,N 

so 

31 

SUNaSUMA A ( J» I ) *0 ( J J 

so 

34 

a 

CONTINUE 

so 

35 

5 

CONTINUE 

so 

36 

BU)aCV(I)-SUM)/A(I,I) 

so 

37 

6 

CONTINUE 

so 

1« 

RETURN 

so 

39 

ENO 

so 

ao» 
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SUBROUTINE  TSOL  CCUPF,  IX,  IP,N) 

T 3 

I 

c 

TS 

2 

c 

TSOL  COMPUTE  3 STRUCTURE  CURRENTS 

TS 

3 

c 

TS 

a 

COMMON  /CONST/  CDT»VEL#DT.NTSTEP 

TS 

5 

Common  /tINC/  EINC(200),E3RC 

TS 

6 

COMMON  /APHAy/  C0T6U00) 

TS 

7 

COMMON  /arm kt/  OH,QI2,Ql3,EC(?00#5)ftQ(200,5),IRET(200),8F»(l5a0 

TS 

9 

10) 

TS 

9 

COMMON  /IOUG/  nCFMX/NCOMX# IOC» IOCU,NTmax,NBOuT,JP1,JP2 

TS 

10 

COMMON  /maTPaR/  NBLO<5>nP8L<*nLAST» INT 

TS 

1 1 

COMMON  /E  SOMC / ESO«C ( 1 020) , I F S T , NT« AN , m TRAN , NTNX , OFRQ , E N I n , £N«D, 1 

TS 

12 

1FEN 

TS 

13 

DIMENSION  ICOFCJ7600).  COF(17600) 

TS 

1« 

EQUIVALENCE  <COF,EC>#  UCOF.EC) 

TS 

15 

N2=N»2 

TS 

16 

nop*2an2t l 

TS 

17 

IC12N2 

TS 

IS 

IC2SN0P 

TS 

19 

1CL=9999999 

TS 

20 

10LIMS9999999 

TS 

21 

LCQ*NCQMx/N2 

TS 

22 

NOCdsLCQ-NTMAX»l 

TS 

23 

IF  (NOCQ.LT.l)  GO  TO  1 

TS 

2a 

IC(JSLC3*N2 

TS 

25 

NOCQsNOCO«N2 

TS 

26 

L 1 a 1 

TS 

27 

L2SNOCO 

TS 

28 

NSFTslCQ-NOCQ 

TS 

29 

PRINT  36 

TS 

30 

GO  TO  2 

TS 

31 

l 

N0C02(LC3/2)*N2 

T5 

32 

IF  (NOcO.LT.N2)  CO  TO  35 

TS 

33 

LCOs2«NOCO 

TS 

3a 

UsNOCOAl 

TS 

35 

*■ 

L2=LC0 

TS 

36 

PRINT  37 

TS 

37 

2 

10CQ*0 

TS 

38 

nBASsO 

TS 

39 

icklg«o 

TS 

aO 

NELOsO 

TS 

a l 

NBOUT  »o 

TS 

U2 

PRINT  38 , NOCO.NTMAX 

TS 

as 

PRINT  19 

TS 

nu 

PRINT  ao 

7"3 

4$ 

DO  1 I o2  , n2  f 2 

TS 

Ub 

3 

CG(I)«0. 

TS 

a? 

TIMEb-DT 

TS 

ad 

c 

TS 

u9 

c 

BEGIN  TIME  LOOP 

TS 

50 

c 

TS 

51 

DO  31  ITIMtsl r nT  STEP 

TS 

52 

TIMEsTIMf aOT 

TS 

53 

ITLOCs(ITIME-l)«N2 

TS 

5a 

c 

TS 

65 

c 

t 

SET  UP  APPLIED  FIELD  AT  PRESENT  TIME 

TS 

56 

A- 59 

non 


c 

TS 

57 

CALL  EINCF  (TIMf) 

T 3 

58 

I F (ITI«t.E3,l)  GO  TO  16 

TS 

59 

I?  (ICFLG.EQ.O)  GO  TO  S 

TS 

60 

c 

TS 

bl 

c 

COMPUTE  FIELD  AT  PRESENT  TIME  DUE  TO  CURRENTS  AT  RETAMOEO  TIMES 

TS 

62 

c 

TS 

63 

n6kS= ( cnTmax-«)«N2»NELO*nOCQ-1 )/NOC0M 

TS 

64 

IF  (NBKS.GT .NdOuT)  NdKSsNBOUT 

TS 

6 5 

n0ASsNELO-N6kS«nOCQ»nOCO 

TS 

6b 

DO  4 Js 1 i NBKS 

TS 

67 

4 

BACKSPACE  11 

TS 

63 

RE*D  (11)  NIN, (C0(J) . Jsl.NOCQ) 

TS 

69 

NSINsl 

TS 

70 

5 

IF  (IQC.EQ.O)  GO  TO  7 

TS 

71 

REMIND  12 

TS 

72 

6 

READ  (12)  ICL.LENO, (COF ( J ) , J* 1 , ICL ) 

TS 

73 

IF  (ITIME.LE.LtND)  GO  TO  6 

TS 

74 

7 

ESCATeo, 

TS 

75 

ICF  * 1 

TS 

76 

1*1 

TS 

77 

s 

IND*=IcOK(ICF) 

TS 

78 

ICF*ICFM 

TS 

79 

IF  (INDX.LT.O)  GO  TO  9 

TS 

60 

ICFsJCFAlNDX/10000 

TS 

81 

GO  TO  6 

TS 

82 

9 

IF  ( ITJME.LE. (•INDX/JOOOO) ) GO  TO  8 

TS 

63 

GO  TU  11 

TS 

84 

JO 

ElNC(I)sEJNC(I)AESCAT 

TS 

85 

ESCATsO. 

TS 

66 

11 

IF  (1NDX.EO.O)  GO  TO  18 

TS 

87 

INDXS-INDX 

TS 

88 

K*INDX/10000 

TS 

89 

I*INOX-K*10000 

TS 

90 

K*ITIME-K 

TS 

91 

12 

ICGB=(K-1)»N2-NBAS 

TS 

92 

IF  (IC0B.LT.10LIM)  GO  TU  14 

TS 

93 

IF  (x8IN.EG.nUks)  GO  TO  14 

TS 

94 

READ  (11)  NIN,(C0{J),J*1,N0C<1) 

TS 

95 

N91N*N9lNt  l 

TS 

96 

NBASsNBAStNOCO 

TS 

97 

GO  TO  12 

TS 

98 

11 

IF  (1CF.LF.ICL)  GO  TO  14 

TS 

99 

READ  (12)  lCL.LtND» (COF (J),J*l, ICL) 

TS 

1 00 

I C F * 1 

TS 

1 0 1 

1« 

I NO  X : I COF (ICF) 

TS 

102 

ICFoICFM 

TS 

103 

IF  (Inox.lE.O)  GO  TO  10 

TS 

104 

JlsINWlOOOO 

TS 

105 

J2sl».0x-Jl»  10000 

TS 

106 

1C  J* ICUBt J2 

TS 

107 

IF  (K.EC.ITIME)  GO  TO  16 

TS 

106 

CO  IS  4*1. Jl 

TS 

109 

ESCAT«ktCAT»COF (ICF)«CO(ICQ) 

TS 

110 

1 CD* IC9+ 1 

TS 

111 

15 

ICFUCFM 

TS 

112 

GO  TO  |S 

TS 

113 

It 

ico*ica*i 

TS 

114 

ICfalCFM 

TS 

115 

DU  17  J*1 » J1 . 2 

TS 

116 

A- 60 


ESCATsFSCAMCOF  (ICF)»C0(ICG> 

TS 

117 

ICQ«IC0+2 

T 3 

1 1 8 

w 

ICF*ICF*2 

T 3 

1 19 

ICF a I C7  • l 

T 3 

120 

CO  TO  1 3 

TS 

121 

c 

TS 

122 

c 

SOLVE  MATRIX  EQUATION  F 08  CURRENTS 

AT  PRESENT 

TIME 

TS 

123 

c 

TS 

52a 

18 

I F ( I NT. F. 0.0)  GO  TO  19 

TS 

125 

CALL  L T SOL v ( CURF , C URF  , N , I x » E INC) 

T 3 

126 

GO  TO  20 

T3 

127 

10 

CALL  SOLVE  (N,CU«E.IP»EINC.N) 

T 5 

128 

20 

PRInT  al.  IT 1ME , TI«fc , Cfc INC ( 1 ) , IsJPl , JP2) 

T 3 

129 

PRINT  a},  ESRC 

TS 

1 10 

IF  (I 1 IMt .GT.nTRan)  GO  TO  21 

TS 

111 

ESURCUTlMtJtESRC 

T 3 

1 12 

21 

ICO=I tluc-i-nelo 

T 5 

111 

DO  22  131, N 

TS 

1 la 

ICO=ICQ *2 

TS 

135 

ICC83ICO+1 

TS 

116 

CQ(ICQ)sEINC(1) 

TS 

117 

22 

C0{IC08l=C0(IC08)tCI3«tINC(I) 

T 3 

118 

IF  (ITIME.EQ.NTSTFP)  GO  TO  31 

TS 

1 19 

IF  UCQ8.lT.LCQ)  GO  TO  26 

T 3 

1 a o 

IF  (10CQ.EU.1)  GO  TO  2J 

TS 

1 U 1 

lOCQsl. 

TS 

1 92 

RE«IN0  11 

TS 

1U3 

IF  (Ll.EQ.t ) GO  TO  23 

TS 

l aa 

wRITf  Cll)  NOC  0,  C C Q ( I ) , 1 3 1 , NOC  Q ) 

TS 

1 95 

N80UT  3N20UT 1 1 

7 S 

ia& 

1 CF  LG* 1 

TS 

1R7 

IOlIRsnocQ 

TS 

1 ae 

23 

wRJTE  (11)  NOCQ, (C0CI)*I=L1#L2) 

TS 

1 a9 

NBOUT  3N800T ♦ 1 

TS 

150 

NELOsNELOwNOCQ 

TS 

151 

IF  ( L 1 . E 3 . 1 ) GO  TO  2« 

TS 

152 

IC13N2-N0CU 

TS 

151 

lC23NaP-N0C0 

TS 

1 5a 

GO  TO  26 

TS 

155 

2« 

DO  26  I»1  , NSF  T 

TS 

156 

I NO  X r I ♦ NOC  Q 

TS 

157 

25 

C3U  JsCQUNOX) 

TS 

158 

NbASsNELO 

TS 

159 

26 

ICQsI TlOC-NELO»N2 

TS 

160 

C 

TS 

161 

C 

COMPOTE  CO*PONtNT  OF  Tint  INTEGRAL 

OF  CURRENT 

AT  NEXT  TIME  STEP 

TS 

162 

C 

THAT  DEPENDS  ON  PRESENT  and  P‘ST  VALUES  OF  CURRENT  And 

TS 

163 

C 

TIhe  Integral  of  current 

TS 

16a 

C 

TS 

165 

IF  (ITlME.tQ.l)  GO  TO  28 

T 3 

166 

DO  27  J*1 ,N 

TS 

167 

1C0*1C0»2 

TS 

168 

ICGxsicO-ICl 

TS 

169 

ICQVUCG- JC2 

TS 

170 

27 

C0UC0)sCQlICGX).CIl»CQ(IC0r)*0I2»ElNCCJ) 

TS 

171 

GO  Tu  10 

T» 

172 

28 

00  29  J31 ,N 

TS 

171 

ICGUCQW2 

TS 

17a 

lCO»*ICQ*ICi 

TS 

175 

20 

CQUCU)sC  i(ICQx)*0I2.£lNCtJ) 

T S 

176 

A-Gl 


10 

IF  (ICl.E0.N2J  IC2SMP 

T 3 

177 

IC1=N2 

TS 

178 

It 

CONTINUE 

T3 

1 79 

IF  (IOCU.EO.0)  CO  TO  14 

TS 

180 

12 

if  (ua.ce.icoa)  go  to  n 

TS 

181 

*91  TE  (tl)  NOCO,  ( C 0 ( I)  . 1=L1  ,L2) 

TS 

182 

NBUUT  =NBOUT ♦ 1 

TS 

1«I 

L 1 BL2* 1 

TS 

184 

L2=L2*N0C0 

TS 

185 

CO  TO  32 

T 3 

1 fib 

31 

Noca=icoa-ut*i 

TS 

1 fi  7 

"SITE  (li)  NCCC, (CC(I).I=L1.ICC3) 

TS 

188 

NSOUT^nBOUT* 1 

TS 

189 

RE  *i  I NO  It 

TS 

190 

10 

IF  (IOC.NE.O)  RE  * I NO  12 

TS 

191 

RETURN 

TS 

192 

15 

PRINT  42 . N 

TS 

191 

STOP 

TS 

194 

C 

TS 

195 

c 

TS 

196 

30 

FORMAT  ( 1 X . 49HS0LUT ION  WILL  USE  CURRENTS  FROM  CORE  STORAGE  ONL 

TS 

197 

IT) 

TS 

198 

37 

FORMAT  ( lX.aUMSOUUTION  RILL  USE  CURRENTS  FROM  FILE  STORAGE) 

TS 

199 

18 

FORMAT  (1X.13HLENGTH  GF  CURRENT  OUTPUT  BLOCK  I S , I 8 . 1 X , 6MW0RDS , , / , 

TS 

200 

tlX,U2HNUM0tR  OF  TIME  STEPS  USED  DURING  SULU T ION= , I 8 , // ) 

TS 

201 

39 

FORMAT  ( 44X.12H-  - - - CURRENT  SOLUTION  - - « -,/) 

TS 

202 

40 

FORMAT  (lX»4*iSTtP.8X,uMTlME,5x,02MSECMENT  CURRENTS  IN  Amps  (READ 

TS 

201 

1 ACROSS ) ....  -,/,2X,lHNO.,5X,6H(SEC.)»«X,b2HLAST  NUMBER  IN 

TS 

204 

2t  *CH  BLOCK  REPRfSENTS  SOURCE  STRENGTH) 

TS 

205 

01 

FORMAT  (/,lX,l4.E12.5,l0Ell.3,/,(l7x,10Elt.3)) 

TS 

206 

42 

FORMAT  ( 1 X, 72h«  “ERROR  - ROUTINE  TSOL--  INSUFFICIENT  STORAGE  FOR 

TS 

207 

1SECMENT  CURRENTS.  Ns.IIO) 

TS 

208 

0} 

FORMAT  (19X, lOMExCITATlON.Et 1.4) 

TS 

209 

END 

TS 

210- 
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SUBROUTINE  KjRt  <XM|,Yw|,Zm1,Xw2,Yw2,Zw2,RaD,NS,ITG) 

-I 

1 

wl 

2 

SUBROUTINE  M I RE  GENERATES  SEGMENT  GEOMETRY  DATA  FOR  A STRAIGHT 

-I 

J 

mI«E  (JF  ^3  segments. 

Wl 

(1 

WI 

5 

COMMON  /DATA/  N,np,X(200)#Y(200).Z(200),$I(200;,UI(200),alR(200), 

wl 

6 

BET (200) . I CONI (200) , I CON2 ( 2 0 C ) , I T A G C 2 0 0 } , IPX, IPY, IP* 

4 1 

7 

DIMENSION  x2(l),  Y2  ( 1 ) » Z2(l) 

MI 

6 

EQUIVALENCE  (X2(l),SI(m,  (Y2(l),»lP(m<  (Z2(l),  BET(tJ) 

*1 

9 

1ST  SNA  1 

Ml 

to 

NSNANS 

MI 

1 I 

NPaN 

Ml 

12 

IF  (NS.LT.I)  RETURN 

«I 

13 

FNSSNS 

Ml 

1« 

XOe(XM2-XWl)/FN3 

Ml 

15 

Y0e(YW2»Vwl )/FNS 

Ml 

16 

ZD=(Z*2-Zwl )/FNS 

Ml 

17 

x S 1 = xn i 

MI 

ia 

Y31=YW1 

MI 

19 

ZS1S2-J 

Ml 

20 

DO  1 I c I S T # N 

Ml 

21 

I T *G  C I ) a I TO 

Ml 

22 

X $2  * X S 1 ♦ XI) 

MI 

23 

Y32*YS1»YD 

MI 

2« 

Z32*  Z S l +ZO 

Ml 

25 

XdlsXSl 

MI 

26 

Y ( I ) e YS  1 

MI 

27 

Z(I)3ZS1 

Wl 

25 

X2(l)=X32 

Ml 

29 

Y2 ( I ) = YS2 

wl 

30 

Z2(l)=Z32 

MI 

31 

BKDsRAO 

Ml 

32 

XS1=XS2 

Mi 

33 

YSl SYS2 

-I 

34 

ZSUZS2 

Ml 

35 

RETURN 

Ml 

36 

END 

MI 

37- 
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APPENDIX  B 
DIAGNOSTICS 

1.  BLOCKING  Nj  N£  N3 
Routine:  FBLOCK 

This  is  printed  when  matrix  factorization  and  solution 
requires  file  storage.  It  does  not  indicate  an  error. 

Nj  = number  of  blocks  on  file  storage 

N^  = number  of  matrix  columns  per  block 

N^  = number  of  matrix  columns  in  last  block  (N^  s N^) 

2.  CHECK  DATA,  PARAMETER  SPECIFYING  SEGMENT 
POSITION  IN  A GROUP  OF  EQUAL  TAGS  MUST  NOT  BE 
ZERO 

Routine:  ISEGNO 

The  number  following  the  tag  number  reference  on  a data 
card  was  zero.  Execution  terminated. 

3 ERROR  FLAG  RETURNED  BY  SUBROUTINE  FORT, 

IFERR  = --- 
Routine:  ITOF 

This  indicates  a program  malfunction.  The  meaning  of 
IFERR  is  explained  in  the  comments  in  subroutine  FORT. 

4.  ‘ ERROR  - ROUTINE  RFLD  --  INCONSISTENT  SEGMENT 

CONNECTION  DATA 
Routine:  RFLD 

This  indicates  a program  malufnetion.  The  connection 
information  in  arrays  ICON1  and  ICON2  was  found  to  be 
contradictory.  A possible  cause  is  an  overflow  of  an 
array. 

5.  ERROR  - ROUTINE  TSOL  --  IT  FFICIENT  STORAGE 

FOR  SEGMENT  CURRENTS.  N = 

Routine:  TSOL 

The  dimension  of  array  CQ  must  be  at  least  4 x N where  N 
is  the  number  of  segments.  This  dimension  is  presently  6400 
for  the  GE  635/645  computer  and  5184  for  the  CDC  3800. 
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6. 


GEOMETRY  DATA  CARD  ERROR 
Routine:  DATAGN 

An  invalid  mnemonic  was  found  on  a card  where  a geometry 
data  card  was  expected. 

7.  INCORRECT  LABEL  FOR  A COMMENT  CARD 
Routine:  MAIN 

An  invalid  mnemonic  was  found  on  a card  where  a comment 
card  was  expected. 

8.  INPUT  ENERGY  NOT  COMPUTED  SINCE  SOURCE  HAS 
BEEN  TRANSFORMED 

Routine:  MAIN 

Request  for  energy  budget  must  preceed  any  requests 
involving  Fourier  transforms  since  the  Fourier  transform 
writes  over  the  record  of  the  source  time  function. 

9.  INTEGRATION  ATTEMPTED  OVER  SINGULARITY. 

B, C = --- 

Singularity  due  to  R going  to  zero  in  range  of  integration 
(R  = j s^  + Bs  + C ^ ).  Possible  cause  is  overlapping 
segments  in  geometry  specification. 

10.  INVALID  DATA  CARD  LABEL  AFTER  SOLUTION 
Routine:  MAIN 

An  invalid  mnemonic  was  found  on  a card  where  a data 
request  card  was  expected.  Valid  mnemonics  are  EB, 

PC,  AT,  RF,  NX,  and  EN 

11.  NO  SEGMENT  HAS  AN  ITAG  OF  -- 
Routine:  ISEGNO 

TAG  number  used  to  refer  to  a segment  docs  not  exist. 

12.  NUMBER  OF  EXCITATION  VALUES  EXCEEDS  ARRAY 
DIMENSION 

Number  of  voltage  sources  cannot  be  greater  than  the 
dimensions  of  EMAG  and  ISRC.  Dimension  is  now  10. 
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OVERFLOW  IN  FILLING  ARRAY  COF 
Routine:  COFS 

This  indicates  a program  malfunction.  May  possibly 
be  corrected  by  increasing  1ST  at  statement  CF  34A 
in  subroutine  COFS  to  cause  more  frequent  dumping  of 
data  to  file. 

OVERFLOW  IN  FILLING  ARRAY  SRT 
Routine  COFS 
Same  as  13. 

PIVOT  ( ) = --- 

Routine:  FACTR  or  LFACTR 

This  will  be  printed  during  Gauss  Doolittle  factoring  of  the 
interaction  matrix  when  a pivot  element  less  than  10  ^ is 
encountered,  and  indicates  that  the  matrix  is  nearly  singular. 
The  number  in  parenthesis  shows  on  which  pass  through  the 
matrix  the  condition  occurred.  This  is  usually  an  abnormal 
condition  although  execution  will  continue.  It  may  result 
from  coinciding  segments  or  a segment  of  zero  length. 
RETARDED  TIMES  FROM  SEGMENT  ---TO  SEGMENTS  --- 

AND---  DIFFER  BY TIME  STEPS 

Routine:  COFS 

Retarded  times  from  a point  to  two  segments  which  are 
connected  to  each  other  cannot  differ  by  more  than  one 
time  step.  This  is  usually  caused  by  a segment  with 
length  greater  than  the  time  increment  multiplied  by  the 
velocity  of  light. 

RETARDED  TIMES  TO  SEGMENTS  -- -AND -- -DIFFER 
BY- --TIME  STEPS 
Routine:  RFLD 


SEGMENT  ---  LIES  IN  OR  BEHIND  SYMMETRY  PLANE 
Routine:  CONECT 

All  segments  must  lie  on  the  positive  side  of  any  symmetry 
plane  s. 

SEGMENT  CONNECTION  ERROR.  J = --- 
Routine:  COFS 

Same  as  4.  Error  detected  at  segment  J. 

ZERO  TIME  STEPS  REQUESTED 
Routine:  MAIN 

Both  TMAX  and  NTSTEP  on  First  data  card  were  zero 
or  blank. 


Appendix  C 

List  of  Scientific  Subroutines  Used  in  TIMDOM 


ABS(X)  = absolute  value  of  X 

AIMAG(Z)  = imaginary  part  of  the  complex  number  Z, 

(result  is  real) 

AINT(X)  = integer  truncation  (result  is  real) 

ALOG(X)  = the  natural  log  of  X 

ALOGIO(X)  = the  log  to  the  base  ten  of  X 

ASIN(X)*  = arcsine  of  X,  result  in  radians 

ATAN2(Xj,  X^)  = arctangent  of  Hx^’  result  1°  radians  covering 

all  four  quadrants. 

2 

CEXP(Z)  = complex  exponential  (e  ) 

CMPLX  (Xj,  X2)  = formation  of  a complex  number  (Z  = Xj  + j X^) 
CONJG(Z)  = conjugate  of  the  complex  number  Z 

COS(X)  = cosine  of  X 

DABS  (X)  = absolute  value  of  a double  precision  number 

(result  is  double  precision) 

DATAN  (X)  = arctangent  of  a double  precision  number  (result 

is  double  precision) 

DLOG  (X)  = natural  log  of  a double  precision  number  (result 

is  double  precision) 

DSQRT  (X)  = square  root  of  a double  precision  number  (result 

is  double  precision) 

_ x 

EXP  (X)  = exponential  function  (e  ) 

FLOAT  (K)  = real  number  equivalent  of  integer  K 

IABS  (I)  = absolute  value  of  integer  I 

REAL(Z)  = real  part  of  the  complex  number  Z 

SIGN  (Xj,  X2)  = sign  of  X2  times  (X  [ 

SIN  (X)  = sine  of  X 

SORT  (X)  = square  root  of  X 

* Note,  this  routine  is  not  found  in  the  Honeywell/GE  scientific 
library;  therefore,  it  is  included  as  a subroutine  in  the  TIMDOM 
deck  for  this  machine. 
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METRIC  SYSTEM 


BASE  UNITS: 

Quantity 


Unit 


SI  Symbol  Formula 


length 

metre 

mass 

kilogram 

time 

second 

electric  currant 

ampere 

thermodynamic  temperature 

kelvin 

amount  of  substance 

mole 

luminous  intensity 

candela 

SUPPLEMENTARY  UNITS: 

plane  angle 

radian 

solid  angle 

steradian 

DERIVED  UNITS: 

Acceleration 

metre  per  second  squared 

activity  (of  a radioactive  source) 

disintegration  per  second 

angular  acceleration 

redian  per  second  squared 

angular  velocity 

radian  per  second 

area 

square  metre 

density 

kilogram  per  cubic  metre 

electric  capacitance 

farad 

electrical  conductance 

siemens 

electric  field  strength 

volt  per  metre 

electric  inductance 

henry 

electric  potential  difference 

volt 

electric  resistance 

ohm 

electromotive  force 

volt 

energy 

joule 

entropy 

joule  per  kelvin 

force 

newton 

frequency 

hertz 

illuminance 

lux 

luminance 

candela  per  square  metre 

luminous  flux 

lumen 

magnetic  field  strength 

ampere  per  metre 

magnetic  flux 

weber 

magnetic  flux  density 

tesla 

magnetomotive  force 

ampere 

power 

watt 

pressure 

pascal 

quantity  of  electricity 

coulomb 

quantity  of  heat 

joule 

radiant  intensity 

watt  per  steradian 

specific  heat 

joule  per  kilogram-kelvin 

stress 

pascal 

thermal  conductivity 

watt  per  metre-kelvin 

velocity 

metre  per  second 

viscosity,  dynamic 

pascal-second 

viscosity,  kinematic 

square  metre  per  second 

voltage 

volt 

volume 

cubic  metre 

wavenumber 

reciprocal  metre 

work 

joule 

SI  PREFIXES 

Multiplication  Factors 

i ooo  ooo  ooo  ooo  = in'1 

1 000  000  000  - 10* 

1 000  000  - 10* 

1 000  - 10' 
too  - in1 
to  » 10' 

0 1 - 10*' 
001  - 10-> 
0001  « 10-* 
0 000  001  > 10-* 
0 000  000  001  - 10*’ 
0 000  000  00(1001  • 10- •• 
o ooo  ooo  ooo  goo  ooi  * in-'* 
0 000  000  000  000  000  001  - 10-'* 

* To  ba  avoidad  whara  poaalbla 


m 

ki 

• 

A 

K 

mol 

cd 


rad 

sr 


m/s 

(disintegration  )/s 
rad/s 

rad/s 

m 

kg/m 

F 

A-s/V 

S 

A/V 

V/m 

H 

V-s/A 

V 

W/A 

V/A 

V 

W/A 

) 

N-m 

1* 

N 

kg-m/s 

Hz 

(cyclaj/s 

lx 

im/m 

ed/m 

Im 

cd-sr 

Aim 

Wb 

V-s 

T 

Wb/m 

A 

W 

i/» 

Pa 

N/m 

C 

1 

A-s 

N-m 

Wlsr 

|/kg-K 

Pa 

N/m 

W/m-K 

mis 

Pa-s 

ms 

V 

W/A 

m 

i 

(wavej/m 

N-m 

Prefix  SI  Symbol 


tnra 

(OH* 

maga 

kilo 

hecto* 

deka* 

dacl* 

rentl* 

mllll 

micro 

nano 

pIco 

inmto 

atto 


T 

(i 

M 

k 

h 

da 

d 

r: 
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